Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,603,724 A→C 5.6% T236P (ACC→CCC)  lsrC → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,7240AC 5.6% 267.0 / 7.1 80T236P (ACC→CCC) lsrCautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base A (34/40);  new base C (2/3);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

CGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCC  >  NC_000913/1603456‑1603981
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                 
cgAACTGAAACAGCTCTCCGCCCCGCTGCTGATTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGGCGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTa                                                                                                                                                                                                                                                             <  2:73578/275‑1 (MQ=255)
    cTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCGGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCg                                                                                                                                                                                                                                                          <  2:110689/274‑1 (MQ=255)
                  cGCCCCGTTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTGATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAg                                                                                                                                                                                                                                           <  2:362631/275‑1 (MQ=255)
                                      tttCAGCAATTGGTTGGTTGGCGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTg                                                                                                                                                                                                                       <  2:448999/275‑1 (MQ=255)
                                      tttCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGGCGCCTGCGTGCTg                                                                                                                                                                                                                       >  2:44559/1‑275 (MQ=255)
                                                                     ggtggCATTAAGTGCCGGGCTGCGGGCAAAGCTGGCGTTTGTACGCGGTTTTTATGCCACGGGCGATAATTTACGGGGCGCGCGGCAACTGGGCGGTCGTACTTAGGCCAGTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                         <  4:147342/274‑1 (MQ=255)
                                                                           aTTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGGTCCGGTTCGATcc                                                                                                                                                                                  >  2:301523/1‑274 (MQ=255)
                                                                                                            tGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGTTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                                                                                                                                  <  2:292428/274‑1 (MQ=255)
                                                                                                                            gCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCgg                                                                                                                                                                                                                                                               <  6:79579/149‑1 (MQ=255)
                                                                                                                              cccGGGGGATATTTTCCAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTgg                                                                                                                               <  2:498819/273‑1 (MQ=255)
                                                                                                                                      aTAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTg                                                                                                                                                                                                                                                     <  5:192625/149‑1 (MQ=255)
                                                                                                                                       tAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTgg                                                                                                                                                                                                                                                    <  8:7020/149‑1 (MQ=255)
                                                                                                                                       tAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTgg                                                                                                                                                                                                                                                    >  7:98228/1‑149 (MQ=255)
                                                                                                                                               aGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                                                                                                                                            <  8:138080/149‑1 (MQ=255)
                                                                                                                                                 ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTgg                                                                                                                                                                                                                                                    >  4:385752/1‑139 (MQ=255)
                                                                                                                                                 ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGc                                                                                                                                                                                                                                          >  5:191117/1‑149 (MQ=255)
                                                                                                                                                 ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCCTTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGCTCGATCGCGTACTGGTGCTGTTGCGCATTCCGGc                                                                                                            >  2:389041/1‑275 (MQ=255)
                                                                                                                                                   ggtgagTCAATTGGGTGTTCGTACGGAGGCCATTTGAGTATTGGAACTTCGGGGCAACGGCTGAGTGGGGGAACGGTTGGGAACGGGGTTTGCTTCGCAGATTGGTTGTATCCCCAACCTGCCCGGGACCAGGCTGGGAATGCAAGCAATTGCAGCCTGCGTGCTGGGCGCAATTAGTTTGCTCGGTGGTTCCGGTGCGATCATCCGTGTGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCa                                                                                                           <  4:370471/269‑1 (MQ=255)
                                                                                                                                                     cTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAt                                                                                                                                                                                                                                                <  2:385757/139‑1 (MQ=255)
                                                                                                                                                     cTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTc                                                                                                                                                                                             >  4:349175/1‑190 (MQ=255)
                                                                                                                                                        gTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                                                                                                                                   >  7:464720/1‑149 (MQ=255)
                                                                                                                                                        gTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                                                                                                                                   >  5:132908/1‑149 (MQ=255)
                                                                                                                                                         tCAACTGGGCGTTCGTACTTAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGCTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAg                                                                                                                                                                                                                                  >  6:85437/1‑149 (MQ=255)
                                                                                                                                                         tCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                         <  2:349180/190‑1 (MQ=255)
                                                                                                                                                              tGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGc                                                                                                                                                                                                                             >  5:10268/1‑149 (MQ=255)
                                                                                                                                                              tGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGc                                                                                                                                                                                                                             <  6:132908/149‑1 (MQ=255)
                                                                                                                                                                        tACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAt                                                                                                                                                                                                                                                <  6:226611/120‑1 (MQ=255)
                                                                                                                                                                        tACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAt                                                                                                                                                                                                                                                >  5:226611/1‑120 (MQ=255)
                                                                                                                                                                           tGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                                                                                                                                >  6:204774/1‑149 (MQ=255)
                                                                                                                                                                             aaGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCAtt                                                                                                                                                                                                              >  6:145912/1‑149 (MQ=255)
                                                                                                                                                                              aGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                                                                                                                                             <  7:119417/149‑1 (MQ=255)
                                                                                                                                                                              aGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                                                                                                                                             <  7:383106/149‑1 (MQ=255)
                                                                                                                                                                               gCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                                                                                                                               >  2:517034/1‑114 (MQ=255)
                                                                                                                                                                                ccATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGt                                                                                                                                                                                                           <  5:161178/149‑1 (MQ=255)
                                                                                                                                                                                 cATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                                                                                                                                          >  8:249200/1‑149 (MQ=255)
                                                                                                                                                                                  aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGttt                                                                                                                                                                                                         <  8:464720/149‑1 (MQ=255)
                                                                                                                                                                                  aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGttt                                                                                                                                                                                                         >  8:187001/1‑149 (MQ=255)
                                                                                                                                                                                   ttCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAg                                                                                                                                                                                                                                           <  4:517029/114‑1 (MQ=255)
                                                                                                                                                                                       catgttgtcATTTTCTTTGAACGGCTTCTTGGCGGCACTGGCGGGAATTCGTTTTGCTTCGCAGATTGGTTTTCTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCCGCCTGCGTGCTGGGCGGCATTCGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                                                                      <  2:267141/267‑1 (MQ=255)
                                                                                                                                                                                          tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtg                                                                                                                                                                                                 <  8:364604/149‑1 (MQ=255)
                                                                                                                                                                                          tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtg                                                                                                                                                                                                 <  5:145912/149‑1 (MQ=255)
                                                                                                                                                                                            tgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTTCTCggtggt                                                                                                                                                                                               >  6:429943/1‑149 (MQ=255)
                                                                                                                                                                                              gCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                                                                                                                                          >  8:400963/1‑136 (MQ=255)
                                                                                                                                                                                              gCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                                                                                                                                          <  7:400963/136‑1 (MQ=255)
                                                                                                                                                                                                  tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGTGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGt                                                                                                                                                                                         <  6:10268/149‑1 (MQ=255)
                                                                                                                                                                                                      gTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGa                                                                                                                                                                                     >  5:296820/1‑149 (MQ=255)
                                                                                                                                                                                                        tGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATc                                                                                                                                                                                   <  6:425119/149‑1 (MQ=255)
                                                                                                                                                                                                          aaCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCAt                                                                                                                                                                                 >  8:128995/1‑149 (MQ=255)
                                                                                                                                                                                                           aCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggtt                                                                                                                                                                                                                                                  <  5:330694/83‑3 (MQ=39)
                                                                                                                                                                                                           aCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggtt                                                                                                                                                                                                                                                  >  6:330694/1‑81 (MQ=39)
                                                                                                                                                                                                                  cATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGGACTGGTGCTGTTGCGCATTCCGGCATGGTGGAAGGATTTTATCGCGGGTGTGGGTTTGCTGGCGGTGGGGGTGTTTGATGGACGCGTgt                                           >  1:93311/1‑274 (MQ=255)
                                                                                                                                                                                                                   aTGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTg                                                                                                                                                                               <  5:87659/142‑1 (MQ=255)
                                                                                                                                                                                                                   aTGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTg                                                                                                                                                                               >  6:87659/1‑142 (MQ=255)
                                                                                                                                                                                                                    tGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTa                                                                                                                                                                       >  5:335137/1‑149 (MQ=255)
                                                                                                                                                                                                                      gcggcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACt                                                                                                                                                                     <  8:366794/149‑1 (MQ=255)
                                                                                                                                                                                                                         gcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                                                                                                                                  <  8:18812/117‑1 (MQ=255)
                                                                                                                                                                                                                         gcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                                                                                                                                  >  7:18812/1‑117 (MQ=255)
                                                                                                                                                                                                                                gggtaattggtgttGTTTCGCCGTGTGGTTTTTGCCCCAACCAGACCGGTACTGGGCTGGGGATGAAGTCAATGTGAGCCTGCGTGCTGGGTGGCATTAGTTGGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGa                              <  4:301519/262‑1 (MQ=255)
                                                                                                                                                                                                                                 gggAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGt                                                                                                                                                          <  7:249200/149‑1 (MQ=255)
                                                                                                                                                                                                                                           tttGCTTCGCAGATTGGCTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCa                                                                                                                                                >  8:295153/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTAcggcg               >  2:196199/1‑275 (MQ=255)
                                                                                                                                                                                                                                               cTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATc                                                                                                                                            <  7:128995/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTc                                                                                                                                                        >  3:47896/1‑133 (MQ=255)
                                                                                                                                                                                                                                                   gCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATa                                                                                                                                        <  5:204774/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGc                                                                                                                                      >  5:18220/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aTTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGa                                                                                                                                                    <  1:47896/133‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                                                                      <  2:498554/211‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttGGTTGTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTa                                                                                                                                   <  6:18220/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACt                                                                                                                                 <  5:85437/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTTTTCGCCCCATCCGGGCCGGTTCGGGTNTGTGGTTGAAAGCACTGGCGGCCAGNTGTCTGGGCGGCATTGGTTTGTTCGGGGGTTCCGGGGTGTTCATTGGTGCGTTGCTCGGCGCATGGTTCCTATCGCAGATNGATAGCGTACTGGGGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGcc  <  3:93311/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCg                                                                                                                                     >  7:280626/1‑144 (MQ=255)
                                                                                                                                                                                                                                                           gTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCg                                                                                                                                     <  8:280626/144‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTg                                                                                                                                >  8:115706/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGt                                                                                                                         >  2:515712/1‑150 (MQ=255)
                                                                                                                                                                                                                                                                  ccccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGt                                                                                                                         <  5:429943/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg                                                                                                                       >  8:45402/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    atctccAGCCCGGTACCGGTCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg                                                                                                                       >  5:4004/5‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgcg                                                                                                                     <  4:515707/150‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGAt                                                                                                                                         >  6:350319/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                         cAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGAt                                                                                                                                         <  5:350319/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCg                                                                                                                                           >  1:22443/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                 
CGAACTGAAACAGCTCTCCGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCC  >  NC_000913/1603456‑1603981

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: