Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,438,542 C→T 6.3% C748Y (TGT→TAT) ‡ ydbK ← putative pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,438,5420CT 6.3% 1118.0 / 37.8 348C748Y (TGT→TAT) ‡ydbKputative pyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (157/169);  new base T (6/16);  total (163/185)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 357 total aligned reads displayed.

GTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAA  >  NC_000913/1438393‑1438691
                                                                                                                                                     |                                                                                                                                                     
 tttGCTACGGTCGATTTCTGGCAGGTTGAGGAAGACATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACa                                                                                                                                                      <  2:735357/150‑1 (MQ=255)
  ttGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACaa                                                                                                                                                     >  1:2253965/1‑150 (MQ=255)
    gCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAAcc                                                                                                                                                   <  2:1626732/150‑1 (MQ=255)
     cTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCg                                                                                                                                                  >  2:2144360/1‑150 (MQ=255)
      tACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCgg                                                                                                                                                 <  1:2338072/150‑1 (MQ=255)
        cGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTg                                                                                                                                               >  2:218139/1‑150 (MQ=255)
          gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCa                                                                                                                                             <  2:747883/150‑1 (MQ=255)
            cGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtc                                                                                                                                           <  1:1114450/150‑4 (MQ=255)
            cGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtc                                                                                                                                           <  2:2443676/150‑4 (MQ=255)
            cGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAAt                                                                                                                                           <  2:474166/150‑1 (MQ=255)
             gATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcg                                                                                                                                          <  2:2460631/150‑5 (MQ=255)
             gATTTCTGGCAGGTGGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcg                                                                                                                                          <  2:2044493/150‑5 (MQ=255)
                 tCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTcc                                                                                                                                      <  2:1892103/150‑1 (MQ=255)
                 tCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTcc                                                                                                                                      >  2:143465/1‑150 (MQ=255)
                  cTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCg                                                                                                                                     >  1:1526189/1‑150 (MQ=255)
                  cTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCg                                                                                                                                     <  2:2318359/150‑1 (MQ=255)
                   tGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCgg                                                                                                                                    <  2:2083578/150‑1 (MQ=255)
                   tGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCgg                                                                                                                                    <  2:2070963/150‑1 (MQ=255)
                     gCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTg                                                                                                                                  <  2:2265601/150‑1 (MQ=255)
                     gCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTg                                                                                                                                  <  2:2447570/150‑1 (MQ=255)
                        ggTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCa                                                                                                                               >  1:2391749/1‑150 (MQ=255)
                          ttGAGGAAGAAATCTTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCAcc                                                                                                                             <  2:1093791/150‑1 (MQ=255)
                           tGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCt                                                                                                                            >  1:2082171/1‑150 (MQ=255)
                           tGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCt                                                                                                                            >  1:1147854/1‑150 (MQ=255)
                            gAGGAAGAAATCGTGATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTg                                                                                                                           >  1:2458199/1‑150 (MQ=255)
                              ggAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCa                                                                                                                         <  2:31778/150‑1 (MQ=255)
                              ggAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCa                                                                                                                         <  1:2160245/150‑1 (MQ=255)
                                aagaaATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAg                                                                                                                       <  2:1473148/150‑1 (MQ=255)
                                  gaaATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGAc                                                                                                                     <  2:277002/150‑1 (MQ=255)
                                  gaaATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTATCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCAcc                                                                                                                             >  1:542066/1‑142 (MQ=255)
                                  gaaATCGTAATTGATTTTCTCTCCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGAc                                                                                                                     <  2:1249625/150‑1 (MQ=255)
                                   aaaTCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACa                                                                                                                    <  1:1829591/150‑1 (MQ=255)
                                    aaTCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATc                                                                                                                                          >  2:613451/1‑127 (MQ=255)
                                     aTCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACata                                                                                                                  >  1:499044/1‑150 (MQ=255)
                                         tAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTc                                                                                                              <  1:1792089/150‑1 (MQ=255)
                                           aTTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTg                                                                                                            <  2:172689/150‑1 (MQ=255)
                                            ttGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGc                                                                                                           <  1:2471472/150‑1 (MQ=255)
                                            ttGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGc                                                                                                           >  2:1758475/1‑150 (MQ=255)
                                                    ctctTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCat                                                                                                   >  2:1142436/1‑150 (MQ=255)
                                                      cttcttCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCatat                                                                                                 >  1:986313/1‑150 (MQ=255)
                                                         cttcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCAc                                                                                              >  1:1401750/1‑150 (MQ=255)
                                                           tctaCATGTTCCAGTGGAGAGATGATATTGTTTGCTTTAATCTCTGGATTCTGAAGGTCTTTCTGCGGGGAAACTTCGACGCAAAGGTTACAACCGGTGCAATCTTCCGGTGGCACCTGCAAGAGATATTTCTGCCCGCTCATATCAcgc                                                                                            <  2:1848673/147‑1 (MQ=255)
                                                           tcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCAcgc                                                                                            <  2:1558165/150‑1 (MQ=255)
                                                             gACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGa                                                                                          >  2:39921/1‑150 (MQ=255)
                                                              aCATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGAt                                                                                         <  2:194919/150‑1 (MQ=255)
                                                               cATGTTCTAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGAtt                                                                                        <  1:627521/150‑1 (MQ=255)
                                                                  gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                                     <  2:1005929/150‑1 (MQ=255)
                                                                  gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                                     >  1:1216890/1‑150 (MQ=255)
                                                                  gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                                     <  2:1333928/150‑1 (MQ=255)
                                                                  gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                                     <  1:212573/150‑1 (MQ=255)
                                                                   ttCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCGGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTcac                                                                                    >  1:1330853/1‑150 (MQ=255)
                                                                    tCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTcaca                                                                                   >  2:44488/1‑150 (MQ=255)
                                                                       aGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATcc                                                                                <  1:1236236/150‑1 (MQ=255)
                                                                             gaCATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGaa                                                                          <  1:22997/150‑1 (MQ=255)
                                                                               catcatATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATg                                                                        >  1:82947/1‑150 (MQ=255)
                                                                               catcatATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCCGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATg                                                                        <  1:2090940/150‑1 (MQ=255)
                                                                               catcatATTGATGGCTTTAATCTCTGGAATCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATg                                                                        >  1:1008753/1‑150 (MQ=255)
                                                                                 tcatATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCa                                                                      >  1:1187563/1‑150 (MQ=255)
                                                                                  catATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAg                                                                     >  2:1221850/1‑150 (MQ=255)
                                                                                          tGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAg                                                             >  2:2381024/1‑150 (MQ=255)
                                                                                            gCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAggg                                                           >  2:683248/1‑150 (MQ=255)
                                                                                             cTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAgggg                                                          <  2:1351035/150‑1 (MQ=255)
                                                                                                 aaTCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtt                                                      >  2:1500791/1‑150 (MQ=255)
                                                                                                 aaTCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtt                                                      >  2:253298/1‑150 (MQ=255)
                                                                                                   tctcTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtttt                                                    >  2:912951/1‑150 (MQ=255)
                                                                                                     tctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtttt                                                    >  1:1227024/1‑148 (MQ=255)
                                                                                                      ctGGATTCTGCCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCa                                                 >  2:971855/1‑150 (MQ=255)
                                                                                                      ctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCa                                                 <  1:194857/150‑1 (MQ=255)
                                                                                                              tGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCa                                         >  1:1739070/1‑150 (MQ=255)
                                                                                                                    tCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTCTCCATCGCTTCAggcggc                                   >  2:600333/1‑150 (MQ=255)
                                                                                                                     cTTTCGCCGGGCAAACTTCGACGCACAGGTTACGACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGca                                  <  1:1875683/150‑1 (MQ=255)
                                                                                                                     cTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGca                                  >  1:2009369/1‑150 (MQ=255)
                                                                                                                          gCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACt                             >  1:1749467/1‑150 (MQ=255)
                                                                                                                          gCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACt                             <  1:360283/150‑1 (MQ=255)
                                                                                                                           ccGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACtt                            >  2:197249/1‑150 (MQ=255)
                                                                                                                           ccGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACtt                            >  1:1449072/1‑150 (MQ=255)
                                                                                                                            cGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACttt                           <  1:1627963/150‑1 (MQ=255)
                                                                                                                               gCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgc                        <  1:1487337/150‑1 (MQ=255)
                                                                                                                                cAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgcg                       >  2:490090/1‑150 (MQ=255)
                                                                                                                                 aaaCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgcgc                      <  1:1042577/150‑1 (MQ=255)
                                                                                                                                 aaaCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgcgc                      >  1:14140/1‑150 (MQ=255)
                                                                                                                                   aCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGa                    <  1:302315/150‑1 (MQ=255)
                                                                                                                                   aCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGa                    <  2:1809410/150‑1 (MQ=255)
                                                                                                                                   aCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGa                    <  2:455127/150‑1 (MQ=255)
                                                                                                                                    cTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa                   <  1:1890391/150‑1 (MQ=255)
                                                                                                                                    cTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa                   >  2:771740/1‑150 (MQ=255)
                                                                                                                                    cTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa                   >  2:846507/1‑150 (MQ=255)
                                                                                                                                      agattgCACCGGTTGTAACCTGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATa                 >  1:735978/6‑150 (MQ=255)
                                                                                                                                       cGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAg                <  1:198010/150‑1 (MQ=255)
                                                                                                                                        gACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGc               <  1:1883514/150‑1 (MQ=255)
                                                                                                                                        attgCACCGGTTGTAACCTGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGc               >  1:654850/4‑150 (MQ=255)
                                                                                                                                          cGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTg             >  1:248578/1‑150 (MQ=255)
                                                                                                                                          cGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTg             <  2:2275542/150‑1 (MQ=255)
                                                                                                                                            cacaGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAg           >  1:790329/1‑150 (MQ=255)
                                                                                                                                            cacaGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCAGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAg           >  2:1108073/1‑150 (MQ=255)
                                                                                                                                                  ttACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTggg     <  1:774165/150‑1 (MQ=255)
                                                                                                                                                  ttACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTggg     <  1:2438335/150‑1 (MQ=255)
                                                                                                                                                  ttACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTggg     <  2:2423519/150‑1 (MQ=255)
                                                                                                                                                    aCAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCa   <  2:780602/150‑1 (MQ=255)
                                                                                                                                                     cAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCaa  >  1:966616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAA  >  NC_000913/1438393‑1438691

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