Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,596,812 C→A 9.3% H275N (CAT→AAT)  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,8120CA 9.3% 534.1 / 3.9 167H275N (CAT→AAT) lgoDL‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (106/39);  new base A (2/20);  total (108/59)
Fisher's exact test for biased strand distribution p-value = 1.01e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

GTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTA  >  NC_000913/4596664‑4596957
                                                                                                                                                     |                                                                                                                                                 
gTTTGGTGGTTCGCTGGCGCAGGAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCa                                                                                                                                                  <  1:2872866/150‑1 (MQ=255)
gTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCa                                                                                                                                                  <  2:2790608/150‑1 (MQ=255)
 tttGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCAt                                                                                                                                                 <  1:3217980/150‑1 (MQ=255)
   tggtggTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaa                                                                                                                                               >  2:2382592/1‑150 (MQ=255)
      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAATGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                                            >  1:476033/1‑150 (MQ=255)
      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                                            >  1:3210027/1‑150 (MQ=255)
      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                                            >  2:278858/1‑150 (MQ=255)
       ggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaaa                                                                                                                                           <  2:2171193/150‑1 (MQ=255)
         ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                         <  1:1128395/150‑1 (MQ=255)
         ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                         <  1:1128167/150‑1 (MQ=255)
           cGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGaaa                                                                                                                                       <  1:2899420/150‑1 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAGACGACGATGATGGGCa                                                                                                                        >  2:1164477/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGACACGACGATGATGGGCa                                                                                                                        >  1:2103622/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:2621136/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:324692/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:1164464/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:118710/1‑150 (MQ=255)
                             tCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGcc                                                                                                                     >  2:3091904/1‑150 (MQ=255)
                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                   >  2:318387/1‑150 (MQ=255)
                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                   >  2:318415/1‑150 (MQ=255)
                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                   >  2:152057/1‑150 (MQ=255)
                                  gcgACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAac                                                                                                                >  1:222393/1‑150 (MQ=255)
                                  gcgACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAac                                                                                                                >  2:1028244/1‑150 (MQ=255)
                                  gcgACAGGTAATAAACTTGCTATGAATAACACCGTGAATTTTTTTCGTCACGGCGGCAG‑GGTGGTATTTGTCGGCCTGTTTAACGGTGAGTTGCATTTCTCCGATCCGGAATTCCATAAAAAATAAACGACGATGATGGGCAGCCGCAac                                                                                                                <  2:1750150/150‑1 (MQ=255)
                                         gaaaTCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGACGGCCTGTTTAAAGGTGAGTTGCAGTTCTCAGATCCGGAATTCAATAAAAAAGAAACGACGATGATGGGAAGCCGCAacgccacgc                                                                                                         >  2:2812196/3‑150 (MQ=255)
                                          tAATCACCTTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgcc                                                                                                        <  2:2135588/150‑1 (MQ=255)
                                          tAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGAGGATGGGCAGCCGCAacgccacgcc                                                                                                        >  1:1316310/1‑150 (MQ=255)
                                            aTCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCgg                                                                                                      >  2:717944/1‑150 (MQ=255)
                                               aaCATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAg                                                                                                   >  1:1255678/1‑150 (MQ=255)
                                                aCATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAATGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGa                                                                                                  >  1:2111492/1‑150 (MQ=255)
                                                  aTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGAtt                                                                                                >  2:101960/1‑150 (MQ=255)
                                                     cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                             >  2:1501644/1‑150 (MQ=255)
                                                     cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                             >  2:1498886/1‑150 (MQ=255)
                                                     cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTCGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                             >  2:2488675/1‑150 (MQ=255)
                                                      gATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGc                                                                                            >  2:42946/1‑150 (MQ=255)
                                                       aTGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTCTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCt                                                                                           >  1:1808334/1‑150 (MQ=255)
                                                       aTGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCt                                                                                           >  1:2232891/1‑150 (MQ=255)
                                                        tGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTa                                                                                          >  1:626914/1‑150 (MQ=255)
                                                          aaTAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTaaa                                                                                        >  1:2441512/1‑150 (MQ=255)
                                                          aaTAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGGGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTaaa                                                                                        >  1:1899404/1‑150 (MQ=255)
                                                                   gTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                               >  1:1535230/1‑150 (MQ=255)
                                                                   gTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                               >  2:1171112/1‑150 (MQ=255)
                                                                   gTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                               >  1:542059/1‑150 (MQ=255)
                                                                        tttGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGAt                                                                          >  1:1755713/1‑150 (MQ=255)
                                                                            aTTCTTCACGGC‑GCACTGTTGTTATTTCTCGGAATTTTTAAAGGTTTTTTCAAGTACACCTAACCGTAATACAAAAAAAAAAAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCg                                                                      <  1:429183/150‑1 (MQ=255)
                                                                            aTTCGTCCGGGGGGCCG‑GGTGGTTTTTGTCGGCCTGTTTAATGGTGCTTGGCTGTACTCCCACCCGTAATTCAATAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCg                                                                      <  1:2573209/150‑1 (MQ=255)
                                                                             ttCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAACGGTGAGTTGCAGTTCTCCGAGCCGGAATTCCACAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCgg                                                                     <  2:753008/150‑1 (MQ=255)
                                                                                 tCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaagg                                                                 >  1:1917670/1‑150 (MQ=255)
                                                                                         gcAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                                                         >  1:1399373/1‑150 (MQ=255)
                                                                                         gcAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                                                         >  2:1504844/1‑150 (MQ=255)
                                                                                              ggtggtATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgc                                                     >  2:7424/1‑150 (MQ=255)
                                                                                                 ttttttttTCGGCCTTTTTAACGGTTTGTTGCTTTTATCCGATCGGTTTTACCAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                  <  1:2014820/146‑1 (MQ=255)
                                                                                                 gttttttGTCGTCCTTTTTAACGGTTCGTTGCAGTCCACCTCTCCGTAATACAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                  <  1:2731477/146‑1 (MQ=255)
                                                                                                 gttttttGTCGGCCTTTTTAATGTTGTTTTGCAGTTCCCCTACCCGTAATACAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                  <  1:717944/146‑1 (MQ=255)
                                                                                                 gttttttGGCGGCCTTTTTAAAGGTGAGTTGCCGTTCTCCTATCCGGAATACCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                  <  2:922868/146‑1 (MQ=255)
                                                                                                 ggtTTTTGTCGGCCTGTTTAACGGTGAGTGGCAGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                  <  2:2011775/150‑1 (MQ=255)
                                                                                                 ggtTTTTGTCGCCCTTTTTAATGGTGTGTTGCTGTTCTCCGATCCGTATTTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                  <  2:388353/150‑1 (MQ=255)
                                                                                                    aTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACat                                               >  1:1978311/1‑150 (MQ=255)
                                                                                                    aTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACat                                               >  1:333666/1‑150 (MQ=255)
                                                                                                     tttGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatg                                              >  2:2856609/1‑150 (MQ=255)
                                                                                                      ttGTCGGCCTGTTTAACGGTGAGTTGCAGTTCTCCGATCGGGAATCCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                             <  2:610959/150‑1 (MQ=255)
                                                                                                      ttGTCGGCATGTTTAAGGGTGTGTTGCAGTTCTCCGATCCGTAATTCAATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                             <  1:2912437/150‑1 (MQ=255)
                                                                                                      ttGTCGCCCTTTTTAAAGGTGTTTTGCTGTTCCCCTACCCGTAATTCAAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                             <  1:83566/150‑1 (MQ=255)
                                                                                                       tGTCGGCCTTTTAAATGGTGTTTTGAATTTCACCGATCCGTAATTCAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatgat                                            <  1:1366756/150‑1 (MQ=255)
                                                                                                            gCCTGTTTAATGGTGTGTTGCTGTTCTCCGCTCCGGAATTCCAAAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTaa                                       <  1:3171358/150‑1 (MQ=255)
                                                                                                            gCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTaa                                       >  1:638403/1‑150 (MQ=255)
                                                                                                            gCCTGTTAAACGGTGTGTTGCAGTCCTCCGATCCGGAATTCAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTaa                                       <  1:771886/150‑1 (MQ=255)
                                                                                                             ccTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc                                      <  2:2165649/150‑1 (MQ=255)
                                                                                                             ccTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAAGTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc                                      <  2:528363/150‑1 (MQ=255)
                                                                                                             cattaatAAAGTTGAGTTGCAGTACTCCGATCCCTAATACAAAAAAAAACAAACGACCATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc                                      <  2:476033/144‑1 (MQ=255)
                                                                                                              ctttttAATGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAcc                                     <  1:1893638/147‑1 (MQ=255)
                                                                                                              cTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAcc                                     >  1:892967/1‑150 (MQ=255)
                                                                                                              cTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAcc                                     >  1:892978/1‑150 (MQ=255)
                                                                                                                 tttAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTTCTAAATTCGGTCGACTGATGGCGGAAGGAAAAATAACTGATGACATGATGTTAACCCAt                                  >  2:1314525/1‑150 (MQ=255)
                                                                                                                 tttAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCAt                                  >  1:2549802/1‑150 (MQ=255)
                                                                                                                 tttAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCAt                                  >  2:1137183/1‑150 (MQ=255)
                                                                                                                 tttAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCAt                                  >  2:1314515/1‑150 (MQ=255)
                                                                                                                        gTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATc                           >  1:2404083/1‑150 (MQ=255)
                                                                                                                          gAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCg                         >  1:2675821/1‑150 (MQ=255)
                                                                                                                          gAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCg                         >  1:1144030/1‑150 (MQ=255)
                                                                                                                             ttgcttttCGCCGTTCCGTAATCCAAAAAAAAATACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTc                      <  2:577698/144‑1 (MQ=255)
                                                                                                                             ttGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTc                      >  1:774494/1‑150 (MQ=255)
                                                                                                                             ttGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTc                      >  2:2017152/1‑150 (MQ=255)
                                                                                                                              tGCTGTTCTCCGATCCGGATTTCAATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCg                     <  2:5789/150‑1 (MQ=255)
                                                                                                                              tGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCg                     >  2:2756441/1‑150 (MQ=255)
                                                                                                                              tGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCg                     >  1:267081/1‑150 (MQ=255)
                                                                                                                               gCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGc                    <  2:1899404/150‑1 (MQ=255)
                                                                                                                               gCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGc                    >  2:653651/1‑150 (MQ=255)
                                                                                                                               gCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGc                    >  2:2406495/1‑150 (MQ=255)
                                                                                                                                    tctcCGATCCGGAATTCAATAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGc               <  1:2603322/150‑1 (MQ=255)
                                                                                                                                    tctcCGACCCGTAATACCAAAAAAAACAATCGACGATGATGGGCAGCCGCAAGGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGTCTGATGGCGGATGGAAAAATCACTGCTGACATGATGTTTACCCATCGCTTTCCGTTCGCCACGc               <  1:3028265/150‑1 (MQ=255)
                                                                                                                                      tcCGATCCTTAATTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTg             <  2:626914/150‑1 (MQ=255)
                                                                                                                                       ccGATCCGGAATTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTgg            <  2:1928820/150‑1 (MQ=255)
                                                                                                                                          aTCCGTATTTCAATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAg         <  1:1697606/150‑1 (MQ=255)
                                                                                                                                          aTCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAg         >  1:1869755/1‑150 (MQ=255)
                                                                                                                                          aTCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAg         >  1:1749769/1‑150 (MQ=255)
                                                                                                                                            ccGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaa       >  1:3168282/1‑150 (MQ=255)
                                                                                                                                             cGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaaa      >  2:2717377/1‑150 (MQ=255)
                                                                                                                                             cGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaaa      >  2:2878722/1‑150 (MQ=255)
                                                                                                                                                 cttCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTa  >  2:2354014/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAC‑GGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTA  >  NC_000913/4596664‑4596957

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: ‑