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,5780CT95.0% 119.5 / ‑3.8 40N30N (AAC→AAT) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (27/11);  total (28/12)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAC‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATT  >  NC_000913/17352‑17817
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
aTCTCACCGCTGATGGCGCAAATTCTTCACTGGATCGTAAAAAACGAGGTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATCGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGCTTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTAt                                                                                                                                                                                                                         >  1:336392/1‑251 (MQ=255)
       cGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGCGGCATTATTCTTATCATTGCCGCTATCCTGGCGTTGATTATGTCCAAT‑AGCGGCGCAACCAGTggtagggagcccgact                                                                                                                                                                                                                  >  2:507176/1‑237 (MQ=255)
            aTGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTg                                                                                                                                                                                                             >  1:132114/1‑251 (MQ=255)
                                             gagttattCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTc                                                                                                                                                                            <  1:150418/248‑1 (MQ=255)
                                                   ttCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATc                                                                                                                                                                      >  2:408826/1‑251 (MQ=255)
                                                                  cTGATTATAACGATGAATTGTGAGAGGAAACATAGTAAACAAGATCTATTCAACTGAAAGCGAAATAAAAAGTGAAACATCAGCAACAATTCTTTAGCAGTGATGCATCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGCTGGTATCACGTCTTTCTGGAGAAGCCGGTTAAGCTCCGAGATGGTTCACTCGCAATccgagaaaacatcatg                                                                                                                                                       >  2:382493/1‑236 (MQ=255)
                                                                    gATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGACACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATGGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAATCATGCTGtt                                                                                                                                                     <  2:225794/251‑1 (MQ=255)
                                                                       tttAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGGTGGTTCACTCGAAATCAACAAAAACATGCTGTTATg                                                                                                                                                  >  1:202592/1‑251 (MQ=255)
                                                                             gatgatTTGTGAGGGGTAAAATAGCAAAAACGATCTATTCACCTGAAAGAGAAATAACAAGTGAAACATCTACATCAATTCTTTAGAAGTGCTGCCTCGGGAGGAATTATTCTTCTCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGCTGGTATCACGACTTTCTGAAGACACAGGTTCCGCTCCGGGTTGTTTACCTCACAAAAAACAAAAACATGCTGTTATGGaaaa                                                                                                                                             >  2:562080/1‑247 (MQ=255)
                                                                              atgatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGCTGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAAt                                                                                                                                           >  1:86043/1‑251 (MQ=255)
                                                                                              aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGGAATCAACAAAAACATGCTGt                                                                                                                                                      >  2:196482/1‑224 (MQ=255)
                                                                                              aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGt                                                                                                                                                      <  1:196482/224‑1 (MQ=255)
                                                                                              aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAC‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCg                                                                                                                            <  1:12808/250‑1 (MQ=255)
                                                                                              aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCACCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGACGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCgg                                                                                                                           <  2:66460/251‑1 (MQ=255)
                                                                                               aaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGt                                                                                                                          <  1:8731/251‑1 (MQ=255)
                                                                                                 aTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTAAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCggcat                                                                                                                        >  1:326710/1‑248 (MQ=255)
                                                                                                                 aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTAACATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTa                                                                                                             <  2:502311/246‑1 (MQ=255)
                                                                                                                    cACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCt                                                                                                     >  1:326576/1‑251 (MQ=255)
                                                                                                                    cACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCt                                                                                                     >  1:293227/1‑251 (MQ=255)
                                                                                                                     aCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTAGTGAGGCGTTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACAACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGGATTTTTCCTGTTAGTCGGTCTg                                                                                                    >  1:176102/1‑251 (MQ=255)
                                                                                                                        tgcaaGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATACCTCGGGAGGCATTATTCTTATCATTGACGCTATCCTGGCGATGATTATGGC‑AACTAGCGGCGCAACCAGTGGATGGTATAACGAATTTCTGGAGACGCCGGTTCAGCTACGGGTTGGTTAACTCGAAGTCAACAAAAACATGCGGTTATGGATAAATGCCGCGCTGATGGCGGAATTTTTCCTGTTCGTCTGTcgtga                                                                                                  >  2:603276/4‑246 (MQ=255)
                                                                                                                        tGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGGATTTTTCCTGTTAGGCCGTCTggc                                                                                                  >  1:422586/1‑249 (MQ=255)
                                                                                                                                                    gcaccGATTCTTTAGCAGTGATGCCACGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATc                                                                       <  2:132114/245‑1 (MQ=255)
                                                                                                                                                      aTCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCt                                                                    >  1:472212/1‑250 (MQ=255)
                                                                                                                                                            tCTTTAGCAGTGATGCATCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGGATTTTTCCTTTTAGTAGGTCTGGGAGTTAAACGTGAACCGATGAACGGATCGCTAGCCAGc                                                             >  2:85317/1‑251 (MQ=255)
                                                                                                                                                                                ggAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTcc                                         <  1:308093/251‑1 (MQ=255)
                                                                                                                                                                                  aGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAg                                       <  1:408826/251‑1 (MQ=255)
                                                                                                                                                                                               aTCATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGGATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAAATGATGCAAGGAACGCTAGCCAGCTTACGCCAGGCCGGATTTCCAGTTATAGCCGCTAt                          >  1:367840/1‑251 (MQ=255)
                                                                                                                                                                                                 cATTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCCGGCCGCATTTCCAGTTATCGCCGCTATtg                        >  1:469232/1‑251 (MQ=255)
                                                                                                                                                                                                 cATTGCCGCTATCCTGGCGATGATTATGGCAAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGCTAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTAGGTATGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTAttt                        >  1:585238/1‑250 (MQ=255)
                                                                                                                                                                                                  aTTGCCGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATtgg                       <  1:138053/251‑1 (MQ=255)
                                                                                                                                                                                                       cGCTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGGCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTTGTgggg                  >  1:357446/1‑250 (MQ=255)
                                                                                                                                                                                                         cTATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCAATCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGAATTTTTCCTGTTAGTCGGGCTGGAAATTAAACGGGAACTGATGCAAGGATCGATAACCAGCTTACGCCAGGCACCATTTCCAGTTATCACCGCTTTTGGTggggtg                >  2:647145/1‑248 (MQ=255)
                                                                                                                                                                                                          tATCCTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGGGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGgatg                >  1:487534/1‑250 (MQ=255)
                                                                                                                                                                                                             ccTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAAAGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATtg            >  1:26699/1‑251 (MQ=255)
                                                                                                                                                                                                             ccTGGCGATGATTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGGTAGTCGGTCTGGAAGGTAAACGTGAAATGATGGAAGGAACGCTAGCCAGCTTACGCCCGGCCGCATTTCCAGTTATCGCCGCTATTgggggggtggttg            >  2:302690/1‑244 (MQ=255)
                                                                                                                                                                                                                      gatTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAt   <  1:302690/251‑1 (MQ=255)
                                                                                                                                                                                                                      gatTATGGCCAAT‑AGCGGCGCAACCAGAGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGCATTTTTCCTGTTAGTCGGTATGGAAGTTAAACGTGAACTGATGCAAGGGTCGCAAGCCCGCTTACGCCCGGCCGCCTTTCCCATTATCTCCGCTTTTGGGGGGGTGGTTGTTTCGGCAt   >  2:489911/1‑251 (MQ=255)
                                                                                                                                                                                                                       atTATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGGAGTTAAACGTGAACTGATGCAAGGTTCGCTAGCCAGCTTTCGCCAGGCCGCCTTTCCAGTTATCGCCGCTATTTGGGGGATGTGTGTTCCGGCAtt  >  2:8056/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttATGGCCAAT‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCAtt  >  1:85633/1‑250 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
ATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAC‑AGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATT  >  NC_000913/17352‑17817

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