Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,349,220 (G)7→9 intergenic (‑157/+211) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase (NADH)/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,349,2131.G90.4% 226.1 / 20.1 83intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase (NADH)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/6);  new base G (17/58);  total (19/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 = 5.72e-01
*NC_0009131,349,2132.G81.9% 203.8 / 46.9 83intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase (NADH)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (5/10);  new base G (14/54);  total (19/64)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 115 total aligned reads displayed.

TAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCA  >  NC_000913/1349122‑1349292
                                                                                                 ||                                                                                        
tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g                                                                                         <  1:6056945/93‑1 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑g                                                                                      >  1:2415666/1‑93 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑g                                                                                      >  2:3387342/1‑93 (MQ=255)
 aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑g                                                                                      >  2:5388314/1‑93 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg                                                                                    <  1:4704316/93‑1 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg                                                                                    <  1:1221315/93‑1 (MQ=255)
   tAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg                                                                                    >  1:783208/1‑93 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑g                                                                                 <  1:4597534/93‑1 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑g                                                                                 <  2:3373784/93‑1 (MQ=255)
    aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑g                                                                                 >  1:1006415/1‑93 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                <  2:2121958/93‑1 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                <  2:1236626/93‑1 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                >  2:5626174/1‑93 (MQ=255)
     aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gg                                                                                >  2:5998091/1‑93 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  1:3449201/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  1:616245/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  1:1135277/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               <  2:5545302/93‑1 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               >  2:736367/1‑93 (MQ=255)
      aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCGGCG‑‑GCA‑‑g‑‑ggg‑‑ggg                                                                               >  1:1177168/1‑93 (MQ=255)
       cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gggg                                                                              >  1:4029118/1‑92 (MQ=255)
        aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggggg                                                                             >  1:4044465/1‑91 (MQ=255)
        aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑ggggg                                                                             >  1:5832396/1‑91 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑g‑‑ggg‑‑gggggg                                                                            >  2:5249082/1‑90 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggc                                                                          <  2:2429656/93‑1 (MQ=255)
         gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggc                                                                          >  2:642843/1‑93 (MQ=255)
           aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggcgg                                                                        >  1:29413/1‑91 (MQ=255)
             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggcgctt                                                                      >  2:237707/1‑89 (MQ=255)
             gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGg‑‑ggcgctt                                                                      >  2:4956103/1‑89 (MQ=255)
              aaCGGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCtttt                                                                     >  1:942791/1‑93 (MQ=255)
                aaGATAAACGGTTATTATAATCACCCTGGCTGTCAGTAGTTATAGTTGCCAGGTCGGACCGGAGCAGGATGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTc                                                                    <  2:2853589/92‑1 (MQ=255)
                  gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTccc                                                                  >  1:4095644/1‑93 (MQ=255)
                   aTAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTtccc                                                                 >  1:979596/1‑93 (MQ=255)
                    tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCt                                                                <  1:1352558/93‑1 (MQ=255)
                     aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCTc                                                               >  2:771180/1‑93 (MQ=255)
                      aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTccccccc                                                              >  1:2040629/1‑90 (MQ=255)
                       aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCTCAc                                                             >  2:1813759/1‑93 (MQ=255)
                        cGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCAcc                                                            >  2:756075/1‑93 (MQ=255)
                         ggTTATTATAGTCACCCTGGCTGTGAGTAGTTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCAccc                                                           <  1:4288771/93‑1 (MQ=255)
                         ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCAccc                                                           <  2:5645069/93‑1 (MQ=255)
                         ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCt                                                          >  1:4233108/1‑93 (MQ=255)
                          gTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCccccccccc                                                          >  1:5659498/1‑92 (MQ=255)
                           tttttaaAATCAACCTGGGCTGTAATAGCTAAAGTTCCCCGGTCCGGCCGGGACCGGCCTGGG‑‑GAGG‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTa                                                         <  2:3362634/90‑1 (MQ=255)
                           ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTa                                                         <  1:2971150/93‑1 (MQ=255)
                            tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCccccccccccc                                                        >  1:1538658/1‑90 (MQ=255)
                             attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAAc                                                     <  1:4472311/93‑1 (MQ=255)
                             attatAATCAACCTGACTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAAc                                                     <  2:5962353/93‑1 (MQ=255)
                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCct                                                  >  1:3177827/1‑93 (MQ=255)
                               tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCct                                                  <  2:3024143/93‑1 (MQ=255)
                                 taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCctc                                                 <  2:3142863/93‑1 (MQ=255)
                                   aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCtctc                                               <  1:4562490/93‑1 (MQ=39)
                                    tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTcc                                              <  2:4425908/93‑1 (MQ=39)
                                     cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTccc                                             >  1:2143565/1‑93 (MQ=39)
                                         cTGGCTGTGATTAGCTATAGTTCCCGGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                         <  2:2759138/93‑1 (MQ=37)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTTCCCCCCCCCCACCCCCTccccccc                                         >  1:3820194/1‑90 (MQ=255)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCAAACCCCTCCCCcga                                         >  2:3030242/1‑93 (MQ=255)
                                         cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                         <  2:2504697/93‑1 (MQ=38)
                                            ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg                                       <  1:3120308/93‑1 (MQ=38)
                                            ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgg                                       <  1:5911470/93‑1 (MQ=38)
                                             gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                      <  1:3390771/93‑1 (MQ=38)
                                              cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                     <  1:778632/93‑1 (MQ=37)
                                              cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                     <  2:3580948/93‑1 (MQ=37)
                                              cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                     <  1:5869653/93‑1 (MQ=37)
                                               tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCTAACCCTCCTCCCACCgggg                                     >  2:5825638/1‑93 (MQ=255)
                                               tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCCCACCCCCACCCCCTCCCCCCAggggg                                    >  2:2987704/1‑92 (MQ=255)
                                                gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCg                                   <  1:4909907/93‑1 (MQ=37)
                                                 tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                  <  2:5278214/93‑1 (MQ=37)
                                                 tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                  <  1:2081604/93‑1 (MQ=37)
                                                 tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                  <  2:3472245/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  1:6144376/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  1:1256274/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  1:1196156/93‑1 (MQ=37)
                                                  gAGTAGCTATAGTTGACAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                 <  2:2944012/93‑1 (MQ=25)
                                                   aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                <  1:4738844/93‑1 (MQ=37)
                                                   aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                <  1:1943569/93‑1 (MQ=37)
                                                    gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCGGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAggg                               <  2:3151092/93‑1 (MQ=25)
                                                     tggttataTTTGCCAGTCCCGACGGGGGCAGCCTGCGCCGGG‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgggg                              <  2:2480520/89‑1 (MQ=9)
                                                     tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgggg                              <  2:5475066/93‑1 (MQ=25)
                                                      aGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGa                             <  2:4054790/93‑1 (MQ=25)
                                                       gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCCAACCCCTCCCCCCCGgggggggggggg                              >  1:5031659/1‑93 (MQ=255)
                                                       gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAc                            <  2:768334/93‑1 (MQ=25)
                                                        gtataGTTGCCAGGGCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAcc                           <  2:5496529/92‑1 (MQ=21)
                                                        cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAcc                           <  1:1589656/93‑1 (MQ=25)
                                                          ataGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTTCCCCCCCCCCCACCCCCCCCCCCCAGAGggggggggggcgc                           >  1:2789774/1‑90 (MQ=255)
                                                            aGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAt                       <  1:2362574/93‑1 (MQ=25)
                                                             gTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                      <  2:818272/93‑1 (MQ=25)
                                                              ttGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtg                     <  1:3725250/93‑1 (MQ=25)
                                                               tGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgt                    <  2:326072/93‑1 (MQ=25)
                                                                      tCCGACCGGAGCAGGCTGCG‑‑GCAGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat            <  2:3909769/93‑1 (MQ=255)
                                                                      cgcGACCGGAGCAGGGTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCGCCCTAACCCCGCCCCCAGAGGGGCGAGGGGACCTGACTGGGCAAata             <  2:5790339/91‑1 (MQ=16)
                                                                       ccGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat            <  1:334463/93‑1 (MQ=25)
                                                                       ccGACCGGAGCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat            <  1:1620376/93‑1 (MQ=25)
                                                                         caCCGGAGCAGTCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTg          <  2:6308122/92‑1 (MQ=21)
                                                                           ccGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGGGTGGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt        <  1:342607/93‑1 (MQ=25)
                                                                               gAGCAGGCTGCG‑‑GCA‑‑GGAGGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcc     <  2:3517896/93‑1 (MQ=25)
                                                                                 gCAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcccc   <  2:84701/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTATCCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  2:3026622/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  2:5783191/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCAGGG‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  1:5772779/93‑1 (MQ=25)
                                                                                  cAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGGGCCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa  <  2:5436536/93‑1 (MQ=25)
                                                                                                 ||                                                                                        
TAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCG‑‑GCA‑‑G‑‑GGG‑‑GGG‑‑CGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCA  >  NC_000913/1349122‑1349292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: