Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,578 C→T N30N (AAC→AAT)  nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,5780CT97.9% 138.7 / NA 47N30N (AAC→AAT) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCCGTCAAAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCAT  >  NC_000913/17335‑17816
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
aCAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATGTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGTGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAgcggc                                                                                                                                                                                                                                          <  2:291636/251‑1 (MQ=255)
          aaaa‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGCTTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGACATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTgg                                                                                                                                                                                                                               >  2:612077/1‑251 (MQ=255)
                       aCCGCTGATGGCGCAAATTCTTCAATGGCTCGTCAAAAACGAGTTGTTCGTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGGCAAATAGTAACAACGATCTATTCACCTGAAAGAGAAATACAAAATGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCCACCAGTGGATGGTATCACgc                                                                                                                                                                                                                   >  2:34682/1‑250 (MQ=255)
                                     aaaTTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACATGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGGGACGc                                                                                                                                                                                                     >  1:510549/1‑251 (MQ=255)
                                                 ggatCGTAAAAAACGAGTTAATCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTAGTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGAGTCTTTAGCAGTGATGCCTCGGGAGGCAGTATTCTTATCATTGCCGCTATCCTGGCGATGATGATGTCCAATAGCGGCGCAACCAGTGTATGGTATTACGTCTTTCTGGATACGcgcgtgcagctc                                                                                                                                                                                          <  2:95923/247‑12 (MQ=255)
                                                   atCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCgg                                                                                                                                                                                       <  1:613947/250‑1 (MQ=255)
                                                     cGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCggg                                                                                                                                                                                      <  2:386136/250‑1 (MQ=255)
                                                                    attCCTACACTATAATCTGAGTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACTCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                                                                                                                                                                       <  1:476368/250‑1 (MQ=255)
                                                                          acacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGCTTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaaca                                                                                                                                                                 >  2:497519/1‑250 (MQ=255)
                                                                          acacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                                                                                                                                                                <  1:124606/251‑1 (MQ=255)
                                                                          acacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaaca                                                                                                                                                                 >  1:158800/1‑250 (MQ=255)
                                                                           cacTATCATATGATTTTCACGATGAATTGTGAGCGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGCGAAATCAAAAGTGCAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGCTGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGCCGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                                                                                                                                                                >  2:566583/1‑250 (MQ=255)
                                                                           cacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGACCCATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCATTACCGGCGCACCCAGTGGCTGTTATCACGCCTTTCTGAACCCGCCGGTTCAGCTCCGGGTTGTTTCACTCgcattcaccaa                                                                                                                                                                >  1:79715/1‑240 (MQ=255)
                                                                               ataATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAATAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACaaaaa                                                                                                                                                             <  2:491399/249‑1 (MQ=255)
                                                                                  aactgaatttaaagatgaattgggaggggGAAAATAGTAAAAAAGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCGGCATCGATTCTTTAGCAGGGATGCCTCGGGAGGCAGTATTCTTATCATTGCCGCTATCCTGGCGATGAGTATGGCCAATAGCGGCGCAACCAGTGGATGGTAGCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACGCGAAAGCAACAAAAACATGc                                                                                                                                                        <  2:79715/226‑1 (MQ=255)
                                                                                              cGATGATTTGTGCGGGGTACAATAGTAAAAACGATCTATTCACCTGAAAGTGAAATAAAAAGTGAAACATCTGCATAGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAATTGCTGTTATGGATa                                                                                                                                             >  1:304943/1‑250 (MQ=255)
                                                                                                atgatTTGTGCGGGGTAAAATAGTAAGAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAAt                                                                                                                                          >  2:116922/1‑251 (MQ=255)
                                                                                                atgatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAAt                                                                                                                                          <  1:347926/251‑1 (MQ=255)
                                                                                                  gatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGACATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGCAATCAACAAAAACATGCTGTTATGGATAAATg                                                                                                                                         >  2:486971/1‑250 (MQ=255)
                                                                                                  gatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATg                                                                                                                                         >  1:219454/1‑250 (MQ=255)
                                                                                                   aatggtgcagtgTAAAATAGTAAAAACTATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGAATCGATTCTTGAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGTCTATGATTATGGCCAATAGCGGCGCAACCAGTGGGTTGTATCACGACTTTCTGGAGACGCCTGTTCATCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATgtc                                                                                                                                       <  2:5222/240‑3 (MQ=255)
                                                                                                         gCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCt                                                                                                                                  <  2:497876/250‑1 (MQ=255)
                                                                                                                aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCAGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCCGTTCAGCCCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCgg                                                                                                                          <  1:34682/251‑1 (MQ=255)
                                                                                                                      gTAAAAACGCTCTCTTCACCTGAAAGTGACATAAACACTGAAACATCTGCATCGCTTCTTTAGCAGACATGACTCGGGAGGCATTATTCTTATCATTGACGCCATCCTGGCGATGATTATGGCAAATAGCGGCGCAACCCGCGGATGCTATCAAGCCTTTCTGGAAACGCAGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAACAACATGATGTTATGGCTAAATGACGCGCTAATGGCAGtttttt                                                                                                                    >  2:269607/1‑251 (MQ=255)
                                                                                                                                      ccccttaaaGAGAAATAAAAAGTGAAACATCTGCATCTATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGTCCAATAGCGGCGCAACCAGTGGATGGTATCACTACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATCCTGTTa                                                                                                                                                   <  2:265917/198‑1 (MQ=255)
                                                                                                                                      cACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGAAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTTACTCGAAATCAACAAAAACATGCTGTTa                                                                                                                                                   >  1:265917/1‑204 (MQ=255)
                                                                                                                                        ccTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTgg                                                                                                  <  2:24580/251‑1 (MQ=255)
                                                                                                                                          tGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCACGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTTGTCGGTCTGGa                                                                                                 >  2:109828/1‑250 (MQ=255)
                                                                                                                                                          aGTGAAAAATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGCTGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATAACGACTTTCTGGGGACGCCGGTTCAGCTCCGGGTGGGTTCACTCGACATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTTGTCGGTCTGGCAGTTTAACGTGTACta                                                                                 >  1:468936/1‑249 (MQ=255)
                                                                                                                                                           gTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGAt                                                                               <  1:531747/251‑1 (MQ=255)
                                                                                                                                                                 cATCTGCATAGTTTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCCTTGCCGCTATCCTGGCGATGATTATGGCCAATAGCAGCGCAACCAGTGGATGGTATCCCGACTTTCTGGAGACGCCGGTACAGCTCCGGGTTGGTTCACTCGAACTCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGATAGTCGGTTTGGAAGTTAAACGGGAACTGATGCAAgg                                                                         >  2:307827/1‑251 (MQ=255)
                                                                                                                                                                        aTCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCt                                                                   >  1:282333/1‑250 (MQ=255)
                                                                                                                                                                           gATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGcc                                                               >  1:470273/1‑251 (MQ=255)
                                                                                                                                                                                                         ttattCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGAAGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATc                                 >  1:207497/1‑251 (MQ=255)
                                                                                                                                                                                                          tattCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATc                                 <  1:486971/250‑1 (MQ=255)
                                                                                                                                                                                                                tATCCTTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGACATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGCTCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCtt                          >  2:517061/1‑250 (MQ=255)
                                                                                                                                                                                                                tATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCt                           >  1:454934/1‑250 (MQ=255)
                                                                                                                                                                                                                   cATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTTCGCCAGGCCGCATTTCCCGTTATCGCCGCTAttt                       >  2:555463/1‑250 (MQ=255)
                                                                                                                                                                                                                    aTTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATtgg                      <  2:219454/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATtgg                      <  1:356044/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccGCGATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGTATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGTCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTggg                  <  2:119527/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGATGGGTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCCAGGGTCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGga                 >  2:346816/1‑247 (MQ=255)
                                                                                                                                                                                                                               ccTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATtg           <  2:408485/251‑1 (MQ=255)
                                                                                                                                                                                                                               ccTGGCGATGATTATGGCCAACAGCGGCTCAACAAGGGTATGGTATCACGCCTTTCTGGTGACGACCGTTCAGCTCCGGGTTTGTTCACTCGAACTCAATAAAAACATGCTGTTATGGGTAAATGTCGCGCTGATGGCTGTAGTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATtg           <  1:329045/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gatgatTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTTGTTCACTCGAAATCAACAACAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTATTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGTATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGTGATGATTGTTCCgg     <  2:78231/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gatgatTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCgg     <  1:420499/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gatTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTTATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAt  <  2:39926/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
ACAGCCGTCAAAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCAT  >  NC_000913/17335‑17816

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