Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 17,528 G→C 5.5% G14R (GGA→CGA)  nhaA → sodium‑proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,5280GC 5.5% 383.4 / 10.3 108G14R (GGA→CGA) nhaAsodium‑proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (46/56);  new base C (3/3);  total (49/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 114 total aligned reads displayed.

ATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGA  >  NC_000913/17381‑17675
                                                                                                                                                   |                                                                                                                                                   
aTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCgg                                                                                                                                                     <  2:657260/148‑1 (MQ=255)
 tAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                    <  2:136475/148‑1 (MQ=255)
 tAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                    >  2:764820/1‑148 (MQ=255)
 tAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGa                                                                                                                                                   <  4:121855/149‑1 (MQ=255)
  aGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCATCGGGa                                                                                                                                                   <  2:529047/148‑1 (MQ=255)
      cGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCa                                                                                                                                               >  1:399991/1‑148 (MQ=255)
       gTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCCCGAGGCatc                                                                                                                                             <  3:621142/149‑2 (MQ=255)
        tAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCatt                                                                                                                                             >  1:269517/1‑148 (MQ=255)
         aaaaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattat                                                                                                                                           >  3:906281/1‑149 (MQ=255)
         aaaaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCatta                                                                                                                                            <  2:874326/148‑1 (MQ=255)
           aaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattatt                                                                                                                                          <  2:793449/148‑1 (MQ=255)
                gAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                                                                                                                    >  4:847986/1‑149 (MQ=255)
                   ttattCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCAtt                                                                                                                                 <  3:523267/149‑1 (MQ=255)
                     attCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTg                                                                                                                                <  1:125365/148‑1 (MQ=255)
                      ttCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGc                                                                                                                               <  1:46866/148‑1 (MQ=255)
                        ccTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCg                                                                                                                             >  1:826587/1‑148 (MQ=255)
                         cTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATTTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCt                                                                                                                           >  3:235980/1‑149 (MQ=255)
                           acacTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTAt                                                                                                                         >  4:608771/1‑149 (MQ=255)
                               tataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTc                                                                                                                                         >  4:912197/1‑129 (MQ=255)
                               tataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTc                                                                                                                                         <  3:912197/129‑1 (MQ=255)
                               tataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCt                                                                                                                      <  2:204942/148‑1 (MQ=255)
                               tataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTg                                                                                                                     >  3:43010/1‑149 (MQ=255)
                                ataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTgg                                                                                                                    <  4:906281/149‑1 (MQ=255)
                                     cTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgat                                                                                                                >  1:839518/1‑148 (MQ=255)
                                      tGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatg                                                                                                               >  1:847546/1‑148 (MQ=255)
                                      tGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatg                                                                                                               <  1:698327/148‑1 (MQ=255)
                                       gATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatgat                                                                                                             <  3:373422/149‑1 (MQ=255)
                                        aTTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatgat                                                                                                             <  2:852911/148‑1 (MQ=255)
                                        aTTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGAtt                                                                                                            >  4:209110/1‑149 (MQ=255)
                                                gatgatTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCaa                                                                                                    >  4:80011/1‑149 (MQ=255)
                                                 atgatTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAAc                                                                                                   <  3:790526/149‑1 (MQ=255)
                                                     ttCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATc                                                                                                                        >  2:319744/1‑124 (MQ=255)
                                                     ttCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATc                                                                                                                        <  1:319744/124‑1 (MQ=255)
                                                     ttCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcg                                                                                               <  4:410526/149‑1 (MQ=255)
                                                       cGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTGTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcggc                                                                                             <  3:699237/149‑1 (MQ=255)
                                                             gggTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAAcc                                                                                       >  3:473024/1‑149 (MQ=255)
                                                                tAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGt                                                                                    >  3:683425/1‑149 (MQ=255)
                                                                tAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGt                                                                                    <  3:902067/149‑1 (MQ=255)
                                                                 aaaaTAGTGAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTg                                                                                   <  3:234141/149‑1 (MQ=255)
                                                                 aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTg                                                                                   <  3:830453/149‑1 (MQ=255)
                                                                   aaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAAc                                                                                        <  1:507556/142‑1 (MQ=255)
                                                                   aaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAAc                                                                                        >  2:507556/1‑142 (MQ=255)
                                                                   aaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTgg                                                                                  <  1:703864/148‑1 (MQ=255)
                                                                    aTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGAt                                                                                <  3:661140/149‑1 (MQ=255)
                                                                      aGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATg                                                                               <  2:630116/148‑1 (MQ=255)
                                                                         aaaaaCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTAt                                                                           >  3:223182/1‑149 (MQ=255)
                                                                           aaaCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATc                                                                          >  1:483763/1‑148 (MQ=255)
                                                                           aaaCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATc                                                                          >  1:429647/1‑148 (MQ=255)
                                                                               gATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGAc                                                                     >  3:126354/1‑149 (MQ=255)
                                                                               gATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGAc                                                                     <  4:79290/149‑1 (MQ=255)
                                                                                  cTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACtt                                                                   >  2:194440/1‑148 (MQ=255)
                                                                                    aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTc                                                                 <  1:675730/148‑1 (MQ=255)
                                                                                        aCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGga                                                             >  2:298862/1‑148 (MQ=255)
                                                                                         ccTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGgag                                                            >  1:966109/1‑148 (MQ=255)
                                                                                            gAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACg                                                         >  1:17883/1‑148 (MQ=255)
                                                                                            gAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGc                                                        <  3:516322/149‑1 (MQ=255)
                                                                                              aaGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatga                                                                                                              >  3:527314/1‑93 (MQ=255)
                                                                                              aaGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatga                                                                                                              <  4:527314/93‑1 (MQ=255)
                                                                                               agagAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCAc                                                                        <  4:151790/130‑1 (MQ=255)
                                                                                               agagAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCg                                                      >  1:833740/1‑148 (MQ=255)
                                                                                                gagaAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGt                                                    <  4:857038/149‑1 (MQ=255)
                                                                                                  gaAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTc                                                  >  3:643640/1‑149 (MQ=255)
                                                                                                   aaaTAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCa                                                 >  4:861815/1‑149 (MQ=255)
                                                                                                        aaaaGTGATACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGAGGAATATGGCCAACCGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTc                                             <  2:826587/148‑1 (MQ=255)
                                                                                                          aaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCg                                           <  1:194440/148‑1 (MQ=255)
                                                                                                           aGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATc                                                                          >  4:27605/1‑116 (MQ=255)
                                                                                                           aGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATc                                                                          <  3:27605/116‑1 (MQ=255)
                                                                                                                 aCATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGggttggtt                                   >  3:5595/1‑149 (MQ=255)
                                                                                                                   aTCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCa                                 >  3:343485/1‑149 (MQ=255)
                                                                                                                    tCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCa                                 <  1:461885/148‑1 (MQ=255)
                                                                                                                     cTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACt                               <  4:235980/149‑1 (MQ=255)
                                                                                                                          tCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                          <  4:390115/149‑1 (MQ=255)
                                                                                                                           cGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                          >  1:178113/1‑148 (MQ=255)
                                                                                                                            gATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCgg                                                     <  3:421119/120‑1 (MQ=255)
                                                                                                                            gATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCgg                                                     >  4:421119/1‑120 (MQ=255)
                                                                                                                            gATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                          >  4:70951/1‑147 (MQ=255)
                                                                                                                            gATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                          <  3:70951/147‑1 (MQ=255)
                                                                                                                            gATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATc                        >  4:575447/1‑149 (MQ=255)
                                                                                                                              ttCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATca                       <  2:839518/148‑1 (MQ=255)
                                                                                                                              ttCTTTAGCAGTAATGCCTCGGGCGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATca                       <  1:301965/148‑1 (MQ=255)
                                                                                                                               tCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaa                      >  1:518640/1‑148 (MQ=255)
                                                                                                                                 tttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                          <  3:787912/142‑1 (MQ=255)
                                                                                                                                     gCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAc               <  4:279495/149‑1 (MQ=255)
                                                                                                                                       aataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGAc                                                          <  3:169070/100‑1 (MQ=255)
                                                                                                                                       aataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGAc                                                          >  4:169070/5‑104 (MQ=255)
                                                                                                                                         taaTGCCTCCCGAGGCATTATTCTTATCATTGACGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAAGAAAAACATGc           >  4:610881/3‑149 (MQ=255)
                                                                                                                                         tGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAATCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATg            <  2:655149/148‑1 (MQ=255)
                                                                                                                                          aaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACt                               >  1:953225/2‑128 (MQ=255)
                                                                                                                                          aaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACt                               <  2:953225/127‑1 (MQ=255)
                                                                                                                                             gCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGTTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCgg                                          <  4:598414/114‑1 (MQ=255)
                                                                                                                                             gCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCgg                                          >  3:598414/1‑114 (MQ=255)
                                                                                                                                              ccTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTa      <  4:643674/149‑1 (MQ=255)
                                                                                                                                               cTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTAt     <  3:462663/149‑1 (MQ=255)
                                                                                                                                                tCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCa                                 >  4:376296/1‑120 (MQ=255)
                                                                                                                                                tCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATg    <  4:473024/149‑1 (MQ=255)
                                                                                                                                                 cGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATgg   <  4:618918/149‑1 (MQ=255)
                                                                                                                                                 cGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATg    >  1:112078/1‑148 (MQ=255)
                                                                                                                                                  gggAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATgg   >  2:682866/1‑148 (MQ=255)
                                                                                                                                                  gggAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGa  <  4:5595/149‑1 (MQ=255)
                                                                                                                                                  gggAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGa  <  3:861815/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGA  >  NC_000913/17381‑17675

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