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,7930CT95.3% 106.1 / ‑2.3 44R89R (CGC→CGTinsA7IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (20/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

GATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATG  >  NC_000913/4518546‑4519039
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
gATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCTCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGt                                                                                                                                                                                                                                                     >  1:89836/1‑251 (MQ=255)
 aTGTCCTTCCTGCCCCGCTACTGAAGGGGTGGTGCGTAACGGCAAAAGCCCTGCCGGCCATCAGCGCTATCTCTGCCCTCATTGCCGTAAAACAGGGCAACTACAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTc                                                                                                                                                                                                                                                    <  1:196850/251‑1 (MQ=17)
  gttccttcctGCCCCGCTCCTGAAGGGGTGGTGGTTAACGGCAAAACCCCTGCGGACCACCAGCGCCTCCCCTCCTCTATTGGCCTTAAAACATGGCAACTACCGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCg                                                                                                                                                                                                                                                   <  2:87401/249‑1 (MQ=17)
                    ggcgAAGGCGTGGTGCGTAGCGTCAAAATCTCGGCCTAACATCAGCGCTACCTCTGTTCCTTTTGCCGTAAAACATGGCAATATCATTTCACTTACACCGCCTCTCAGCCCGGTACGCGCCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACTCGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACaa                                                                                                                                                                                                                                 <  2:219962/248‑1 (MQ=17)
                    aCTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTGGGCCTCAACACGGTTTTACGTCACTTCAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACa                                                                                                                                                                                                                                  >  1:175323/1‑250 (MQ=255)
                     cTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTAATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACAAGGTTTTCCGTCACTTAAAAAACTCAGGCCGGAGCCGGTAACCTCGCGCATACAAc                                                                                                                                                                                                                                >  2:328472/1‑251 (MQ=17)
                       acaGGCGTGGGGCGTAACGTCAAAACTCGTGCCGGACATCGGCGCTATCCCTGCCCTCATTGCCTAAAAACTTGGAACCTCCGTGTCCCTGCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCg                                                                                                                                                                                                                              <  1:105729/249‑1 (MQ=17)
                         aGGCGTGGTACGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGGCCTCAACACGGATTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCggg                                                                                                                                                                                                                            >  2:103567/1‑251 (MQ=255)
                            cGTGGCGCGTAACGGCAAAAGCCAGGCCGTACATCAGCGCTATCTCTGCGCCCATGGCCGAAAAACAGGGCAACAACAGTTCATTGACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTGAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAg                                                                                                                                                                                                                         <  2:185920/251‑1 (MQ=17)
                                           cAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTg                                                                                                                                                                                                          >  1:162458/1‑251 (MQ=255)
                                            aaaaGCACTGCCGGACACCAGCGCTATCTCGGCTCTCATTGCCTTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGTCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACTGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgc                                                                                                                                                                                                         <  1:16594/251‑1 (MQ=18)
                                                      gCCGGACATCAGCGCTACCTCTGCTCCCATTGCCGTAAAACAGGGCAACTACAGTCCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGaaa                                                                                                                                                                                                  <  1:73632/249‑1 (MQ=18)
                                                                    gcTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTggg                                                                                                                                                                                   >  1:281601/1‑251 (MQ=255)
                                                                                ctcattGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACCGTGGGGCTACGTCGGTg                                                                                                                                                                       >  1:319473/5‑251 (MQ=255)
                                                                                         atAAAACATGGCAACTCCAGCTCTATTTAGCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGTCATGAATGGCGTCGGATTGCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCAc                                                                                                                                                              <  2:119323/250‑1 (MQ=18)
                                                                                             aaCATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATTGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCACTGGACGCACAGTGGGGCTACGTCGGTGCTAAATCACGTc                                                                                                                                                           >  1:255578/1‑250 (MQ=255)
                                                                                                  tgCAACTACATTTCACTTACACCGCCTCTCAGCCCGGTACGCACAAGAAAAACATTGATATGGCCCTGAATAGAGTCGGCTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTAAACCCGGTTTTACGTCACGTAAAAAACTCAGGCCGTAGTCGGTACCCTCGCGCATACAACCAGACAGTGATGTTATTGTCTGCGCTGAAATGGGCGAACCGTGGGGCTACGTCGGTGCTAAAGCACGTAAGCgct                                                                                                                                                     >  2:377061/2‑251 (MQ=255)
                                                                                                         aCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAACGCTGGCTGtt                                                                                                                                              >  1:215413/1‑251 (MQ=255)
                                                                                                           cAGTTCATTTACCCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCTCGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGttt                                                                                                                                             <  1:280505/251‑1 (MQ=255)
                                                                                                                tCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGCCCCCCCAGTGCCCGCCTTCTTGGCGTTGGCCTCAACACGGTTTTACGTCACTTACAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGCGCTAAATCACGTCAGCGCTGGCTGTTTTAcg                                                                                                                                         >  1:53999/1‑251 (MQ=17)
                                                                                                                   aCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCCTTGATATGGCCATGACTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCATCAACACGGTTTTACGTCCCTTAAAAAACTCAGGCCGTAGACGGTAACCTCGCGCATACAACCGGGAAGTGATGTGCTTGTCTGCGCTGAAATGGACGAACAGTGGGGCGTCGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTAcggg                                                                                                                                       >  2:143875/1‑249 (MQ=255)
                                                                                                                     ttACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGt                                                                                                                                      >  1:150524/1‑250 (MQ=255)
                                                                                                                              ctctCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGACATGGACGAACCGTAGGGCTACGTCGGGGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAgg                                                                                                                              >  2:74220/1‑249 (MQ=255)
                                                                                                                                  cagcacggtacgaacacgaacaTCATTGAGGTGGCCATGAATGGCGTCGGATGTCGCGCCAGTGGACGCATTATGGGCGTTGGCCTCAACACGGTTTTAAGTCACTTAAAAAACTCAGGCAGTAGGCGGTAACCTCGCGCATACACACGGTCCGTGATGTGAGTGTCTGCGCTGAAATGGACGAACAGTGGTGCAACGTAGGTGCAAAATCACATCCGAGCAGGCTGTATTAAGCGCATGACGGGATACgg                                                                                                                        >  2:76845/22‑251 (MQ=255)
                                                                                                                                                       aTCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTc                                                                                                   >  1:34138/1‑251 (MQ=255)
                                                                                                                                                             gATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAAACGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTCAATCACGTAAGCGCTGGCTGTTTTACGCGTCTGACAGGAAGCGGAGGACGGTTGTGGCGCAAGTCTTCGGt                                                                                             >  2:310504/1‑251 (MQ=255)
                                                                                                                                                                                       tCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTGAGCGTATGACAGGATCCGGCGGACGGTTGTGGAGCACGTCTTCGGTGGACGCACTCTGGCCGAACTGGAGCg                                                                   >  2:98194/1‑251 (MQ=255)
                                                                                                                                                                                        cgcgCCAGTGCACGCATTATGGGCGTTGGCCGCAACACGGTTTTACGTCACTTAAAAACCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCGGGCCACACTGGAGCGt                                                                  <  2:150524/251‑1 (MQ=255)
                                                                                                                                                                                                     gCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATgg                                                                                                                                                                                               <  1:27336/113‑1 (MQ=255)
                                                                                                                                                                                                     gCATTATGGGCGTTGGACTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATgg                                                                                                                                                                                               >  2:27336/1‑113 (MQ=255)
                                                                                                                                                                                                                  tGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGGGAACGCACTCTGGCCACACTGGAGCTGCTTCTGAGCCTTCTGTCGGCCTTTgg                                        >  2:144519/1‑250 (MQ=255)
                                                                                                                                                                                                                    gCCTCAACACCGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGGCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACTCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAgg                                      <  2:162458/251‑1 (MQ=255)
                                                                                                                                                                                                                          acctggTTTTACGTCCCTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTTCGCGGAAATGGACGAACAGTGGGGCTACGGCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTGCGGTGAACGCACTCGGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTgg                                <  2:300932/247‑1 (MQ=255)
                                                                                                                                                                                                                           cacGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGt                               <  1:103567/251‑1 (MQ=255)
                                                                                                                                                                                                                           cacGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCGGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTCACCCTTTGacgtggtgat                               <  1:278305/251‑10 (MQ=11)
                                                                                                                                                                                                                            acGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGta                              <  1:144519/251‑1 (MQ=255)
                                                                                                                                                                                                                              tgTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGAGATATGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGGGCTAAATCACGTCAGAGCTGGCTGTTGTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTATTCGGCGGACGCACTCGGGCCACACCGGAGAGTCTTCTGAGCCTGTTGGCGGCCTTTGAGGTCGAGGtgt                             >  2:291770/2‑248 (MQ=255)
                                                                                                                                                                                                                              ggTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGAt                                                                                                                            <  1:71821/155‑1 (MQ=255)
                                                                                                                                                                                                                              ggTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtata                            <  1:215496/251‑1 (MQ=255)
                                                                                                                                                                                                                                     catCACTTAAAAAACTCAGGCCGTAGTCGGAAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCAATGGACGAACAGTGGGGCTACGTAGGTGCTAAATCACGTCCGCCATGGCTGTTTTACGCGTAAGACCGGATACGGGGGACGGTTTTGGCGCACGTCTTCGGGGCACGCAATCTGGCCCCACGGGAGCGTCTTCCGAGCCTGCTGTCGGCCTTTgcggctgtgggttcgatgc                       >  2:139741/3‑231 (MQ=255)
                                                                                                                                                                                                                                            tAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGAtggc               <  1:98194/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aaaaaaCTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGAtggct              <  1:291770/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     cAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGCCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCt       <  1:139741/249‑1 (MQ=255)
                                                                                                                                                                                                                                                        gCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGACTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATg  <  1:118840/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 < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: