Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,681 A→C K52T (AAG→ACG)  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,6810AC98.4% 183.0 / NA 63K52T (AAG→ACG) insA7IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (35/27);  total (36/27)
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.92e-01

TCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGG  >  NC_000913/4518443‑4518921
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
tCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAAt                                                                                                                                                                                                                                      >  1:337143/1‑251 (MQ=255)
 ccTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAAt                                                                                                                                                                                                                                      <  2:155763/250‑1 (MQ=255)
         cccTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGatat                                                                                                                                                                                                                             >  1:488526/1‑251 (MQ=255)
         cccTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGatat                                                                                                                                                                                                                             >  2:378346/1‑251 (MQ=255)
               aGGATGGGGATTTCGGTAATGCTGCCAACTTACTGCTTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCAt                                                                                                                                                                                                                       >  1:113946/1‑251 (MQ=255)
                                    gctcttccgatcTGATTTAGTGGATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACGCCGCCTCTCAGCCCGGTACGCACCAGAAAATAAATGATATGGCCATGAATGGCGTCGGATGTcaggc                                                                                                                                                                                                  >  2:130900/12‑247 (MQ=18)
                                                  aTTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGCAAATCATTGATATGGCCATGAATGGCGTCGGATGGCGCGCCAGAGCACGCAtt                                                                                                                                                                                     >  2:277702/1‑250 (MQ=255)
                                                                ggTGATTTTAAGGTGCTTGCGTGTCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATTAATGTCGTTGGATGCCGGGCAACCGCCCGCATTATGGGCGTTGGCCt                                                                                                                                                                       <  1:142990/250‑1 (MQ=11)
                                                                        tAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGt                                                                                                                                                              >  2:418207/1‑251 (MQ=255)
                                                                              gCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCATCAACACGGTTTTACg                                                                                                                                                        >  1:461965/1‑251 (MQ=255)
                                                                                 tGCGTGGCTTCCATTTCCATCAGATTTCCTTCCTGCCCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACAGGGCAACTCCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGt                                                                                                                                                       <  1:25033/249‑1 (MQ=255)
                                                                                   cGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACt                                                                                                                                                   >  2:149074/1‑251 (MQ=255)
                                                                                     tGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACtt                                                                                                                                                  >  2:411366/1‑250 (MQ=255)
                                                                                     tGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGGTTTACGTCACtt                                                                                                                                                  >  2:207713/1‑250 (MQ=255)
                                                                                                tccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                                                                                      >  1:179293/1‑251 (MQ=255)
                                                                                                  catCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                                                                                      >  2:124520/1‑249 (MQ=255)
                                                                                                  catCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                                                                                      <  1:124520/249‑1 (MQ=255)
                                                                                                       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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
TCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGG  >  NC_000913/4518443‑4518921

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: