Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,034 A→T 5.4% intergenic (‑43/+27) rnd ← / ← fadD ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,888,0340AT 5.4% 371.2 / 7.4 134intergenic (‑43/+27)rnd/fadDribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base A (65/60);  new base T (4/3);  total (70/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 139 total aligned reads displayed.

GTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGT  >  NC_000913/1887885‑1888183
                                                                                                                                                      |                                                                                                                                                     
gTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTaa                                                                                                                                                       >  2:1838532/1‑150 (MQ=255)
gTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTaa                                                                                                                                                       >  2:623291/1‑150 (MQ=255)
  gcgcGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAAcc                                                                                                                                                     >  2:1349277/1‑150 (MQ=255)
     cgTACGAACAAATTCAGTATCCACGGCTATCGCCGGAAACGCACGGACGGCTTCACACAAAGAACCCAGCGCATCGTCCGTGGTAAACATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGCCAAAAAAAACGCCGGATTAACCGGc                                                                                                                                                  <  2:748072/150‑1 (MQ=255)
      gTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAGACGCCGGATTAACCGGCg                                                                                                                                                 >  2:112352/1‑150 (MQ=255)
        aCGAACAAATTCAGTATCAAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGtt                                                                                                                                               <  2:352880/150‑2 (MQ=255)
           aaCAAATTCAGTATCCAGGGCTATCGCCGGGAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTga                                                                                                                                            >  1:790889/1‑150 (MQ=255)
           aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGAATAACCGGCGTCTga                                                                                                                                            >  2:536688/1‑150 (MQ=255)
             cAAATTCAGTATTCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTATgacg                                                                                                                                          >  2:1792612/1‑150 (MQ=255)
              aaaTTCAGTATCCACGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGTTAATCCGGCGtttttttt                                                                                                                                         <  1:688607/150‑8 (MQ=255)
                 ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                                                      >  2:1122339/1‑150 (MQ=255)
                   cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgac                                                                                                                                    >  2:1579360/1‑150 (MQ=255)
                    cgTATCCCGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgact                                                                                                                                   <  2:445446/149‑1 (MQ=255)
                    aGTATCCAGGTCTATCGACGGAAAGGCACGGACGCCTTCACACAAAGAAGCCATCGCTGCGTCGGTGGTAATCATTTGGTACTTCAAATTTTTTTCTCTTTAGTGGGCGACAAAAAAAACGCCGGATAACCCGGCGTCTGACgactgact                                                                                                                                   <  2:525930/150‑1 (MQ=255)
                    aGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTATTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGAGTAACCGGCGTCTGACgactgact                                                                                                                                   >  1:370814/1‑150 (MQ=255)
                    aGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACgactgact                                                                                                                                   >  1:1758645/1‑150 (MQ=255)
                        tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAc                                                                                                                               <  2:618664/150‑1 (MQ=255)
                        tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAc                                                                                                                               >  1:1098869/1‑150 (MQ=255)
                        accAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAc                                                                                                                               <  1:678600/149‑1 (MQ=255)
                         ccAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACg                                                                                                                              >  1:1430722/1‑150 (MQ=255)
                          cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGACGCCCGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGc                                                                                                                             >  2:1848579/1‑150 (MQ=255)
                          cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGCCGACTGAGTTAACGc                                                                                                                             >  2:64840/1‑150 (MQ=255)
                          cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGc                                                                                                                             <  2:827988/150‑1 (MQ=255)
                           agtgCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCt                                                                                                                            <  1:1029317/147‑1 (MQ=255)
                            gggCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTc                                                                                                                           <  1:224959/150‑1 (MQ=255)
                                   ctccGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                                                    <  2:416975/148‑1 (MQ=255)
                                   cGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                                                    >  2:265163/1‑150 (MQ=255)
                                        gAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                               >  1:1577224/1‑150 (MQ=255)
                                          aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGGCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                             <  2:902152/150‑1 (MQ=255)
                                           aGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGCCTGACTTAACGCTCAGGCTTTATTGTACa                                                                                                            >  2:1896391/1‑150 (MQ=255)
                                            ggCACGGCCGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCAc                                                                                                           >  2:1075303/1‑150 (MQ=255)
                                             gCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACt                                                                                                          >  1:1743196/1‑150 (MQ=255)
                                               acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACttt                                                                                                        >  1:266915/1‑150 (MQ=255)
                                                cggatggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAACCGCCGGATTAACCGGCGTCTGACGCCTGACTTAACGCTCAGGCGTTATTGTCCACTTTGc                                                                                                      >  2:376047/1‑150 (MQ=255)
                                                    cggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                                   <  2:1751516/150‑1 (MQ=255)
                                                      gCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgc                                                                                                 >  1:1944649/1‑150 (MQ=255)
                                                      gCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgc                                                                                                 <  2:1006844/150‑1 (MQ=255)
                                                         tCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACAGGCGTCTGACGACTGACTTAACGCGCAGGCTTTATTGTCCACTTTGACGCGCGCt                                                                                              >  2:1912492/1‑150 (MQ=255)
                                                           acacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                            >  1:1979286/1‑150 (MQ=255)
                                                            cacaAAGAAGCCAGCGCATCGTCTGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                            >  2:511769/1‑150 (MQ=255)
                                                            cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTATTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGCAATAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                           >  2:1236219/1‑150 (MQ=255)
                                                            cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGCTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                           >  2:638933/1‑150 (MQ=255)
                                                              caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                         >  2:1163080/1‑150 (MQ=255)
                                                              caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                         >  1:230644/1‑150 (MQ=255)
                                                                 agaagCCAGCGCATCGTCCGTGGTAATCATTTGGAAATTCAATTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                      <  2:1895421/150‑1 (MQ=255)
                                                                  gaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGt                                                                                     >  1:646899/1‑150 (MQ=255)
                                                                        aGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtct                                                                               >  2:1872572/1‑150 (MQ=255)
                                                                           gcATCGTGCGTGGTAATTATTTGGTTATTCAAATTGTTTTCTCTTTAGTGGGGGTCAATTAAAACGCCGTTTTAACCGTCGTCTGAGGACTGACTTAACGCTCAGGCTTTATTGTCCAGTTTGCCGCGCGCTTCGTCACGTAATtctctt                                                                            <  2:714920/150‑3 (MQ=255)
                                                                            cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtc                                                                           <  1:1147163/150‑1 (MQ=255)
                                                                              tCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGc                                                                         <  2:902676/150‑1 (MQ=255)
                                                                                      ggTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGCCGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGTGCGCTTCGTCACGTAATTCTCGTCGCAAAAtttt                                                                 >  1:1100942/1‑150 (MQ=255)
                                                                                       gTAATCTTTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                                <  1:1023810/150‑1 (MQ=255)
                                                                                       gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGTGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                                <  1:1664862/150‑1 (MQ=255)
                                                                                       gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                                <  1:1427593/150‑1 (MQ=255)
                                                                                         aaTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                                                              <  2:1403015/150‑1 (MQ=255)
                                                                                             aTTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                          >  1:423190/1‑150 (MQ=255)
                                                                                             aTTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACAGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCAGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                          >  1:122233/1‑150 (MQ=255)
                                                                                                  gTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGa                                                     <  2:1979286/150‑1 (MQ=255)
                                                                                                   tAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGGAAAAATTTTGCGACgagaatt                                                    <  2:1477984/150‑7 (MQ=255)
                                                                                                      ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTc                                                 <  1:1967550/150‑1 (MQ=255)
                                                                                                       tCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCg                                                <  2:1693387/150‑1 (MQ=255)
                                                                                                        ctaaTTGTTTTCTCTTTAGTGGGCGTCAAAATAAACGCCGGATTAACCGGCGTCTGACGAATGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGACACCTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                               <  1:1199591/148‑1 (MQ=255)
                                                                                                        cAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                               <  1:1602656/150‑1 (MQ=255)
                                                                                                        cAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGCCTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                              >  2:848945/1‑150 (MQ=255)
                                                                                                        cAAATTGTTTTCACTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                               <  1:1851712/150‑1 (MQ=255)
                                                                                                        aaaaTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCgg                                               <  2:1838255/149‑1 (MQ=255)
                                                                                                         aaaTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAGCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTGTCGTCGCAAAATTTTTGCGACGTTAGATTTCGGt                                              >  2:172295/1‑150 (MQ=255)
                                                                                                         aaaTTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGAATGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAAGTTAGATTTCGGt                                              >  1:1002810/1‑150 (MQ=255)
                                                                                                            ttGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCTCTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCTGTAAc                                           <  2:1735239/150‑1 (MQ=255)
                                                                                                              gTTTTCTCTTTTGTGGGCGTCAAAAAAAACGCCGGATTAACCGTCGACTGACGACTGACTTAAATCTCAGGCTTTATTGTCCACTTAGACGCGCGCTTCATCACGTAATTCTCGTCGCAAAAGTTTTACGACGTTAGATTTCGGTAAAtc                                         >  1:1216701/1‑150 (MQ=255)
                                                                                                              gTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                         >  2:1032026/1‑150 (MQ=255)
                                                                                                                 ttCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcata                                      >  2:874471/1‑149 (MQ=255)
                                                                                                                  tctcTTTAGTGGGCGTCAAAAATAACGCGGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatca                                     <  2:642282/150‑1 (MQ=255)
                                                                                                                  tctcTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatca                                     <  1:1349277/150‑1 (MQ=255)
                                                                                                                   ctctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                    <  1:1665356/150‑1 (MQ=255)
                                                                                                                   ctctTTAGTGGGCGGCAAAAAAAACGCCGGATTAAGCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                                    <  1:864247/150‑1 (MQ=255)
                                                                                                                      tttAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                                 <  1:1579360/150‑1 (MQ=255)
                                                                                                                       ttAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAATTTTTTCCGACGTTAGATTTCGGTACCTCAGCACGaaa                                >  1:1953613/1‑150 (MQ=255)
                                                                                                                           tGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTcc                            >  1:1223426/1‑150 (MQ=255)
                                                                                                                              gCGTCAAAAAATACGCCGGAGTAACCGTCGTCTGACGTCTGACTTAACGCTCAGGCTTTATTTTGCACTTTGCCGCGCGCTTCGTCACGTAAGTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAAGTCATCACGATACTccacc                         <  1:314495/150‑1 (MQ=255)
                                                                                                                               cGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacca                        <  2:214219/150‑1 (MQ=255)
                                                                                                                                 tCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGc                      >  1:470875/1‑150 (MQ=255)
                                                                                                                                     cgtcagaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                   >  1:1723489/7‑150 (MQ=255)
                                                                                                                                      gtcagaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg                  >  2:1658824/6‑150 (MQ=255)
                                                                                                                                        aaaaCGCCGGATTAACCGGCGTCTTACGTCTGACTTAACGCTCAGTCTTTATTGTCCACTTTGGCGCGCGCATATTCACGTAATTCTCGTCGCAAAATTTTTCCGAGGTTAGATTTCGGTAACTCATCACGAAACTCCAGCAGCTTCGGt                <  2:160674/150‑1 (MQ=255)
                                                                                                                                         aaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa               >  1:886645/1‑150 (MQ=255)
                                                                                                                                          gaCGCCGGTTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTAc              >  1:1518179/2‑150 (MQ=255)
                                                                                                                                            cGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACtt            >  2:188190/1‑150 (MQ=255)
                                                                                                                                             gCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACttt           <  1:1354534/150‑1 (MQ=255)
                                                                                                                                              ccGGATTAACCGGCGTCTGCCGCCTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGAATAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACGTTg          >  2:1637246/1‑150 (MQ=255)
                                                                                                                                              ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCCCGAAACTCCACCAGCTTCGGTACTTTg          >  1:1184412/1‑150 (MQ=255)
                                                                                                                                              ccGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTg          <  2:370814/150‑1 (MQ=255)
                                                                                                                                               cGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCCCGAAACTCCACCAGCTTCGGTACTTTGt         >  1:205977/1‑150 (MQ=255)
                                                                                                                                               cGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGt         >  1:1871391/1‑150 (MQ=255)
                                                                                                                                                ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTa        <  2:1430722/150‑1 (MQ=255)
                                                                                                                                                ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTa        >  1:521005/1‑150 (MQ=255)
                                                                                                                                                   ttAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATcc     <  1:638933/150‑1 (MQ=255)
                                                                                                                                                     aaCCGGCGTATGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCg   <  2:1950218/150‑1 (MQ=255)
                                                                                                                                                      aCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGt  >  2:237975/1‑150 (MQ=255)
                                                                                                                                                      aCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGt  >  1:193040/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGT  >  NC_000913/1887885‑1888183

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