Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,618 T→C Q114Q (CAA→CAG yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,6180TC98.0% 143.3 / ‑4.5 51Q114Q (CAA→CAGyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/21);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

ATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAA  >  NC_000913/7382‑7859
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
atcgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGGCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCAGCACTTTCTTTATAAAGCTGTGCGAGGGAACa                                                                                                                                                                                                                                        <  1:76969/248‑1 (MQ=255)
  cgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACAt                                                                                                                                                                                                                                       <  1:563648/249‑1 (MQ=255)
   gCCAGCAGAGCAAAGACAGCGAGAATAATACCTATCACCAGCGGGGGAAAATCAACAGAAAAACTCAGGGCTCGCGCAACGGCGATCGCCTGAATTCCGCTGAAAATTATGACGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATACCCAGCCCGCGCGCCATATACCATGCCGTTCCGACACGCAACTGCCCATTGGCATACCGATTTTTATAAAGCTGTGCGAGGGAACATTc                                                                                                                                                                                                                                     >  2:676572/1‑251 (MQ=255)
     cagcagAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCgg                                                                                                                                                                                                                                   >  1:274353/1‑251 (MQ=255)
           agcaaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGCCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCaaa                                                                                                                                                                                                                               >  1:283527/4‑249 (MQ=255)
             caaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACg                                                                                                                                                                                                                             >  1:444787/2‑249 (MQ=255)
              aaaGACAGCGAGAATAATAACTGTCAACAGCGGGGGAAAGTCAAAAGAAAAACTCAGGGATCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAAAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCAATTGGCATCACGTTATTTATAAAGCTGTGCGAGGGAACCTTCGGCAAACGAgg                                                                                                                                                                                                                          >  1:588138/1‑251 (MQ=255)
                    aGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGcc                                                                                                                                                                                                                     >  1:675250/1‑250 (MQ=255)
                        agaATAATACCTGTCAACAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGAGCAACGGCGGTCGCCTGAATTCCGCTGAAAATTATGCAGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCAGCCACGACACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGCACATTCGGAAAACGAGGTCGcactgcc                                                                                                                                                                                                                >  2:468907/1‑245 (MQ=255)
                           ataataCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCACGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCCTATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGCGGTCGCCATTCCGa                                                                                                                                                                                                              >  2:330703/1‑250 (MQ=255)
                             aataccTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATa                                                                                                                                                                                                            >  1:617306/5‑250 (MQ=255)
                                accTGTCACCAGCGGGGGAAAATCAAAAGGAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAGTTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGATGAGCGAGGGAACATTCGGCAATCGAGGCCGCCATGCCGATAAAcg                                                                                                                                                                                                        >  2:46497/2‑251 (MQ=255)
                                   ggTCACCAGCGGGGGAAAAGCAAAAGAAAAACGCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCAGGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCgg                                                                                                                                                                                                     <  1:86498/250‑1 (MQ=255)
                                     tCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGc                                                                                                                                                                                                    >  1:140285/1‑250 (MQ=255)
                                     tCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGAAAACGAGGGCGCCCTGCCGATAAACGCGGc                                                                                                                                                                                                    >  2:508498/1‑250 (MQ=255)
                                             gggggAACATCAAACGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGCAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCCTATACCATGCCGGTCCGCCCCGCAACTGCCCCTGGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGACAAACGAGGACGCCCTGCCGATAAACGCGGCCACCcacac                                                                                                                                                                                           >  2:154775/1‑250 (MQ=255)
                                                  aaactCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCCCGCAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGCAGCCCTGCAGATAAACGCGTAAACCCACATcc                                                                                                                                                                                         >  2:195973/5‑248 (MQ=255)
                                                       cAAAAGAAAAACTCAGGGCACGCGAAACGGCGTTCGCTTGAAGTCCGCTGAAAATTATGCCGCAGGCGATGAGCAAAAAGACGGCGAAAAGAACCCCCATCCAGCGCACCCCCAGCCCGCGCGCCTTATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGAc                                                                                                                                                                                  <  2:58691/250‑1 (MQ=255)
                                                        aaaaGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGCGCAAAAGGACGGCGAACAGAACGCCCCTCCAGCGCATCACCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGACACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGTGGTCGCCATGCCGCTAAACGCAGCAACCCACATCCAGAAGACgg                                                                                                                                                                                >  1:320716/1‑251 (MQ=255)
                                                                             cccaccggcgtccgcctaaatcccGCTAAAAATTAGGCCGTGGGCGATGCGCAAAAAGCCGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGGGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTaa                                                                                                                                                           <  2:672594/229‑1 (MQ=255)
                                                                               cAACGGCGTTCGCCTGAATTCCGATGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAAATCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGGGGACGCCATGCAGTTAAACGCGGAAACCCACATCCAGACGACTGCCACAGGTCCACccgcggcaaat                                                                                                                                                         >  2:657629/1‑246 (MQ=255)
                                                                                aaCGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCACCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAttg                                                                                                                                                        <  1:334784/251‑3 (MQ=255)
                                                                                aaCGGCGTTCGCCTGAATCCCGCGGAAAATTATGCCGTAGGCGATGACCAAAAAGACGGCGAACAGAACCCCCATCCAGCGCTCCCCCAGCCCGCGCGCCATATCCCAGCCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAttg                                                                                                                                                        <  1:263093/251‑3 (MQ=255)
                                                                                      gatCGCCTGAATTCCGCTGAAAATTATGCCGAAGGCGATGCGCAAAAAGACGGCGAACAGCACGCCCAACCGGCGCATCCCCAGCCCGCGCGCCATATACCATGACGGGCCGCCACGAAACTTCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTAGGCAACAGAAGTAGCCCTGCCGAACAACGCGACAACACACATACAGAAGCAGGCTCCACATCCACcgccgggaaaatacagc                                                                                                                                                   >  2:376219/3‑239 (MQ=255)
                                                                                           cctGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTTCCCGCGCAAcc                                                                                                                                            >  1:382952/3‑250 (MQ=255)
                                                                                           cctGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGGAAtttcccacgcaacc                                                                                                                                             >  1:632394/3‑238 (MQ=255)
                                                                                              gaattCCGCCGAAAATTATTCCGTCGGCGATGATCAAAAAGCCGGCGAACAGAACCCCCATCCAGCTCATCCCCAGCCCGCGCGCCATATCCCAAGCCGGGCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCTAGCTCAACGCCg                                                                                                                                         <  1:166098/247‑1 (MQ=255)
                                                                                                           aaTTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGGTGCCGc                                                                                                                            >  1:67150/1‑251 (MQ=255)
                                                                                                            aTTATGCCGTAGGCAATGAGAAAAAAGCCGGCGAACAGACCGCCCAGCCAGCGCATCCCCAACCCGCGCGCCATATACCAGGCCGGTCCGCCACGAAACTGCCCAGTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGACGCCCTGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCt                                                                                                                           <  1:46497/251‑1 (MQ=255)
                                                                                                                         cGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGcgcc                                                                                                              <  1:154775/251‑1 (MQ=255)
                                                                                                                          gATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGACCAACGCCGGCCAGGGTGCCCCCACCCACGCGcgccc                                                                                                             >  1:284770/1‑250 (MQ=255)
                                                                                                                                                     ccccctCCAGTGCACCCCCAGCCCGCGCGCCATATCCCACTCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAACCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATg                                                                                 <  1:565219/246‑1 (MQ=255)
                                                                                                                                                          ccccgccaacccccgccccGCGCCCCATATCCCATCACGGTCCGCCAAGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGGCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAGTGCCAGCGCAACGCCGGCCAGGGTGCCGCGACCCACGCGCGTCGCAAGACTGGTACACAATGACTGAAATGAGGtt                                                                           <  2:544714/235‑1 (MQ=255)
                                                                                                                                                           ccAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCCGCGCAACGCCCGCCCGGGTGCCGCTACCCACGCGCGCCGCCATGCTGGGACACAATGaagatcggaagcgcc                                                                           >  1:480280/1‑237 (MQ=255)
                                                                                                                                                               caccacccccccccGCGCCCCATATCCCATGCGGGTCCGCAACAAAACTGCCCATTGGCACCACGTGCTTTATAAAGCTGTGGGAGGGACCATCCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACAGCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAAcc                                                                      <  2:629151/240‑1 (MQ=255)
                                                                                                                                                                gcATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCCACGCAACGCCGGCCCGGGTGCCGCTAACCACGCGCGCCGCCAGACTGGGACACAAAGACTGAAATGAGGTTaaaacc                                                                     >  1:20508/1‑248 (MQ=255)
                                                                                                                                                                       cAGCCCGCGCGCCATATCCCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGTTACACAATGACTGAAATGAGGTTAAACCGCCTGGCt                                                              <  1:571211/251‑1 (MQ=255)
                                                                                                                                                                         ccccgcgccccaaaacacATGCCGGTCCACCACAAAGCTGCCCGTTGGCATCAGTTGCTTTATAAGGCGGTGCGGGGGAACATCCGGAAAACGAGGTCGCCATGCCGATAAACGCGGCACCCCACATCCAGAAGACGGCACCAGGTCAACCGGCGGTAATTGCCAGCGTAACGCCGGCCAGGTGGCCGCTCCCCACGCGCGCCGAAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCtgt                                                            <  2:107214/234‑1 (MQ=255)
                                                                                                                                                                                  ccATATCCCATGCCGGCCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGAAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATg                                                   <  1:603899/251‑1 (MQ=255)
                                                                                                                                                                                   cATAGACCATGCCGGGCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATGCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGAc                                                                                                       <  1:86826/200‑1 (MQ=255)
                                                                                                                                                                                      ataCCATGCCGTCCCCCCACGAACCTGCCCATGGGCATCACGTTCTTAATAAACCGGTGCGAGGTACCATCCGCCAAACGAGGTCGCCTTGCCGATAAACGCGGCAACCCACACCCAGAATACGGCTCCAGGTCCACCGGCGGTAATGGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCCCGCGCGCCGTAAGACTTGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTAt                                               <  1:340879/251‑1 (MQ=255)
                                                                                                                                                                                          cATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGAAAATTCCAGCGCCAACCCGGACAGGGTGCCGCCACCCCCACGGCCCGCAAGGCTGGGCAACACTGAATGAAATGAGGGTAAACCCCCTGGCTGTGGATGAATGCTCttttt                                           >  2:415830/1‑251 (MQ=255)
                                                                                                                                                                                               cGGTCCGCCACGAACCTGCCCATTGGCATCCCTTCCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGAc                                      <  1:144357/251‑1 (MQ=255)
                                                                                                                                                                                               cGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTa                                                                          <  1:668087/215‑1 (MQ=255)
                                                                                                                                                                                                    cGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGAAATTGCCAGCGCAAAGCCGGCCAGGATGCCGCTACCCACGCGCGCCGCCAGACTGGTACACCATGACTGGAATGAGGGTAAACCGCCTGGCTTTGGGTGTATGCTATTTTTAAggcctttg                                 >  1:480929/1‑244 (MQ=255)
                                                                                                                                                                                                     gCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCCTCGCAACGCCGGCCAGGGTGCCGCCACCCCCCCGCGCCGCATGACTGGGACACAATGGCTGAACAGAGGTTAAACCGCACGGCTTGGGATGAATGATCTTTTTAAAACtttttc                                >  2:526969/1‑249 (MQ=255)
                                                                                                                                                                                                       cACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCCGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGGACACCATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGGTGAATGCTATTTTTAAGACTTTTGccc                              >  1:364508/1‑250 (MQ=255)
                                                                                                                                                                                                          caaaaccgcccattggcatcccttacTTTATAAAGCTGTGCGAGGGAACATGCGGCAAATGGGGTCGCCAGGCCGAAAAACGGGGCAACCCACACCCAGAAGACGGCTCCTGGTCCACGGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAGGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCaaa                            <  2:675250/226‑1 (MQ=255)
                                                                                                                                                                                                             aCTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGGAATTGCCAGCGCAAAGCCGGCCAGGGTGCCGCTACCCCCGCGCGCCGCAAGGCTGGTACACAATGACTGAAATGAGGTTAAACCGCCGGGGTGTGGGAGAATGCTATTTTTAAACCTTTTGCCAAACTgg                        >  1:227743/1‑251 (MQ=255)
                                                                                                                                                                                                                              cGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCCGCCCAACGCCGGCCAGGTTGCCGCTACCCCCGCTCGCAGCAAGACTGGACAACAATGCTTGAAATTGGGTTTAACCGCCTCGATGTGGGTGAATGCTCATTTTAACACTTTTTCCAAACTGGCggggctagaggaaaa        >  2:7721/1‑237 (MQ=255)
                                                                                                                                                                                                                                     ttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTTATAGCCAGCCCCACCCCGGCCAGGGTGCCGCCACCCACCCGCGCCGCCAGAACGGGACACAATAGCTGCAATTAGGGTTAAACCCCCGGCTGGGGAGGAAGGCGATTTTTAAAAAATTTGACCAACTGGCgggggggaaggaaacacacac  >  2:120823/1‑231 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
ATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAA  >  NC_000913/7382‑7859

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