Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,603,735 G→T 5.3% G239G (GGG→GGT 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,7350GT 5.3% 273.4 / 5.1 80G239G (GGG→GGTlsrCautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base G (36/39);  new base T (2/3);  total (38/42)
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.02e-01

CGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTT  >  NC_000913/1603474‑1604008
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                  
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)
                                                                                                            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                                                                                                                                                                                                                                                         <  6:132908/149‑1 (MQ=255)
                                                                                                                                            tGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGc                                                                                                                                                                                                                                                         >  5:10268/1‑149 (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:187001/1‑149 (MQ=255)
                                                                                                                                                                aTTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGttt                                                                                                                                                                                                                                     <  8:464720/149‑1 (MQ=255)
                                                                                                                                                                 ttCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAg                                                                                                                                                                                                                                                                       <  4:517029/114‑1 (MQ=255)
                                                                                                                                                                     catgttgtcATTTTCTTTGAACGGCTTCTTGGCGGCACTGGCGGGAATTCGTTTTGCTTCGCAGATTGGTTTTCTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCCGCCTGCGTGCTGGGCGGCATTCGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                                                                                                  <  2:267141/267‑1 (MQ=255)
                                                                                                                                                                        tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtg                                                                                                                                                                                                                             <  5:145912/149‑1 (MQ=255)
                                                                                                                                                                        tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtg                                                                                                                                                                                                                             <  8:364604/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                                                                                                                                                                                                           >  6:87659/1‑142 (MQ=255)
                                                                                                                                                                                                 aTGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTg                                                                                                                                                                                                           <  5:87659/142‑1 (MQ=255)
                                                                                                                                                                                                  tGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTa                                                                                                                                                                                                   >  5:335137/1‑149 (MQ=255)
                                                                                                                                                                                                    gcggcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACt                                                                                                                                                                                                 <  8:366794/149‑1 (MQ=255)
                                                                                                                                                                                                       gcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                                                                                                                                                              >  7:18812/1‑117 (MQ=255)
                                                                                                                                                                                                       gcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggt                                                                                                                                                                                                                              <  8:18812/117‑1 (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                                                                                                                                                                     <  5:350319/126‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGAt                                                                                                                                                                     >  6:350319/1‑126 (MQ=255)
                                                                                                                                                                                                                                                          aCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCg                                                                                                                                                                       >  1:22443/1‑121 (MQ=255)
                                                                                                                                                                                                                                                              gTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAg                                                                                                                                                                   <  3:22443/121‑1 (MQ=255)
                                                                                                                                                                                                                                                               tACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTCGTGCGGTACTCGGCGCATGGTTCCTGACGCACCTCGATAGCGTACTGGTGCTGTTGCGCACTCCGGCATGGTGGAATGATTGTATAGCTGGTCTGTTTCTGCTTGCGGTTCTGGTGTTCGATGGATCCTGGCTTTGTGATCGGGAACGTTGTCTCTGCGGGAAAAACAATGCCCCCTTACTGAcacttctcgtc        >  4:62073/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                   gggCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGAAAAAATATGcc                             <  2:62074/249‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggCTGGAGATGAAAGAAATTGCAGCCTGCGTGCTGGGCGTCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGGGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGtt  <  1:4358/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                  
CGCCCCGCTGCTGCTTGGCGTTTCAGCAATTGGTTGGTTGACGATAATTCTGGTGGCATTTATGGCCTGGCTGCTGGCAAAGACGGCGTTTGGACGCAGTTTTTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGCGGCAAAAATATGCCCGCTTTATGACGCCACCGCCATCCGTT  >  NC_000913/1603474‑1604008

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