Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,635 C→T 8.1% A353T (GCC→ACC)  tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,971,6350CT 8.1% 161.5 / 2.3 59A353T (GCC→ACC) tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/3);  total (28/31)
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.89e-01

GTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACG  >  NC_000913/1971340‑1971933
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                          
gtcaCGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGaacc                                                                                                                                                                                                                                                                                                       >  2:124427/1‑298 (MQ=255)
  caCGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGAGAACTCa                                                                                                                                                                                                                                                                                                     >  2:525125/1‑301 (MQ=255)
                        atatTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCaa                                                                                                                                                                                                                                                                               >  2:201136/1‑300 (MQ=255)
                               tCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATccaccaa                                                                                                                                                                                                                                                                        >  2:336729/1‑300 (MQ=255)
                                cATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCAATCCACCACt                                                                                                                                                                                                                                                                       >  2:384472/1‑301 (MQ=255)
                                   gttttCCCGTGCGTTTTCGCTCGGTACCGTACCGGTATCAACGCGTGAGACGGGGTCTTCATTGATGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTg                                                                                                                                                                                                                                                                    <  1:232995/299‑1 (MQ=255)
                                    ttttcgccgggtcttgtggccaggtcctaaCCGGTATCAAGTCGTGGTACGGAGTCTTTAATGAGGGCTTTGATCTCTTTGGCCGCCGGGGCGCGGCGACGGGTAAGATTACGTACTGCCCCCGCTACCACGGCAAAACCACTGCCCTGTTCACCCGCACGCGCGGCTTCAACCGAGGCATTCAGCGCGAGGATATGAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTgc                                                                                                                                                                                                                                                                   <  1:444599/273‑1 (MQ=255)
                                         cctggtgcgttCGACTAGCACCGAACGTGTATCAACGCGTGAGACGGAGTCTTTAATGAGGGCTTTGATCTCGTTTGCCGCCTGGGCGCGGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGccgccg                                                                                                                                                                                                                                                              <  1:180585/292‑1 (MQ=255)
                                            ggCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGACGCCGTGc                                                                                                                                                                                                                                                           >  1:411984/1‑301 (MQ=255)
                                                 tttCGACCAGCACCGAACCGGTATCAAGGGTGGAGACTGAGTCTTCAATGGGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGc                                                                                                                                                                                                                                                      <  1:200893/301‑1 (MQ=255)
                                                        cAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCg                                                                                                                                                                                                                                               >  2:592019/1‑301 (MQ=255)
                                                           cACCGAACCGGTATCAACGCGTGTGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAg                                                                                                                                                                                                                                            <  2:128177/301‑1 (MQ=255)
                                                                          aaCGCGGGGGACGGAGTCTTTAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGt                                                                                                                                                                                                                             <  1:194422/301‑1 (MQ=255)
                                                                              cgTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCCGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATCATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCCTCGTTTTCACTACCCCATCCCCCACTTTGCCGCCGTGCTGGCCGGTATCGGAGGCCCTTTGCGCCCGTTGcc                                                                                                                                                                                                                         >  1:482766/1‑300 (MQ=255)
                                                                                         ttCTTCAATGAGGGCTTTTATCTCTTTTGCCGCCAGGGCGCTGCGACTGGCAAGATTACGCACTTCCACCGCCACCATGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCAGTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACGTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAggcctcgcaactggcgc                                                                                                                                                                                                              <  2:630739/300‑17 (MQ=255)
                                                                                                   ggggTTTTGCTCTCTTTTGCCGCCTGGGCGCTGCGACTGGCACTATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCgg                                                                                                                                                                                                    <  1:524453/300‑1 (MQ=255)
                                                                                                         ttggtctctgttgtgtcctgtgcgtgtcGACTGGCAAGATTACGCACTTCACCCTCTATGACGGTAAAACCACGGCCCTGTTCGCCCGCACGCGCTGCTTCTTCCGCGGCAGTGTGAGCTGTGATTATAGGCTGGATGGTTGTATTTTCGATAACGTTGATAATGTCGGCAATCTTCTTCGAACTATTTGTGATCGCATGAATCGTTTTCACTACGGCAGTCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGGACTTGGCGCCAGTTGCGAGGGCTGGCGGGCGTTATCGGCGgttt                                                                                                                                                                                              <  1:191526/274‑1 (MQ=255)
                                                                                                         ttGATTTCTTTTGCCGCCTGGTCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGtttt                                                                                                                                                                                              <  2:228608/301‑1 (MQ=255)
                                                                                                           gATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGtttacg                                                                                                                                                                                            >  1:228365/1‑298 (MQ=255)
                                                                                                           gATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGc                                                                                                                                                                                            >  2:276151/1‑301 (MQ=255)
                                                                                                                ctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCAc                                                                                                                                                                                       >  2:125536/1‑301 (MQ=255)
                                                                                                                     gccgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTcg                                                                                                                                                                                  >  2:55789/1‑301 (MQ=255)
                                                                                                                        gccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCgg                                                                                                                                                                               <  2:291333/301‑1 (MQ=255)
                                                                                                                           ggggctgtGCGTCTTGCATGATTGCGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGa                                                                                                                                                                            <  1:222091/294‑1 (MQ=255)
                                                                                                                                ccgggtattGGCAAGATTACGCACTTCACCCGCCACCGCGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAgctgc                                                                                                                                                                       <  1:217583/293‑1 (MQ=255)
                                                                                                                                 cTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGTGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTATCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCCGCCGTTTTCTTCAccttcccgcgtatgcgcc                                                                                                                                                                      >  1:531953/1‑284 (MQ=255)
                                                                                                                                   gCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCCCTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCCCTTTGCGCCAGTTGCGCGGCCTGGAGGGCGTTAACGGCGTTTTTCTTCCCTGTCCCGGttctcctgccc                                                                                                                                                                    >  1:326607/1‑291 (MQ=255)
                                                                                                                                       cTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGc                                                                                                                                                                >  2:458279/1‑301 (MQ=255)
                                                                                                                                           cAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCTCCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGGTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCCCCATGCCggc                                                                                                                                                            >  1:94494/1‑301 (MQ=255)
                                                                                                                                                   cGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCCGTGTCGGCGGCACTTTGCGCCCGTTCCCAGGCCTGGCGGCCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTggcagaacttt                                                                                                                                                    >  1:167990/1‑293 (MQ=255)
                                                                                                                                                        ttCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTcc                                                                                                                                               <  2:412547/301‑1 (MQ=255)
                                                                                                                                                            ccgccccccacGACAAAACCACGGCCGGTGTCCGCCGGACGCGCGGCTTCTACCGCGGCATGCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGTTGATAATGTCGGCAATTTGCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAgcg                                                                                                                                           <  1:336324/294‑1 (MQ=255)
                                                                                                                                                               cccaccaCGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCgg                                                                                                                                        <  2:327011/300‑1 (MQ=255)
                                                                                                                                                                  atcgcggtaaagtcggggCCGTGTTCATCCGCTCTGTCGGCTTCAGCCTCGGCATTCAGCGCGTGGATATTAGGCGGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACGTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATg                                                                                                                                     <  1:141895/285‑1 (MQ=255)
                                                                                                                                                                                             cccGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGaa                                                                                                          >  2:338697/1‑300 (MQ=255)
                                                                                                                                                                                              cggCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGAGAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATc                                                                                                         <  1:55721/299‑1 (MQ=255)
                                                                                                                                                                                                         gttgaagccGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCCGTACGGGCGGCAAGAACGGGGTTGcccc                                                                                              >  1:634048/8‑300 (MQ=255)
                                                                                                                                                                                                           ttCAACCGCGGCATTCAGCCCGAGGCTATTAGTCTTGACTGCACTACCGTCGATAACGCTGATACTGTCGGCAATTTTCTTCGAACTATCGGCGCTCTCATGCATCGTTTTCACTCCGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGCGCTGCTCCATGCTGGCGGCAGTTTCTTCCAACGCGGATACCTGCTGCTCTGTACGGGAGCAAAGACAGGGGCTGCccgcc                                                                                            >  1:371822/1‑301 (MQ=255)
                                                                                                                                                                                                            tCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCTCCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCccgccg                                                                                           <  2:94618/301‑1 (MQ=255)
                                                                                                                                                                                                                cggggggaTTCAGGGGGAGGGTATTAGTTGGGAAGGCAATTCCGTCGATAAGGTTGGTAATGTCGGCAATTTTCTGCGAACGATGGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACGTGGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTGATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAt                                                                                       <  1:591258/294‑1 (MQ=255)
                                                                                                                                                                                                                ccGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCaaa                                                                                       >  2:359305/1‑300 (MQ=255)
                                                                                                                                                                                                                      cATTCAGCGCGAGGTTATTTGTTTGGAAGGCAATACCGTCGATACCGCTGTTAATGTCGGCAATTTTCTTCGAACTTTCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACg                                                                                 <  1:275848/301‑1 (MQ=255)
                                                                                                                                                                                                                      cATTCAGCGCGAGGATATTATTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACg                                                                                 <  2:483367/301‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTAGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCTGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGTCGGCAGTTTCTGCCAGCGCGGAGTCCTGCTGTTCCGTCAGGGCGGAAAGATCGGTGATACCCGCCTCAATTTCACGGGTACCggact                                                                      >  1:521208/1‑298 (MQ=255)
                                                                                                                                                                                                                                     ttgctcttccgatctcAATACTGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGTGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGAt                                                                  <  2:388104/286‑1 (MQ=255)
                                                                                                                                                                                                                                       ttAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCCGTACGGGAGGAACGATCGGTGTTGCCCGCCGCAATTTCACGGGTACAGGCATaactgg                                                                >  1:544965/1‑296 (MQ=255)
                                                                                                                                                                                                                                                    cAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGTGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGACATTTCACGGGTACAGGCATAGAtagaccggaagagcgg                                                   >  1:387590/1‑286 (MQ=255)
                                                                                                                                                                                                                                                     aaTACCGTCGATAACGCTGATAATGTCGGGAATTTTTTTCGAACTGTCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTcg                                                  <  1:125412/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        gtcgttgttAACGCTTATAATGGTGGCAATTTGCTTCGAACTGTCGGCGATCTCATGCATCGTTTGCACGACTCCAGCCACCACTTGGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTGGCGCCAGTGGCGAGGCCTGGCGGGCGTTAGCGGCGGTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCGTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCgg                                               <  1:84334/293‑1 (MQ=255)
                                                                                                                                                                                                                                                             cGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTCTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGCTTCTTCCAGCCCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCCCGGCTACCGGCATAGATGGCATCTGAACCTCCGCGGACAtt                                          >  1:384438/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                     ttcagacgtgtgctcttccgatctgAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACgg                                  <  2:140559/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttccctacctgacgctcgtccgcttttCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTGTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTgtc                              <  1:420962/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGGTTTCACTACGCCAGCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTGCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTgtc                              <  2:372318/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aaTTTTCTTCGAACTATCGGCGAGTTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAagag                    <  2:326518/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         tttcGAATTATCGGCGATCTCATGCATCGGTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGc              <  2:634925/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             gAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACggaaacaggaagagcgacgtgaaggg          >  1:140434/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                                                 tATCGGTGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGa      <  2:532626/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   tCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTgtcagatcggaagagcacacgtctgaaca    >  2:421525/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                                                     ggCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGaaaag  >  2:232556/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                          
GTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGCATCGTTTTCACTACGCCATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACG  >  NC_000913/1971340‑1971933

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