Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,828,535 T→A 5.7% intergenic (+19/‑86) srlD → / → srlM sorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,828,5350TA 5.7% 801.5 / 24.6 224intergenic (+19/‑86)srlD/srlMsorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (115/95);  new base A (2/11);  total (118/106)
Fisher's exact test for biased strand distribution p-value = 8.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

ACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCT  >  NC_000913/2828386‑2828684
                                                                                                                                                     |                                                                                                                                                     
aCAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCAtt                                                                                                                                                       <  1:652133/150‑1 (MQ=255)
aCAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCAtt                                                                                                                                                       <  1:4269190/150‑1 (MQ=255)
  aaaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGAGGAGATCCATTaa                                                                                                                                                     >  2:2607974/1‑150 (MQ=255)
   aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAg                                                                                                                                                    <  2:3652216/150‑1 (MQ=255)
   aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAg                                                                                                                                                    <  1:1862366/150‑1 (MQ=255)
      gaCCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGAt                                                                                                                                                 <  1:4150264/149‑1 (MQ=255)
        ccGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGAtct                                                                                                                                               <  1:2860536/150‑1 (MQ=255)
        ccGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGAtct                                                                                                                                               <  2:1056913/150‑1 (MQ=255)
         cGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGAtctc                                                                                                                                              <  2:903217/150‑1 (MQ=255)
         cGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGAtctc                                                                                                                                              <  2:1348540/150‑1 (MQ=255)
         cGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGAtctc                                                                                                                                              >  1:1594376/1‑150 (MQ=255)
         cGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGAtctc                                                                                                                                              >  1:4037882/1‑150 (MQ=255)
           cTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCg                                                                                                                                            <  1:2190111/150‑1 (MQ=255)
            tCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGt                                                                                                                                           <  1:2505736/150‑1 (MQ=255)
              acaCGCGGCTGCGATTCTCCAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGATCTCCgctg                                                                                                                                         <  2:3533424/148‑4 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCTACAGCGGAGATCCTTAATGGATCTCCgctg                                                                                                                                         <  2:3045099/150‑4 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGATCTCCgctg                                                                                                                                         <  2:267126/150‑4 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGATCTCCgctg                                                                                                                                         <  2:591061/150‑4 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTga                                                                                                                                         <  1:950011/150‑1 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTga                                                                                                                                         >  1:2725348/1‑150 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTga                                                                                                                                         <  2:1910165/150‑1 (MQ=255)
              aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTga                                                                                                                                         <  2:476568/150‑1 (MQ=255)
                aCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGCTCTCCGTgaga                                                                                                                                       <  1:3339156/150‑1 (MQ=255)
                aCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTgaga                                                                                                                                       <  1:1901575/150‑1 (MQ=255)
                    ggCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACtat                                                                                                                                   <  2:3374745/150‑1 (MQ=255)
                         cGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAAt                                                                                                                              <  2:2856374/150‑1 (MQ=255)
                          gATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATg                                                                                                                             >  2:1417683/1‑150 (MQ=255)
                              aTCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGTCTg                                                                                                                         >  1:3757232/1‑150 (MQ=255)
                               tCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGa                                                                                                                        <  1:2252459/150‑1 (MQ=255)
                                cAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGAt                                                                                                                       >  2:3487835/1‑150 (MQ=255)
                                cAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGAt                                                                                                                       >  1:3533752/1‑150 (MQ=255)
                                  aGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATgc                                                                                                                     >  2:1935391/1‑150 (MQ=255)
                                    aTGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATgcgc                                                                                                                   <  1:14054/150‑1 (MQ=255)
                                    aTGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATgcgc                                                                                                                   <  1:3981973/150‑1 (MQ=255)
                                       tgCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTAc                                                                                                                >  1:708355/1‑150 (MQ=255)
                                         cTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATgagcgtagc                                                                                                              <  1:2024411/150‑6 (MQ=255)
                                                 gctgctGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCTATGTCACCGGCGGTCAGTTGATGTTCTGATCAACTGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATTCGCTACGCTCATCAgg                                                                                                      <  2:4260600/150‑1 (MQ=255)
                                                 gctgctGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAgg                                                                                                      <  1:2938115/150‑1 (MQ=255)
                                                 gctgctGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAgg                                                                                                      <  1:4170539/150‑1 (MQ=255)
                                                  ctgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                     <  1:3176204/150‑1 (MQ=255)
                                                        ttCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAg                                                                                               <  1:2422549/150‑1 (MQ=255)
                                                          cTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGa                                                                                             <  1:1217992/150‑1 (MQ=255)
                                                           tACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGAc                                                                                            >  2:2216305/1‑150 (MQ=255)
                                                                   tCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCa                                                                                    <  1:786347/150‑1 (MQ=255)
                                                                    ccTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCTCTACGCTCATCAGGCATACAGGACTTCCGCCAc                                                                                   >  2:602896/1‑150 (MQ=255)
                                                                      tatggCGTCGTACTGCACCGGACAGTCGATCTATGTCACCGGCGGTGAGGTGATGTTCTGATAAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                 <  2:3871410/147‑1 (MQ=255)
                                                                      tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                 >  2:1201361/1‑150 (MQ=255)
                                                                      tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                 >  2:2937332/1‑150 (MQ=255)
                                                                      tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                 <  2:3799975/150‑1 (MQ=255)
                                                                      tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                 >  2:1201291/1‑150 (MQ=255)
                                                                      tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                 >  2:811844/1‑150 (MQ=255)
                                                                       aaGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTAc                                                                                <  1:1173189/150‑1 (MQ=255)
                                                                         ggCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACAt                                                                              <  1:2041725/150‑1 (MQ=255)
                                                                          gCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACAtt                                                                             >  1:1806856/1‑150 (MQ=255)
                                                                             tcgtACTGCACCGGACAGTCGATCAATGGCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAg                                                                          <  2:1222625/150‑1 (MQ=255)
                                                                               gtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGa                                                                        >  2:211200/1‑150 (MQ=255)
                                                                                  cTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                     <  2:434911/150‑1 (MQ=255)
                                                                                   tGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAg                                                                    >  2:518869/1‑150 (MQ=255)
                                                                                    gCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGt                                                                   >  1:161133/1‑150 (MQ=255)
                                                                                     cACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGtt                                                                  >  2:3909673/1‑150 (MQ=255)
                                                                                       ccGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAt                                                                >  2:3186891/1‑150 (MQ=255)
                                                                                        cGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATg                                                               >  1:3041707/1‑150 (MQ=255)
                                                                                         ggACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATgg                                                              <  2:981664/150‑1 (MQ=255)
                                                                                             aGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATc                                                          <  1:4085623/150‑1 (MQ=255)
                                                                                              gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATcc                                                         >  2:1635103/1‑150 (MQ=255)
                                                                                               tCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCg                                                        >  1:2811571/1‑150 (MQ=255)
                                                                                               tCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCg                                                        >  1:2811246/1‑150 (MQ=255)
                                                                                                     aaTGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtca                                                  >  2:3306782/1‑150 (MQ=255)
                                                                                                         tCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCAc                                              >  1:1023570/1‑150 (MQ=255)
                                                                                                            ccGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGt                                           <  1:3701441/150‑1 (MQ=255)
                                                                                                               gcggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGc                                        <  2:2309782/150‑1 (MQ=255)
                                                                                                                  gTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGt                                     >  1:1460578/1‑150 (MQ=255)
                                                                                                                     aGGTGATGTTCTGCTCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGCATGCCTGAGGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACCTTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTAt                                  <  1:2903102/150‑1 (MQ=255)
                                                                                                                     aGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTAt                                  <  1:1578597/150‑1 (MQ=255)
                                                                                                                     aGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTAt                                  <  1:3265802/150‑1 (MQ=255)
                                                                                                                       gTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCg                                >  2:2759538/1‑150 (MQ=255)
                                                                                                                       gTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCg                                >  1:2236900/1‑150 (MQ=255)
                                                                                                                       gTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCg                                >  2:1517514/1‑150 (MQ=255)
                                                                                                                        tGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGc                               >  2:2343604/1‑150 (MQ=255)
                                                                                                                         gATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGcc                              <  2:771692/150‑1 (MQ=255)
                                                                                                                         gATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGcc                              <  2:771523/150‑1 (MQ=255)
                                                                                                                         gATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGcc                              <  1:1683831/150‑1 (MQ=255)
                                                                                                                          aTGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                             <  1:2306169/150‑1 (MQ=255)
                                                                                                                          aTGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                             <  1:2390063/150‑1 (MQ=255)
                                                                                                                          aTGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                             <  2:1689196/150‑1 (MQ=255)
                                                                                                                           tGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTg                            <  2:2630021/150‑1 (MQ=255)
                                                                                                                           tGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTg                            <  2:2305784/150‑1 (MQ=255)
                                                                                                                             ttCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgt                          >  1:3994763/1‑150 (MQ=255)
                                                                                                                             ttCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgt                          >  2:2554051/1‑150 (MQ=255)
                                                                                                                                 gATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgc                      >  1:3605782/1‑150 (MQ=255)
                                                                                                                                  aTCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgcg                     >  1:3888378/1‑150 (MQ=255)
                                                                                                                                  aTCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgcg                     >  1:4184372/1‑150 (MQ=255)
                                                                                                                                  aTCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgcg                     >  1:1126588/1‑150 (MQ=255)
                                                                                                                                   tCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgcgc                    >  1:2658672/1‑150 (MQ=255)
                                                                                                                                   tCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCAAACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgcgc                    >  1:1940095/1‑150 (MQ=255)
                                                                                                                                    cAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCa                   >  2:422440/1‑150 (MQ=255)
                                                                                                                                         gCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTgg              >  1:881416/1‑150 (MQ=255)
                                                                                                                                               aTCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAg        <  2:1754962/150‑1 (MQ=255)
                                                                                                                                                tCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAgg       >  2:840022/1‑150 (MQ=255)
                                                                                                                                                    ttAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcggc   <  1:3404836/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCT  >  NC_000913/2828386‑2828684

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