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.

GTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTA  >  NC_000913/4596664‑4596957
                                                                                                                                                     |                                                                                                                                                 
gTTTGGTGGTTCGCTGGCGCAGGAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCa                                                                                                                                                  <  1:2872866/150‑1 (MQ=255)
gTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCa                                                                                                                                                  <  2:2790608/150‑1 (MQ=255)
 tttGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCAt                                                                                                                                                 <  1:3217980/150‑1 (MQ=255)
   tggtggTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaa                                                                                                                                               >  2:2382592/1‑150 (MQ=255)
      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAATGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                                            >  1:476033/1‑150 (MQ=255)
      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                                            >  1:3210027/1‑150 (MQ=255)
      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                                            >  2:278858/1‑150 (MQ=255)
       ggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaaa                                                                                                                                           <  2:2171193/150‑1 (MQ=255)
         ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                         <  1:1128395/150‑1 (MQ=255)
         ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                         <  1:1128167/150‑1 (MQ=255)
           cGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGaaa                                                                                                                                       <  1:2899420/150‑1 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAGACGACGATGATGGGCa                                                                                                                        >  2:1164477/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGACACGACGATGATGGGCa                                                                                                                        >  1:2103622/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:2621136/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:324692/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:1164464/1‑150 (MQ=255)
                          tGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                        >  2:118710/1‑150 (MQ=255)
                             tCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGcc                                                                                                                     >  2:3091904/1‑150 (MQ=255)
                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                   >  2:318387/1‑150 (MQ=255)
                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                   >  2:318415/1‑150 (MQ=255)
                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                   >  2:152057/1‑150 (MQ=255)
                                  gcgACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAac                                                                                                                >  1:222393/1‑150 (MQ=255)
                                  gcgACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAac                                                                                                                >  2:1028244/1‑150 (MQ=255)
                                  gcgACAGGTAATAAACTTGCTATGAATAACACCGTGAATTTTTTTCGTCACGGCGGCAGGGTGGTATTTGTCGGCCTGTTTAACGGTGAGTTGCATTTCTCCGATCCGGAATTCCATAAAAAATAAACGACGATGATGGGCAGCCGCAac                                                                                                                <  2:1750150/150‑1 (MQ=255)
                                         gaaaTCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGACGGCCTGTTTAAAGGTGAGTTGCAGTTCTCAGATCCGGAATTCAATAAAAAAGAAACGACGATGATGGGAAGCCGCAacgccacgc                                                                                                         >  2:2812196/3‑150 (MQ=255)
                                          tAATCACCTTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgcc                                                                                                        <  2:2135588/150‑1 (MQ=255)
                                          tAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGAGGATGGGCAGCCGCAacgccacgcc                                                                                                        >  1:1316310/1‑150 (MQ=255)
                                            aTCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCgg                                                                                                      >  2:717944/1‑150 (MQ=255)
                                               aaCATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAg                                                                                                   >  1:1255678/1‑150 (MQ=255)
                                                aCATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAATGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGa                                                                                                  >  1:2111492/1‑150 (MQ=255)
                                                  aTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGAtt                                                                                                >  2:101960/1‑150 (MQ=255)
                                                     cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                             >  2:1501644/1‑150 (MQ=255)
                                                     cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                             >  2:1498886/1‑150 (MQ=255)
                                                     cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTCGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                             >  2:2488675/1‑150 (MQ=255)
                                                      gATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGc                                                                                            >  2:42946/1‑150 (MQ=255)
                                                       aTGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTCTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCt                                                                                           >  1:1808334/1‑150 (MQ=255)
                                                       aTGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCt                                                                                           >  1:2232891/1‑150 (MQ=255)
                                                        tGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTa                                                                                          >  1:626914/1‑150 (MQ=255)
                                                          aaTAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTaaa                                                                                        >  1:2441512/1‑150 (MQ=255)
                                                          aaTAACACCGTGAATTTGATTCGTCACGGCGGCACGGGGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTaaa                                                                                        >  1:1899404/1‑150 (MQ=255)
                                                                   gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                               >  1:1535230/1‑150 (MQ=255)
                                                                   gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                               >  2:1171112/1‑150 (MQ=255)
                                                                   gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                               >  1:542059/1‑150 (MQ=255)
                                                                        tttGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGAt                                                                          >  1:1755713/1‑150 (MQ=255)
                                                                            aTTCTTCACGGCGCACTGTTGTTATTTCTCGGAATTTTTAAAGGTTTTTTCAAGTACACCTAACCGTAATACAAAAAAAAAAAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCg                                                                      <  1:429183/150‑1 (MQ=255)
                                                                            aTTCGTCCGGGGGGCCGGGTGGTTTTTGTCGGCCTGTTTAATGGTGCTTGGCTGTACTCCCACCCGTAATTCAATAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCg                                                                      <  1:2573209/150‑1 (MQ=255)
                                                                             ttCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAACGGTGAGTTGCAGTTCTCCGAGCCGGAATTCCACAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCgg                                                                     <  2:753008/150‑1 (MQ=255)
                                                                                 tCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaagg                                                                 >  1:1917670/1‑150 (MQ=255)
                                                                                         gcACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                                                         >  1:1399373/1‑150 (MQ=255)
                                                                                         gcACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                                                         >  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)
                                                                                                                                                     |                                                                                                                                                 
GTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTA  >  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: