Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,521,502 A→T 8.8% intergenic (+174/+91) lysV → / ← xapR tRNA‑Lys/transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,5020AT 8.8% 203.9 / 11.5 69intergenic (+174/+91)lysV/xapRtRNA‑Lys/transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base T (3/3);  total (33/36)
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.55e-03

CCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAGTTCATTTGTGCGTAGCTGTCCGCTATCAATGTGATACCA  >  NC_000913/2521233‑2521752
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
ccGTCATCACCCCCCCCTTCGGGTCGTTAGCTCAGTTGGTAGTGCAGTTGCCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGAACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGTGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                                                                                                                        <  2:68177/274‑1 (MQ=255)
    catcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCt                                                                                                                                                                                                                                                    <  2:460377/274‑1 (MQ=38)
                                        agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTNAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGGGAGTCATCCNGCACGACCCACAACTGTAACAAAACACCCTAACGCCCCCTTTTTTNTNTCTGCGCTTTTCCAAATTTCCTgtacc                                                                                                                                                                                                                >  3:42958/1‑270 (MQ=11)
                                           gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAATGATTCGAACCCGCAGCAGGTTTGAGTCGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCCAGCGAAGCGAGTCACCCTGCACGACCCACCAATCTAAACAAGCGCCCTAAAGCCGCTTTTTTGCTATCTGCCATTTGCGTAATTGCCTGATGCGCt                                                                                                                                                                                                             >  2:71718/1‑274 (MQ=14)
                                                  aCTTTTAATCAATTGTTCGCAGTTTCCAATCCTGCACTCCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTTAGTTTAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg                                                                                                                                                                                                     <  2:311876/275‑1 (MQ=255)
                                                      ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                                                                                                                                                                                                        >  2:359251/1‑252 (MQ=39)
                                                      ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGAACGACCCCCCAATGTAACAAACCGCCCTCAAGGCGCCTTTTTGCTATCTCGCATTTTCGAAATTGCCTGATGCGCTTACCTTTGCAgc                                                                                                                                                                                                 >  4:183362/1‑274 (MQ=14)
                                                              ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                                                                                                                          >  2:195223/1‑210 (MQ=31)
                                                                           cGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACACACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTc                                                                                                                                                                              >  2:202635/1‑273 (MQ=255)
                                                                              aTCCTGAACAACCCAACAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACGCAAAGATTCGAGCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGtc                                                                                                                                                                         <  2:335312/275‑1 (MQ=255)
                                                                                 cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGtctc                                                                                                                                                                       <  2:313459/274‑1 (MQ=255)
                                                                                  tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcct                                                                                                                                                                     >  2:451129/1‑275 (MQ=255)
                                                                                  tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGCAGCCGCTCGCGGCGACGGCCCGAAGGGCGCGCGAGTCGAGTCATCCGCCACGCACCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTTGCTTCGATTTGCGCAATTGCATGATTCGCTGCACTTTGCAGACTACTCTTTCCCGTACTTCCCGtctctc                                                                                                                                                                     >  4:485133/1‑273 (MQ=14)
                                                                                  tGCACGACCCACAAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcct                                                                                                                                                                     <  2:261836/275‑1 (MQ=255)
                                                                                              aaaTGTAACAAAGCGCCCTACAGGCGCTTTTTTACTATCTGCGATCATCAAAGATTCGACCCTGCAGCAGGTGTGAGTTGAGCACAGCGAAACAACGGAGCCGCTCGCCGCGACGGCCCGAAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTgtgt                                                                                                                                                         <  1:54079/274‑1 (MQ=18)
                                                                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGCCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCCATg                                                                                                                                                    >  2:200424/1‑275 (MQ=255)
                                                                                                         aGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCa                                                                                                                                              <  2:376128/275‑1 (MQ=255)
                                                                                                             cccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTGCAACCTGCAGCAGCTTTTAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCCGCTACGGCCCGACGGGCGAGCGGAGCGAGTCATCCTGACCCATCCGACAATGTAAAAAAGCGCCCTAAAGGAGCTTTTTTGCTATCTGGGATTTGCGAAATTGCATGATTCGCTTCACTTACAAGATTAATCTTTCCGGCCATTCCTGTGTCCTCACCTACTTCTTCAATCCAGCCcaca                                                                                                                                           >  4:65566/1‑274 (MQ=17)
                                                                                                                        cTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCTAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCCCTTTTTTGCTATTTGGGCTTGTCACAATTGCCTGATCCCCTTTCTTTATCCGCCTACTATTTCCGGCAATTTCCGGCCCCTACCCTCTGGTGTCAATGCACCCAACCACTTAACCATc                                                                                                                                >  4:208474/1‑274 (MQ=14)
                                                                                                                        cTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCCGCGATTTGCGAAATTGCCTTATCCGCTTCACTTAACAGACTACTATTTCCGGCAATTCTTGTCTCCCCACCTACTGTGTCAATGCAGCCAACAGCTCAACCaac                                                                                                                                >  4:54068/1‑272 (MQ=255)
                                                                                                                          tttttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGGGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATcgcg                                                                                                                             <  2:54069/275‑1 (MQ=255)
                                                                                                                                cTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTc                                                                                                                       <  2:183364/275‑1 (MQ=255)
                                                                                                                                   tCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTACCCATCGCTGGCCTCAcc                                                                                                                    >  2:393921/1‑275 (MQ=255)
                                                                                                                                   tCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCACTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTCGCAGCCTACTAATTCCGGCAATTCCTGTCTCCTCACCTACCGTTTCAATGCAGCCAACAGCTTAACCATCCCCGGCGTCaac                                                                                                                    >  4:316903/1‑273 (MQ=255)
                                                                                                                                    cTGGGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCAcc                                                                                                                    <  2:485138/274‑1 (MQ=255)
                                                                                                                                                 aGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                                                                                                                                                                                                                   >  4:304385/1‑166 (MQ=255)
                                                                                                                                                   gttCGAACCTGCAGCAGGTTTGAGTTGAGCGCGGCGAAACAACGAAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAGAGGCGCTTTTTTACTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                                                    <  2:340381/274‑1 (MQ=255)
                                                                                                                                                   aTTCGAGCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCTAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                                                    <  2:65567/275‑1 (MQ=255)
                                                                                                                                                     tCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGGGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCACTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAAc                                                                                                   <  2:316907/274‑1 (MQ=255)
                                                                                                                                                     tCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                                                                                                                                                                                                               <  2:304389/166‑1 (MQ=255)
                                                                                                                                                      cGAACCTGCAGAAGGTTTGAGTGAGGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACa                                                                                                  <  2:6553/274‑1 (MQ=255)
                                                                                                                                                             gcggagggggtgagatgagCGCAGCAAAAAAACGGAGGCGCGAGAGGCGACGACCGAAAGGGCAGGCGAAGCGCGTCAGCCTGCCCGACACACCAATGTAAAAAAGCGCCGAAAAGGCGCATTTTTGCGATCTGGGATTTGCGAAAGTGCCTGGTGCGCTTCACTGTACAGACTGATAATACCGACAAGTCCTGTCGACTCACCTGCTGTGTCAAGTCAGCCAACAGCTTAACCATCGCGGGCGTCTCCTGCTGTGTTTCATAAACAATATaaaa                                                                                          <  4:129138/265‑1 (MQ=255)
                                                                                                                                                             gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGCAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCTCCTGCTGTGTTTCATATACAatat                                                                                              >  2:129139/1‑271 (MQ=255)
                                                                                                                                                             gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATaa                                                                                          <  1:42958/275‑1 (MQ=255)
                                                                                                                                                             gcagaaGGTTGGAGTGAAGCGAAGCGAAACACCGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATaa                                                                                          <  2:208477/275‑1 (MQ=255)
                                                                                                                                                              cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                                                                                                                                                                                                       >  6:439318/1‑149 (MQ=255)
                                                                                                                                                                             gggggaaGCGAAACAATGGAGCTGCTCGGGGCGACGGCCAGAGGGGCGGGCGAAGCGAGTCATCCGGACCTACCCGCCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCTCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATgcg                                                                           <  4:200422/268‑1 (MQ=255)
                                                                                                                                                                              agCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTCTCa                                                                                                        >  2:379240/1‑244 (MQ=255)
                                                                                                                                                                                  cAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCa                                                                                                                                              <  2:263656/202‑1 (MQ=255)
                                                                                                                                                                                      gAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGact                                                                                                                                                                                               >  8:95371/1‑149 (MQ=255)
                                                                                                                                                                                                  cgcgccgcgACGACGGCCCGAAGGGCGAGCAACGTGATGCATCCTCCAGGACCCACCATTGTAAAAAAGCGCCCTAAAGGCGCTTTTTGGCTATCTGCGATTTGCGAAATGGCTTGATGAGCTTCACTTAGCAGACTCCTATTTCCGGCAATTCCTGTCCCCTCACCTACTGGGTAAATGCAGCCAACAGCTTAACCTTCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAAt                                                      <  4:451124/269‑1 (MQ=255)
                                                                                                                                                                                                     cTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCa                                                                                                                                                                   >  1:149582/1‑162 (MQ=255)
                                                                                                                                                                                                         cggcgACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCt                                                                                                                                                                            <  7:95371/149‑1 (MQ=255)
                                                                                                                                                                                                         cggcgACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTa                                                                                                                                                               <  3:149582/162‑1 (MQ=255)
                                                                                                                                                                                                               cGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAg                                         <  2:185330/274‑1 (MQ=255)
                                                                                                                                                                                                                            gagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAAt                                                                                                                                                                                >  1:66209/1‑126 (MQ=255)
                                                                                                                                                                                                                             agcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAACTTGCCTGATGCGCTACACATAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGGCCCCTGCTGTGTTTCATAAACAATATATAACCTTGCAGGGATCCGCTGTTCGAGCGGTCGGAAAATGACACCCGGCACTTTCATTTGTGCGTAg                          >  4:408221/1‑275 (MQ=38)
                                                                                                                                                                                                                                gaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCt                                                                                                                                                                            <  3:66209/126‑1 (MQ=255)
                                                                                                                                                                                                                                         cATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGCCGCTTTTTTGCGATCTGCGATTTGCGAAATTGCCGGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAGTTCATTTGTGCGTAGCTGTCCGCTAt               <  2:283865/274‑1 (MQ=255)
                                                                                                                                                                                                                                           tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAg                                                                                                                                          >  7:138557/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCAc                                                                                                                                                                  >  3:24463/1‑115 (MQ=255)
                                                                                                                                                                                                                                                     cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCGAGCCAGTTCATTTTGGCGTAGCTGTCCGCTATCAATGTGATACCa  >  2:58115/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATcgc                                                                                                                              >  6:108890/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTc                                                                                                                                                                                       >  2:516534/1‑91 (MQ=255)
                                                                                                                                                                                                                                                         ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTAc                                                                                                                                                              <  1:24463/115‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGc                                                                                                                                                                                   <  4:516529/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcct                                                                                                                                                                     >  5:209368/1‑101 (MQ=39)
                                                                                                                                                                                                                                                                aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcct                                                                                                                                                                     <  6:209368/101‑1 (MQ=39)
                                                                                                                                                                                                                                                                 aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCTGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTaa                                                                                                                                     >  5:216449/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCTGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTaa                                                                                                                                     <  6:216449/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCAcc                                                                                                                    >  7:193759/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACct                                                                                                                   >  5:71261/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACct                                                                                                                   <  5:439318/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaGCGCCCTAAAGGCGCTTTTTCGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTc                                                                                                                                                                              >  4:508753/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                     aGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACctgct                                                                                                                <  6:71261/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gcgcCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGt                                                                                                                                                                          <  2:508758/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cgcCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                                                                                                                                                                                                         >  1:142660/1‑58 (MQ=255)
                                                                                                                                                                                                                                                                         cccTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCggg                                                                                                                           >  4:231943/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                           cTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg                                                                                                                                                                                                     <  3:142660/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTc                                                                                                         <  8:138557/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aaaGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTc                                                                                                                       <  2:231947/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
CCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAGTTCATTTGTGCGTAGCTGTCCGCTATCAATGTGATACCA  >  NC_000913/2521233‑2521752

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