Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,793 C→T R89R (CGC→CGT insA7 → IS1 family repressor protein InsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,7930CT98.6% 211.2 / NA 73R89R (CGC→CGTinsA7IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/33);  total (39/34)
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.81e-01

GATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATG  >  NC_000913/4518546‑4519039
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
gctGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGACCATGAATGGCGTCGGATTTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAACAAACTCAGGCCGTAGt                                                                                                                                                                                                                                                      >  2:487692/3‑251 (MQ=18)
  tGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGAGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGCTATGGCCATGAATGGCGTCGGATGTCGCGCCCGTGCACGCATTATGGGCGTTGGCATCAACCCGGTTTTAAGTCACTTAACAAACTCAGGCCGTAGTCg                                                                                                                                                                                                                                                    >  2:400151/1‑251 (MQ=18)
      cttcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGCCATCAGCGCTATCTCTGCTCTCATAGCCGTAAACCATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATAATTGATATGGCCATGAATGGCGTCGTATGTCGCGCCAGTGCACGCCTTATGGGCGTTGGCCTCAACACGGTTTTACGTCCCTTAAAAAACTCAGGCCGTAGTCGGTa                                                                                                                                                                                                                                                 >  1:474378/1‑250 (MQ=255)
      cttcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCCGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCCTTGTTATGGCCATGAATGGCGTCGGATGTCGCGCCAGAGCCCGCATTCTGGGCGTTGGCCTCGACACGGTATTACGTCACTTAAAACACTCAGGCCGTAGACGGTa                                                                                                                                                                                                                                                 >  2:272014/1‑250 (MQ=17)
       ttcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACACCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTa                                                                                                                                                                                                                                                 <  1:393799/249‑1 (MQ=255)
           tGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTcg                                                                                                                                                                                                                                           >  1:304860/1‑251 (MQ=255)
                       gagggcgtgtgtcgaaaCGACAAAAGCCCTCCCGAACACCCGCGCTATCTCTGCTCTCATTGCCGAAAAACAGGGCAACTCAAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAAc                                                                                                                                                                                                                                 <  2:285727/234‑1 (MQ=17)
                       gAAGGCGTGGTTCGTAACGGCAAAAGCACTGCCGCAAACCAGCGCTATCTCTCCTCTCACCGCCGTAAAACTTGGCAACTACAGTTCACTTATCCCGCCTCTCAGCCCGCTACGCACCAGAAAATCATTGATATGTCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCg                                                                                                                                                                                                                               <  1:143141/251‑1 (MQ=255)
                          ggCGTGGTGCTAAACGAAAAAATCACTCCCGGACATCCGCGCTATCTCTGCTCCCATTGCCGTAAAACGGGGCAAATACAGTTCACTTACCCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCTCTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGc                                                                                                                                                                                                                            <  2:129714/251‑1 (MQ=17)
                                gtgCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTAGCGCCAGTGCACGCATTCTGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAAATCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgat                                                                                                                                                                                                                      >  2:461387/1‑251 (MQ=255)
                                   cGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTCCGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatg                                                                                                                                                                                                                     >  1:207950/1‑249 (MQ=255)
                                         ggCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTc                                                                                                                                                                                                             >  2:459227/1‑251 (MQ=255)
                                          gCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACGACAGTTCACTTACACCGCCTCTCAGCCCGGGACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGACAGTGCACGCATTATGTGCGTTGGCCTCAACACGGTTTTACGACACTTAAAACACTCAGGCCGTAGTCGGTAACCTCGCACATACCACCGAGCAGGGATATGATTGTCt                                                                                                                                                                                                            >  1:149026/1‑251 (MQ=17)
                                            aaaaGCACTGCCGGAAATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgc                                                                                                                                                                                                          <  1:445878/251‑1 (MQ=255)
                                             aaaGCACTGCCGAAGATCAGCCCTTGCCCTGCTCTCATGGCCGTAAAACATGGCACCTACAGTTCCCTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgcg                                                                                                                                                                                                         <  2:180956/251‑1 (MQ=18)
                                                        gcaCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACAACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCAACAGAAAATAATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACAAGGTTTTACGTCACTTCAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAGTGGAc                                                                                                                                                                                              >  2:357526/3‑251 (MQ=255)
                                                               aGCGCTATCTCTGCTCCCATTGCCGAAAAACATGGCAACTACATTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGt                                                                                                                                                                                       <  1:300614/251‑1 (MQ=18)
                                                                    tATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGc                                                                                                                                                                                  >  1:397448/1‑251 (MQ=255)
                                                                      cctctGCTCTCATTGCCGAAAAACATGGCAACTCCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTa                                                                                                                                                                                <  1:409456/250‑1 (MQ=255)
                                                                         ctGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTa                                                                                                                                                                                <  1:174137/248‑1 (MQ=255)
                                                                           cctctcatggCCGTAAAACAGGGCAACTACAGTCCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGTGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTc                                                                                                                                                                            <  1:272014/241‑1 (MQ=18)
                                                                                cattGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCAATGGCCGAACAGTGGGGCTACGTCGGTGCt                                                                                                                                                                      >  2:59822/3‑251 (MQ=255)
                                                                                       gTAAAACATGGCACCTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCAc                                                                                                                                                               <  1:186652/251‑1 (MQ=255)
                                                                                           aacatgacacatactgTTCTCTTACATCGCCTCTCAGCCCGATACGCACAAGTAAATCATTGATAGTGCCATGAATGGCGTCGGATGGCGCGCCCGTGCACGCATTATGTGCGTTGTCCTCAACACGGTTTTACGTCACTTAACAAACGCATGCCGGAGTCGGTAACCGCGCGCATACCTCCGTGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGTTGCTAAATCAcctcc                                                                                                                                                           <  2:57091/236‑5 (MQ=255)
                                                                                                gtaaactccatttccatacaacCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAACGGCGTCGGATGTCGCGTAAGTGCACGCATTATGGGCGTTGTCCCCAACACGGTTTTACGTCACTTAAAAAACTCAGGTCGTAGTCGGTAACCTCGCTCATACTACCGGGCCGTTATGTGATTGTCTGCGCTGAAATGGACGACCAGTGGGGCTACGTCTGTCCTAAATCACGCCAgcgcc                                                                                                                                                      <  2:452838/231‑2 (MQ=17)
                                                                                                ggCACCTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCgct                                                                                                                                                      <  1:104692/251‑1 (MQ=255)
                                                                                                ggCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATACCGTCAGCgct                                                                                                                                                      >  2:390688/1‑251 (MQ=255)
                                                                                                  caattacgttcgctacaccccGCCTCTCAGCCCTGTACGCACCAGAAAATCATGGATATGGCCATGAATTGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTATGTCGGGGCTAAATCGCGCCCGCgctgg                                                                                                                                                    <  2:484040/232‑1 (MQ=17)
                                                                                                    aataaaGTTCACTTACACCGCCTCCCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCgctggct                                                                                                                                                  <  1:105382/246‑1 (MQ=255)
                                                                                                      tACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTAGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGt                                                                                                                                                >  1:205320/1‑251 (MQ=255)
                                                                                                       aCAGTTCACTTACCCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGtt                                                                                                                                               <  1:57025/251‑1 (MQ=255)
                                                                                                       aCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatgtga                                                                                                                                                                                                                  >  1:367776/1‑184 (MQ=255)
                                                                                                         aGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGCTTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGtttt                                                                                                                                             >  2:203597/1‑251 (MQ=255)
                                                                                                            tCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTAc                                                                                                                                           <  1:418207/250‑1 (MQ=255)
                                                                                                               cTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGt                                                                                                                                       >  1:121035/1‑251 (MQ=255)
                                                                                                                 tACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTAt                                                                                                                                     >  1:491237/1‑251 (MQ=255)
                                                                                                                   cacCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATg                                                                                                                                    >  1:327777/1‑250 (MQ=255)
                                                                                                                    acCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGAc                                                                                                                                  >  1:374342/1‑251 (MQ=255)
                                                                                                                     ccGCCTCTCAGCCCGTGACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCTTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGTCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACa                                                                                                                                 <  1:335866/251‑1 (MQ=255)
                                                                                                                     ccGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGGGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACa                                                                                                                                 <  2:113946/251‑1 (MQ=255)
                                                                                                                        ccTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGa                                                                                                                              <  2:501433/251‑1 (MQ=255)
                                                                                                                         ctctCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACCGGAt                                                                                                                             >  2:453519/1‑251 (MQ=255)
                                                                                                                         ctctCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGAt                                                                                                                             >  2:192907/1‑251 (MQ=255)
                                                                                                                         ctctCAGCCCGGTAAGCACCAGAAAATCATTGATCTGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCTCGTCAGCGCTGGCTGTTTTTCGCGTATGACAGGAt                                                                                                                             >  2:69524/1‑251 (MQ=255)
                                                                                                                             cAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACgg                                                                                                                         <  1:210797/251‑1 (MQ=255)
                                                                                                                             cAGCCCGGTACGCACCAGAAAATCATTGATATGGACATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATa                                                                                                                            <  2:38325/248‑1 (MQ=255)
                                                                                                                                      acgCACCAGAATATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTGTGGGCTTTTGCCTCAACACGGTTTTACTTCAATTAAAAAACTCAGTCCGTAGTCGGTAACCTCGCGCCTACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGTGGCTACGTCGGTGCTAAATCTCGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGtt                                                                                                               <  2:169204/249‑1 (MQ=18)
                                                                                                                                       cgccccAGAAAAACATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAATCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTtg                                                                                                              <  1:275872/247‑1 (MQ=255)
                                                                                                                                        gCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTtgt                                                                                                             <  1:459227/251‑1 (MQ=255)
                                                                                                                                         cACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTtgtg                                                                                                            >  1:45449/1‑251 (MQ=255)
                                                                                                                                                 aaTCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCCAAATCACGTCAGCGCTGGCTGTTTTACGCGTACGACAGGATACGGAGGACGGTTGTGGCGCACGt                                                                                                    >  1:408729/1‑251 (MQ=255)
                                                                                                                                                           gtGGCCATGAATGGCGTCGGCTGTCGAGCCAGTGTACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGAGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCAATGGACGAACAGTGGGGCTACGGCGGGGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGGATGACAGGGTACGGAGGACGGTTGGGGCGCACGGCTTCGGTGa                                                                                           >  2:106686/2‑250 (MQ=255)
                                                                                                                                                                 aTGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCAct                                                                                    <  2:374342/251‑1 (MQ=255)
                                                                                                                                                                    aaTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGGGGACGGTTGTGGCGCACGTCTTCGGTGAACGCAccc                                                                                   >  2:237312/1‑247 (MQ=255)
                                                                                                                                                                     aTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGACATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGCTACGGAGACAGGCTGTGGCGACCTTCTTCGGTGAACGCCctcttg                                                                                >  2:416175/1‑249 (MQ=21)
                                                                                                                                                                         cGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGACCGCACTCTGGCcac                                                                            >  2:185537/1‑251 (MQ=255)
                                                                                                                                                                              gatgtccctcccgtgccctcATTATGGGCGTTGGCCTCAACACGGTTTTCCGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATTCAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATTGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTgg                                                                       <  2:179293/232‑1 (MQ=18)
                                                                                                                                                                                 gTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAg                                                                     <  1:411366/250‑1 (MQ=255)
                                                                                                                                                                                      gcCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACCGGATACGGAGGACGGTTGTGGCGCCCGTCTTCGGTGAACGCACTCTGGCCACCCTGGAGAGtctt                                                               >  2:128950/1‑251 (MQ=255)
                                                                                                                                                                                        cAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGtcttct                                                             <  1:277702/251‑1 (MQ=255)
                                                                                                                                                                                          gTGCACGAATTCTGGGCGTTGGCCTCACCACTGTTTTACGTCACTTAAACAACTCAGGCCGTAGTCGGTAACCTCACGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGAAGAACAGTGGGGCTACGTCGGTTCTAAATAAAGTCAGCGCTGGCTGTTTTCCGCGTATGCAAGGTAACGGAGGAAGGTTGCGGGCCAAGTCTTAGGTGAAAGACCTCTAGCCACCCTGGAGCAtcttatgg                                                           >  2:339594/1‑247 (MQ=255)
                                                                                                                                                                                           tGCACGCATTATGGGCGTTGGCCTCAACACGGCTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGGGATGTGATTGTCTGCGCTGAAATGGACGAGAAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTATGCGTATGACAGGATACGGAGGGCGGTTGGGGCGAACGACATCGGAGAAGGCACTCTGGCCACCCTggcacgacgtctg                                                            >  2:1422/1‑238 (MQ=255)
                                                                                                                                                                                            gCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGc                                                         <  1:285932/251‑1 (MQ=255)
                                                                                                                                                                                             cACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCCCGTCGTCGGAGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGcc                                                        >  2:434199/1‑251 (MQ=255)
                                                                                                                                                                                              aCGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAACCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCAAGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCct                                                       >  1:376345/1‑251 (MQ=255)
                                                                                                                                                                                                         gCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCCGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGAAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGGGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGGGGCGCACGTCTTCGGTGAACGCACTCTGGCCACAATGTGGCGTCTTCTGGGCCTGCTGTCGGCCt                                            >  2:340989/1‑251 (MQ=21)
                                                                                                                                                                                                          cGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCtt                                           >  1:368003/1‑251 (MQ=255)
                                                                                                                                                                                                              ggCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGa                                        <  1:360936/250‑1 (MQ=255)
                                                                                                                                                                                                                ccTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGAAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTCTGACAGGATACGGCGGACGGTTGTGGCGCACGTCTTCGGTGACCGCACTCTGGCCACACTGGAGCGTCTTCTGTGCCTGCTGTCGGCCTTTGAggg                                     >  2:49093/1‑250 (MQ=255)
                                                                                                                                                                                                                      cacGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGt                               <  1:453519/251‑1 (MQ=255)
                                                                                                                                                                                                                                        acaaaaCTCAGGCAGTAGTCGGTAACCTCGCGCATACAACCGGGAAGAGTTGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGACCGCACTATGGCCACACTGGAGCGTCTTCTGAACCTGCTGTCGGCCTTTTAGGACGTGCTATGGATTGCGAAtggctg             >  2:131512/3‑251 (MQ=18)
                                                                                                                                                                                                                                          aaaaCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGAtggctggc           >  1:217101/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   gCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATg  <  1:149074/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATG  >  NC_000913/4518546‑4519039

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