Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 870,892 C→A 6.2% T494T (ACC→ACA) ‡ gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913870,8920CA 6.2% 379.1 / 9.9 108T494T (ACC→ACA) ‡gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (48/53);  new base A (5/2);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 112 total aligned reads displayed.

TGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTT  >  NC_000913/870745‑871040
                                                                                                                                                   |                                                                                                                                                    
tGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCg                                                                                                                                                     <  3:595839/149‑1 (MQ=255)
  aattCCTGGCTTAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGt                                                                                                                                                   <  4:372517/147‑1 (MQ=255)
   cTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                                  >  4:331676/1‑149 (MQ=255)
    ttCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCGGTCAGGttt                                                                                                                                                 >  1:359628/1‑149 (MQ=255)
      ccTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                               <  1:180271/149‑1 (MQ=255)
      ccTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                               <  1:58606/149‑1 (MQ=255)
        tGGCTCAGGCCGTGAAAACTAATGATACGGCGGAAACGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                             >  2:296810/1‑149 (MQ=255)
        tGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                             >  4:820361/1‑149 (MQ=255)
        tGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                             >  1:473579/1‑149 (MQ=255)
         ggCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCGGTCAGGtttttact                                                                                                                                            <  2:359628/149‑4 (MQ=255)
         ggCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGAGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTATAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                            <  4:768077/149‑1 (MQ=255)
          gCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGAt                                                                                                                                           <  3:424239/149‑1 (MQ=255)
           cTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGAGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATc                                                                                                                                          <  4:116451/149‑1 (MQ=255)
             cAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCAt                                                                                                                                        <  4:509179/149‑1 (MQ=255)
             cAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCAt                                                                                                                                        <  1:183206/149‑1 (MQ=255)
              aGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                                       >  3:766354/1‑149 (MQ=255)
               ggCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGc                                                                                                                                      <  1:457638/149‑1 (MQ=255)
                gCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGcc                                                                                                                                     >  4:724071/1‑149 (MQ=255)
                 cACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCa                                                                                                                                    <  3:309471/149‑1 (MQ=255)
                 cACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCa                                                                                                                                    <  3:562952/149‑1 (MQ=255)
                    tGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGac                                                                                                                                 >  3:842781/1‑149 (MQ=255)
                     gAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGaca                                                                                                                                >  1:190795/1‑149 (MQ=255)
                     gAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGaca                                                                                                                                >  3:919196/1‑149 (MQ=255)
                       aaaCTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACAcc                                                                                                                              <  3:403541/149‑1 (MQ=255)
                        aaCTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCg                                                                                                                             >  3:582796/1‑149 (MQ=255)
                                aTCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGc                                                                                                                     >  2:351374/1‑149 (MQ=255)
                                    ggcggAAAGGGCCCGCTTATGTAAGGCGGCGCCGTATAGCATCGGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAAAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGCCTTCAGCTTTg                                                                                                                 >  4:441848/1‑149 (MQ=255)
                                    ggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTg                                                                                                                 >  2:486432/1‑149 (MQ=255)
                                     gcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGa                                                                                                                >  4:8964/1‑149 (MQ=255)
                                      cggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGaa                                                                                                               <  4:718184/149‑1 (MQ=255)
                                        gAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGa                                                                                                             <  4:391390/149‑1 (MQ=255)
                                        gAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGa                                                                                                             >  2:640935/1‑149 (MQ=255)
                                         aaaaGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                               <  4:813663/114‑1 (MQ=255)
                                           aaGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAcgc                                                                                                          <  3:627483/149‑1 (MQ=255)
                                           aaGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAcgc                                                                                                          >  1:334785/1‑149 (MQ=255)
                                           aaGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAcgc                                                                                                          <  3:688245/149‑1 (MQ=255)
                                             gACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAACAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCgg                                                                                                        <  2:571728/149‑1 (MQ=255)
                                             gACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCgg                                                                                                        >  1:632959/1‑149 (MQ=255)
                                                ccGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAtt                                                                                                     >  2:378746/1‑149 (MQ=255)
                                                ccGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAtt                                                                                                     >  2:808340/1‑149 (MQ=255)
                                                  gCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTg                                                                                                   <  3:8964/149‑1 (MQ=255)
                                                  gCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTg                                                                                                   <  4:515589/149‑1 (MQ=255)
                                                       tataAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATa                                                                                              >  3:177455/1‑149 (MQ=255)
                                                        ataAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATaa                                                                                             <  2:244095/149‑1 (MQ=255)
                                                                  gcAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAgg                                                                                   <  1:239766/149‑1 (MQ=255)
                                                                      gATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                               >  1:692643/1‑149 (MQ=255)
                                                                       atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGa                                                                                                                <  3:114970/115‑1 (MQ=255)
                                                                       atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGa                                                                                                                >  4:114970/1‑115 (MQ=255)
                                                                       atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTg                                                                              <  1:556944/149‑1 (MQ=255)
                                                                       atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTg                                                                              <  2:741525/149‑1 (MQ=255)
                                                                       atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTg                                                                              <  2:873868/149‑1 (MQ=255)
                                                                        tatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTTCAATAAGCAACGCAGGGAGTgg                                                                             >  1:847582/1‑149 (MQ=255)
                                                                         atcatcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                            >  2:734597/1‑149 (MQ=255)
                                                                          tcatcTGGCAAGCATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGTCGCGGATTTGCAATAAGCAACGCAGGGAGTGGaa                                                                           <  4:582796/149‑1 (MQ=255)
                                                                          tcatcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGaa                                                                           <  1:506398/149‑1 (MQ=255)
                                                                          tcatcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGaa                                                                           >  3:381608/1‑149 (MQ=255)
                                                                            atcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                         >  3:420754/1‑149 (MQ=255)
                                                                               tGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCtt                                                                      <  4:766354/149‑1 (MQ=255)
                                                                                  cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                                   <  1:296810/149‑1 (MQ=255)
                                                                                  cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                                   <  3:820361/149‑1 (MQ=255)
                                                                                  cAAGAATCGCCGTGGATCCCCCTGGTGGTACAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATTATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGGAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                                   <  3:596930/149‑1 (MQ=255)
                                                                                     gAATCGCCGTGGATCCGGCTGGTGGTAGAAAAACTGGTATCGGCACACAGTAAAAACCAGACCGGTTTTTGGATCATGCCAGACACCTGCTTCAGCTTGGAAGACGCGCATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTAc                                                                >  3:777256/1‑149 (MQ=255)
                                                                                     gAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTAc                                                                <  4:381608/149‑1 (MQ=255)
                                                                                      aaTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                               <  4:248866/149‑1 (MQ=255)
                                                                                                  tCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCt                                                   >  3:627276/1‑149 (MQ=255)
                                                                                                  tCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCt                                                   >  2:880712/1‑149 (MQ=255)
                                                                                                    ccGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTa                                                 >  3:520094/1‑149 (MQ=255)
                                                                                                         ggtggtAGAAAGACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTggg                                            >  1:762114/1‑149 (MQ=255)
                                                                                                            ggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGtt                                         <  1:880712/149‑1 (MQ=255)
                                                                                                            ggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGAACATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGCAATGCTTAATTACGTTATCAAACGCTTACTGGGGtt                                         >  4:841582/1‑149 (MQ=255)
                                                                                                             gtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTg                                        >  3:791818/1‑149 (MQ=255)
                                                                                                              tAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGa                                       >  2:432230/1‑149 (MQ=255)
                                                                                                               aGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAt                                      >  1:468479/1‑149 (MQ=255)
                                                                                                                gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAtt                                     >  3:899543/1‑149 (MQ=255)
                                                                                                                   aaaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCg                                  >  2:115755/1‑149 (MQ=255)
                                                                                                                       tGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGc                              <  3:447669/149‑1 (MQ=255)
                                                                                                                       tGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGc                              >  3:157970/1‑149 (MQ=255)
                                                                                                                           gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGtt                          <  4:65797/149‑1 (MQ=255)
                                                                                                                            tCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCaaa                                                       >  3:530707/1‑119 (MQ=255)
                                                                                                                            tCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCaaa                                                       <  4:530707/119‑1 (MQ=255)
                                                                                                                            tCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGttt                         >  2:509573/1‑149 (MQ=255)
                                                                                                                             cGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTa                        <  1:378746/149‑1 (MQ=255)
                                                                                                                               gCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATc                      <  4:842781/149‑1 (MQ=255)
                                                                                                                                 acacaGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGt                    >  1:735562/1‑149 (MQ=255)
                                                                                                                                    caGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTCCTGGGGTTGATTCCGACGCTGTTTATCGtctc                 <  4:885995/149‑1 (MQ=255)
                                                                                                                                       tAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGt              >  4:67025/1‑149 (MQ=255)
                                                                                                                                       caaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTag                                                                >  2:719936/2‑98 (MQ=255)
                                                                                                                                       caaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTAc                                                                <  1:719936/98‑1 (MQ=255)
                                                                                                                                       caaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGt              >  4:539649/2‑149 (MQ=255)
                                                                                                                                        aaaaaCCGGTCATGTTTTTGGATCACGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGGCTCGGTg             >  3:299482/1‑149 (MQ=255)
                                                                                                                                         aaaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGc            >  2:694353/1‑149 (MQ=255)
                                                                                                                                            aCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTg             <  3:232352/145‑1 (MQ=255)
                                                                                                                                                 aCCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttt    <  4:520094/149‑1 (MQ=255)
                                                                                                                                                 aCCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttt    <  2:692643/149‑1 (MQ=255)
                                                                                                                                                 aCCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttt    >  2:104771/1‑149 (MQ=255)
                                                                                                                                                 aCCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttt    <  1:734597/149‑1 (MQ=255)
                                                                                                                                                  ccGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGtttt   <  4:452837/149‑1 (MQ=255)
                                                                                                                                                   cGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttttt  >  1:409606/1‑149 (MQ=255)
                                                                                                                                                   cGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttttt  >  2:720609/1‑149 (MQ=255)
                                                                                                                                                   cGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttttt  >  4:1800/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTT  >  NC_000913/870745‑871040

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