Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,596,823 A→C 7.8% E278D (GAA→GAC lgoD → L‑galactonate oxidoreductase; L‑gulonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,596,8230AC 7.8% 1763.4 / 32.8 584E278D (GAA→GAClgoDL‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (380/144);  new base C (2/57);  total (383/201)
Fisher's exact test for biased strand distribution p-value = 1.01e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 604 total aligned reads displayed.

CGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGAT  >  NC_000913/4596675‑4596972
                                                                                                                                                    |                                                                                                                                                     
 gCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAac                                                                                                                                                     >  2:2088890/1‑150 (MQ=255)
  cTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCGATCCGGAATTCCATAAAAAAGAAacga                                                                                                                                                   >  1:5261266/1‑150 (MQ=255)
     gcgcAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAacgacg                                                                                                                                                 >  2:1575247/1‑150 (MQ=255)
       gcAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACgat                                                                                                                                               >  1:1011271/1‑150 (MQ=255)
       gcAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACgat                                                                                                                                               >  2:4393488/1‑150 (MQ=255)
         aGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACgatga                                                                                                                                             >  1:7508251/1‑150 (MQ=255)
         aGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGACACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACgatga                                                                                                                                             >  2:2532089/1‑150 (MQ=255)
            aaGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAATAGAAACGACGATGATgg                                                                                                                                          >  2:8147214/1‑150 (MQ=255)
             aGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATggg                                                                                                                                         >  1:3703778/1‑150 (MQ=255)
               tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                                       >  2:10201162/1‑150 (MQ=255)
                 aTCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGc                                                                                                                                     >  1:1631121/1‑150 (MQ=255)
                      cgcgACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCaa                                                                                                                                >  2:10044408/1‑150 (MQ=255)
                             ggTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacg                                                                                                                         >  2:4650772/1‑150 (MQ=255)
                              gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                                                                        >  1:2746503/1‑150 (MQ=255)
                              gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                                                                        >  2:3304326/1‑150 (MQ=255)
                              gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                                                                        >  1:45143/1‑150 (MQ=255)
                              gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                                                                        >  2:7669837/1‑150 (MQ=255)
                              gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGGATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAccgccacgc                                                                                                                        >  1:7710098/1‑150 (MQ=255)
                                 aTCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCgg                                                                                                                     >  1:3249527/1‑150 (MQ=255)
                                   cAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGaa                                                                                                                   >  1:7838868/1‑150 (MQ=255)
                                   cAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGaa                                                                                                                   >  1:2247043/1‑150 (MQ=255)
                                    aaCATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAg                                                                                                                  >  1:4825853/1‑150 (MQ=255)
                                                 tAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgt                                                                                                     >  1:6299062/1‑150 (MQ=255)
                                                  aaCACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtc                                                                                                    >  1:10570364/1‑150 (MQ=255)
                                                  aaCACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtc                                                                                                    >  1:1109679/1‑150 (MQ=255)
                                                   acacCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcg                                                                                                   >  1:1286666/1‑150 (MQ=255)
                                                    cacCGTTATTTTGATTCGTCACGGCGGCACGGTGGTTTTTGTCGGCATTTTTAACGTTGAGTTGCATTTCTCCGATCCGGAATTCAAAAAAAAAAAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcgg                                                                                                  <  2:1299542/150‑1 (MQ=255)
                                                      cctttattttGTTTCGTCCGGGGGGCAGGGTTTTTTTTGGGGGCCTTTTTATTGGTTCGTTGCTGTTCCCCCACCCTTAATACAATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcggtc                                                                                                <  1:6609532/143‑1 (MQ=255)
                                                       cGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcggtcg                                                                                               >  1:7891661/1‑150 (MQ=255)
                                                       cGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcggtcg                                                                                               >  1:2463249/1‑150 (MQ=255)
                                                        gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                                              >  1:2357324/1‑150 (MQ=255)
                                                        gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                                              >  1:6448530/1‑150 (MQ=255)
                                                         tGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGAc                                                                                             >  2:5117681/1‑150 (MQ=255)
                                                         tGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGAc                                                                                             >  2:8441661/1‑150 (MQ=255)
                                                                 tttCGGCCCGGCGGCTGGGTTTTTTTTGTCGGCGTTTTTACTGGTGTTTTGCAGTTCCCCGATCCGTAATCCAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCg                                                                                     <  1:1448651/149‑1 (MQ=255)
                                                                   tCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCgga                                                                                   >  2:6518016/1‑150 (MQ=255)
                                                                   tCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCgga                                                                                   >  2:634608/1‑150 (MQ=255)
                                                                    cttCAGGGCGGCCGGGTGTTTTTGGGCGGCCTTTTTAATGTTTCGTGGCCTTTCTCCGATCGGTAATACAAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaa                                                                                  <  1:8840113/148‑1 (MQ=255)
                                                                    cgacACGGCGGCCCGGTGGTTTTTGTCGGCCTTTTTAAAGGTGTGTTGCTGTTCTCCGAGCCGTAATCCCAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaa                                                                                  <  2:9429970/147‑1 (MQ=255)
                                                                      tCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaagg                                                                                >  2:11027336/1‑150 (MQ=255)
                                                                      tCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaagg                                                                                >  1:5334405/1‑150 (MQ=255)
                                                                      tCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGACGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaagg                                                                                >  2:7943716/1‑150 (MQ=255)
                                                                        aCGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaaggaa                                                                              >  2:588569/1‑150 (MQ=255)
                                                                          ggcggcAGGTTGTTTTTTTTCGGCCTTTTTAATGGTGTGTTGCAGTTCTCCTATCCGTATTTCAATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGaaaa                                                                            <  2:9764184/150‑1 (MQ=255)
                                                                           ggggctgggtggtTTTTGTGGGCCTTTTTAATGGTGTGTTGCTGTTCCCCGATCCGTATTTCAATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGaaaaa                                                                           <  1:8139708/143‑1 (MQ=255)
                                                                            cggccgggtggaATTGGCCGGCCGGTTTAATGGTGTGTGGCCGGTCACCGAGCCGGAATTCCAGAACACATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAAt                                                                          <  2:4000610/144‑1 (MQ=255)
                                                                              gccgggtgttTTTTGTCGGCATTTTTAACGGTGTGTGGCTTTACTCCGATCCGGAATCCAATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                                                                        <  1:3926259/146‑1 (MQ=255)
                                                                               ctgggtggtTTTTGTCGGCGTTTTTAATGTTGTGTTGATTTACCCCCAACGGTAATACAATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCAc                                                                       <  1:5648485/147‑1 (MQ=255)
                                                                                cgggtgttTTTTGTCGGCCTGTTTAAAGGTGTTTTGCTGTTCTCCGATCCGTAATACAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCAct                                                                      <  1:6805668/148‑1 (MQ=255)
                                                                                 cGGTGGTTTTTGTGGGCCTGTTTAACGGTGAGTTGCTGTTCTCCGATCCGGACTTCGATAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActg                                                                     <  2:9205754/150‑1 (MQ=255)
                                                                                  ggtggtTTTTGTCGGCCTTTTTAATGGTGTGTTGCCTTTCTCCGATCCGTAATTCAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgc                                                                    <  1:3465566/150‑1 (MQ=255)
                                                                                   gtggtATTTGTCGGCCTGTTTAATGGTGAGTTGCTGTTCCCCGACCCGTAATTCCAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgct                                                                   <  1:10966630/150‑1 (MQ=255)
                                                                                    tgttttttGGCGGCCTTTTTAACGGTGAGTGGCTGTTCTCCGATCCGGATTTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgctg                                                                  <  2:5738174/145‑1 (MQ=255)
                                                                                     ggtTTTTGTCGGCCTGTTTAACGGTGTGTTGCAGTCCTCCGATCCGTAATACCTTAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACGGATGGCGGAAGGAAAAATCACTGCTGa                                                                 <  1:7053370/150‑1 (MQ=255)
                                                                                      ttttttGCCGGCGTGTTTATTGTTTTGTTGCATTTCCCCGATCCGTAATCCAAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGAc                                                                <  1:1903386/147‑1 (MQ=255)
                                                                                         tttGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatg                                                             >  2:2172437/1‑150 (MQ=255)
                                                                                         tttGGCGGCTTTTTTAATGTTTTGTTGCTGTTCTCCGATCCGTAATTCAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatg                                                             <  2:2849997/150‑1 (MQ=255)
                                                                                          ttGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                                            >  1:5943017/1‑150 (MQ=255)
                                                                                            gTCGGCCTTTTTAACTGTGTTTTTCTTTTCACCGATCCTCAATACCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatgatg                                                          <  1:2585505/150‑1 (MQ=255)
                                                                                                gCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCACCGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTaa                                                      >  2:10470849/1‑150 (MQ=255)
                                                                                                gCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTaa                                                      >  2:5578346/1‑150 (MQ=255)
                                                                                                gCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACAGGATGTTaa                                                      >  1:7266387/1‑150 (MQ=255)
                                                                                                   tGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAccc                                                   >  2:3893501/1‑150 (MQ=255)
                                                                                                   tGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAccc                                                   >  2:10340350/1‑150 (MQ=255)
                                                                                                    gTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCa                                                  >  2:1176027/1‑150 (MQ=255)
                                                                                                      ttAATGGTGTGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc                                                <  2:3644466/150‑1 (MQ=255)
                                                                                                      ttAAAGGTGAGTTGCTGTTCTCCGATCCGTATTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc                                                <  2:7306459/150‑1 (MQ=255)
                                                                                                       tAAAGGTTAGTTGCAGTTCTCCGATCCGGATTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCg                                               <  2:2751223/150‑1 (MQ=255)
                                                                                                       tAAAGGTGAGTTGCTGTTCTCCGATCCGGATTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCg                                               <  1:510646/150‑1 (MQ=255)
                                                                                                       tAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCg                                               >  2:6640313/1‑150 (MQ=255)
                                                                                                             tGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATcc                                         >  2:7370895/1‑150 (MQ=255)
                                                                                                              gAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCg                                        >  1:6948660/1‑150 (MQ=255)
                                                                                                                 ttGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTc                                     >  2:6933268/1‑150 (MQ=255)
                                                                                                                       ttCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACg                               >  1:8465925/1‑150 (MQ=255)
                                                                                                                       ttCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACg                               >  1:3070244/1‑150 (MQ=255)
                                                                                                                          tcCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTg                            >  2:862387/1‑150 (MQ=255)
                                                                                                                          tcCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTg                            >  2:10754819/1‑150 (MQ=255)
                                                                                                                          tcCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTg                            <  1:7967494/150‑1 (MQ=255)
                                                                                                                           ccGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTgg                           <  2:7838868/150‑1 (MQ=255)
                                                                                                                            cGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGc                          >  2:190892/1‑150 (MQ=255)
                                                                                                                            cGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGc                          >  2:7153125/1‑150 (MQ=255)
                                                                                                                                ccGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaa                      <  1:2826025/150‑1 (MQ=255)
                                                                                                                                 cGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaaa                     >  1:9761710/1‑150 (MQ=255)
                                                                                                                                 cGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaaa                     >  2:11079394/1‑150 (MQ=255)
                                                                                                                                  ggAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAAc                    >  1:9961207/1‑150 (MQ=255)
                                                                                                                                        ccATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGa              >  1:9460804/1‑150 (MQ=255)
                                                                                                                                        ccATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGa              >  2:10626022/1‑150 (MQ=255)
                                                                                                                                          aTAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAgc            >  1:8692344/1‑150 (MQ=255)
                                                                                                                                          aTAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAgc            >  2:9526286/1‑150 (MQ=255)
                                                                                                                                          aTAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAgc            >  2:4884778/1‑150 (MQ=255)
                                                                                                                                           tAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAgcg           >  1:1078128/1‑150 (MQ=255)
                                                                                                                                           tAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAgcg           >  2:3119681/1‑150 (MQ=255)
                                                                                                                                           aaaaaaaGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAgcg           <  1:9111541/149‑1 (MQ=255)
                                                                                                                                               aaaGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGAt       <  1:9664274/150‑1 (MQ=255)
                                                                                                                                               aaaGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGAt       <  2:4219492/150‑1 (MQ=255)
                                                                                                                                                aaGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGAtg      >  1:4213403/1‑150 (MQ=255)
                                                                                                                                                aaGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGAtg      >  2:1779608/1‑150 (MQ=255)
                                                                                                                                                   aaaCGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGa   <  2:9006909/150‑1 (MQ=255)
                                                                                                                                                   aaaCGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGa   <  1:4110904/150‑1 (MQ=255)
                                                                                                                                                   aaaCGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGa   <  1:10018569/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGAT  >  NC_000913/4596675‑4596972

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