Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,148,957 C→T 5.1% intergenic (+20/+20) gpr → / ← yqhA L‑glyceraldehyde 3‑phosphate reductase/UPF0114 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,148,9570CT 5.1% 254.6 / 2.2 78intergenic (+20/+20)gpr/yqhAL‑glyceraldehyde 3‑phosphate reductase/UPF0114 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (44/30);  new base T (2/2);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCG  >  NC_000913/3148813‑3149105
                                                                                                                                                |                                                                                                                                                    
aGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAAc                                                                                                                                                  <  4:281170/149‑1 (MQ=255)
 gCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAAcc                                                                                                                                                 >  1:182626/1‑149 (MQ=255)
      tAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGcc                                                                                                                                            >  7:191190/1‑149 (MQ=255)
        gaggagAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCGGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGATAATAACGGGCGCGCGAACCGCGccc                                                                                                                                           >  5:174760/1‑148 (MQ=255)
        gaggagAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGccc                                                                                                                                           >  1:139572/1‑148 (MQ=255)
        gaggagAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCg                                                                                                                                          <  8:154583/149‑1 (MQ=255)
         aggagAACGTGCAGGCGCTGAATAATCTTACATTTAGCACCAAGGAGCTGGCGCAGATTGAGCAGCAGATCGCCGATGGCGTGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGGGCCCGa                                                                                                                                         <  6:136417/149‑1 (MQ=255)
          ggagAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGaa                                                                                                                                        >  3:20723/1‑149 (MQ=255)
           gagaACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGTGAACCGCGCCCGAat                                                                                                                                       >  8:316570/1‑149 (MQ=255)
             gaACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAat                                                                                                                                       >  1:163356/1‑147 (MQ=255)
             gaACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAat                                                                                                                                       >  2:55713/1‑147 (MQ=255)
                cGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatcagat                                                                                                                                  <  2:163356/149‑1 (MQ=255)
                 gTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatcagatc                                                                                                                                 <  6:286385/149‑1 (MQ=255)
                    cAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatcagatcag                                                                                                                               <  3:56715/148‑1 (MQ=255)
                             aataatCTGAGATTTAGCACCAAGGAGCTGGCGCAGATTTATCAGCATATCGCCGATGGCGAGCTGAATCTGGGGCAGGCGGCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGa                                                                                                                      <  8:259288/148‑1 (MQ=255)
                             aataatCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGTCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGAc                                                                                                                     <  3:49435/149‑1 (MQ=255)
                              ataatCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCTAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGAc                                                                                                                     <  3:44924/148‑1 (MQ=255)
                                 atCTGACATTTAGCACCAAGGAGCTGGCGCAGTTTGATGAGCATATCGCCGATGGCGAGCTGAATGTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAAGGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGagt                                                                                                                 <  1:141710/149‑1 (MQ=255)
                                 atCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGagt                                                                                                                 >  2:271267/1‑149 (MQ=255)
                                  tCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGagt                                                                                                                 >  2:164540/1‑148 (MQ=255)
                                  tCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGagt                                                                                                                 >  8:350995/1‑148 (MQ=255)
                                    tGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGagtcag                                                                                                              >  1:319649/1‑149 (MQ=255)
                                        aTTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCgg                                                                                                          <  1:271267/149‑1 (MQ=255)
                                         tttAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGt                                                                                                         >  5:6960/1‑149 (MQ=255)
                                         tttAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGt                                                                                                         >  8:20493/1‑149 (MQ=255)
                                         tttAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGt                                                                                                         <  1:55713/149‑1 (MQ=255)
                                         tttAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGt                                                                                                         <  6:254780/149‑1 (MQ=255)
                                              cACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGACCCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTc                                                                                                        >  5:240914/1‑145 (MQ=255)
                                              cACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGa                                                                                                    <  1:102082/149‑1 (MQ=255)
                                               aCCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGAt                                                                                                   <  2:139572/149‑1 (MQ=255)
                                                ccAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGAtt                                                                                                                          >  1:175274/1‑125 (MQ=255)
                                                ccAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGAtt                                                                                                                          <  2:175274/125‑1 (MQ=255)
                                                ccAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGa                                                                                                    >  5:109733/1‑147 (MQ=255)
                                                 cAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGTCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAg                                                                                                 >  8:151866/1‑149 (MQ=255)
                                                  aaGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGc                                                                                                <  1:51698/149‑1 (MQ=255)
                                                     gAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCa                                                                                             >  5:337585/1‑149 (MQ=255)
                                                          ggCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGACCCGCGCCCGAATCAGATCACTGATTATGACGAGTCAGTCGGTCAAGCTAGCCCATcac                                                                                         >  2:299871/1‑148 (MQ=255)
                                                           gcgcAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGTGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATcac                                                                                         <  7:316570/147‑1 (MQ=255)
                                                            cgcAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACaaa                                                                                      <  8:191190/149‑1 (MQ=255)
                                                             gcAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATc                                                                                           <  1:185282/143‑1 (MQ=255)
                                                             gcAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATc                                                                                           >  2:185282/1‑143 (MQ=255)
                                                             gcAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAAt                                                                                     >  7:39163/1‑149 (MQ=255)
                                                             gcAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGGGGCAGGCGGCTTCCGATAATTGACCTGTTTATATCGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAAt                                                                                     <  6:174760/149‑1 (MQ=255)
                                                                     aTCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGTTCAGTGATTATGACGAGTGAGTCGGTCTAGATAGCCGCTCTCAAATGCAGaacg                                                                             >  5:156931/1‑147 (MQ=255)
                                                                      tCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGACCCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAg                                                                             >  6:73370/1‑148 (MQ=255)
                                                                            tatCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAAc                                                                      >  5:286979/1‑149 (MQ=255)
                                                                            tatCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAAc                                                                      <  2:182626/149‑1 (MQ=255)
                                                                              tCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGATCCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGt                                                                    >  1:51607/1‑149 (MQ=255)
                                                                              tCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAg                                                                                                     >  7:100083/1‑116 (MQ=255)
                                                                              tCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAg                                                                                                     <  8:100083/116‑1 (MQ=255)
                                                                                gCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCa                                                                  <  1:164540/149‑1 (MQ=255)
                                                                                 ccGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAg                                                                 <  7:151866/149‑1 (MQ=255)
                                                                                        cGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGAt                                                           >  6:158552/1‑148 (MQ=255)
                                                                                        cGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATa                                                          <  6:6960/149‑1 (MQ=255)
                                                                                         gAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATaa                                                         >  7:8889/1‑149 (MQ=255)
                                                                                              gAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACg                                                    >  3:85225/1‑149 (MQ=255)
                                                                                               aaTCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACg                                                    <  6:109733/148‑1 (MQ=255)
                                                                                                aTCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTa                                                  <  3:28487/149‑1 (MQ=255)
                                                                                                   tgtgGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACca                                               >  4:117976/1‑149 (MQ=255)
                                                                                                     tgGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACcaca                                             >  2:13749/1‑149 (MQ=255)
                                                                                                      ggCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACAt                                            >  8:325046/1‑149 (MQ=255)
                                                                                                       gCAGGCGTCTTCCGATAAATGACATGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACAt                                            >  4:172355/1‑148 (MQ=255)
                                                                                                           gCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGtt                                       >  6:65799/1‑149 (MQ=255)
                                                                                                             gTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTg                                     <  3:117976/149‑1 (MQ=255)
                                                                                                                   cGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATc                                <  2:319649/148‑1 (MQ=255)
                                                                                                                    gATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATcag                              >  5:111575/1‑147 (MQ=255)
                                                                                                                     aTAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCgg                             >  6:60469/1‑149 (MQ=255)
                                                                                                                         aTGACCTGTTAATAACGGGCGCGTGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACa                         >  8:311082/1‑149 (MQ=255)
                                                                                                                              cTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACAttttt                    >  5:69081/1‑149 (MQ=255)
                                                                                                                               tGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTc                   >  2:334645/1‑149 (MQ=255)
                                                                                                                                gTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCg                  >  2:124112/1‑149 (MQ=255)
                                                                                                                                gTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCg                  >  8:58501/1‑149 (MQ=255)
                                                                                                                                  taataaCGGGCGCGTGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGc                 <  7:311082/148‑1 (MQ=255)
                                                                                                                                   aataaCGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCAt               >  1:89165/1‑149 (MQ=255)
                                                                                                                                   aataaCGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAAGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCAt               >  6:248105/1‑149 (MQ=255)
                                                                                                                                     taaCGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATcc             <  7:20493/149‑1 (MQ=255)
                                                                                                                                      aaCGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTAAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTAACCGGGAACTTTTTCGCATCCa            >  3:66939/1‑149 (MQ=255)
                                                                                                                                         gggCGCGCGAACCGCGCCCGAATCAGATCGGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATa          >  5:314029/1‑148 (MQ=255)
                                                                                                                                            cgcgcgAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGa      <  5:65799/149‑1 (MQ=255)
                                                                                                                                             gcgcgAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGAc     <  6:337585/149‑1 (MQ=255)
                                                                                                                                                cgAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGAcgcg  <  5:60469/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCG  >  NC_000913/3148813‑3149105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑