Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 876,917 C→T E373E (GAG→GAA rimO ← ribosomal protein S12 methylthiotransferase RimO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913876,9170CT98.3% 191.6 / ‑3.9 58E373E (GAG→GAArimOribosomal protein S12 methylthiotransferase RimO
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/27);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCA  >  NC_000913/876624‑877212
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                       
gagTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTcggcg                                                                                                                                                                                                                                                                                                   >  1:68212/1‑300 (MQ=255)
         aaaCGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTAGTCCACTTCGTAGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATctgc                                                                                                                                                                                                                                                                                          >  2:352132/1‑300 (MQ=255)
              gCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATctgctgctg                                                                                                                                                                                                                                                                                     <  1:675383/300‑1 (MQ=255)
                        gATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCAGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATg                                                                                                                                                                                                                                                                           >  2:336326/1‑300 (MQ=255)
                            aGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAAc                                                                                                                                                                                                                                                                       <  1:107008/300‑1 (MQ=255)
                            aGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAAc                                                                                                                                                                                                                                                                       <  2:408043/300‑1 (MQ=255)
                                ggCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGACTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGtt                                                                                                                                                                                                                                                                   >  1:186755/1‑300 (MQ=255)
                                    gTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAg                                                                                                                                                                                                                                                               >  1:620826/1‑300 (MQ=255)
                                        tCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCCGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGtt                                                                                                                                                                                                                                                           >  2:67397/1‑300 (MQ=255)
                                         cATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGttc                                                                                                                                                                                                                                                          <  2:438684/300‑1 (MQ=255)
                                          attattCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCAACATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCACACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGttct                                                                                                                                                                                                                                                         >  2:9567/1‑300 (MQ=255)
                                                   aGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAAct                                                                                                                                                                                                                                                <  1:598454/300‑1 (MQ=255)
                                                    ggCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAActt                                                                                                                                                                                                                                               >  1:606909/1‑300 (MQ=255)
                                                    ggCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTATTTAActt                                                                                                                                                                                                                                               >  1:691577/1‑300 (MQ=255)
                                                    ggCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCAACATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAActt                                                                                                                                                                                                                                               >  1:671936/1‑300 (MQ=255)
                                                              aTGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAAc                                                                                                                                                                                                                                     <  1:638333/300‑1 (MQ=255)
                                                                    ttGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGTCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTc                                                                                                                                                                                                                               <  2:380899/300‑1 (MQ=255)
                                                                     tGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCNTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTca                                                                                                                                                                                                                              >  2:152441/1‑300 (MQ=255)
                                                                             cTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGc                                                                                                                                                                                                                      <  2:314773/300‑1 (MQ=255)
                                                                              tGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCa                                                                                                                                                                                                                     <  2:393512/300‑1 (MQ=255)
                                                                               gAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCAt                                                                                                                                                                                                                    >  1:184705/1‑300 (MQ=255)
                                                                                  ttaattaaACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCGCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTgg                                                                                                                                                                                                                 <  1:99252/300‑1 (MQ=255)
                                                                                    aattaaACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCg                                                                                                                                                                                                               <  2:400954/300‑1 (MQ=255)
                                                                                    aattaaACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCg                                                                                                                                                                                                               <  1:182404/300‑1 (MQ=255)
                                                                                           cccGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCAc                                                                                                                                                                                                        <  2:376279/300‑1 (MQ=255)
                                                                                                     ccACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAAAGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTATTCGTCAACTTCGGCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCCTTGGCGTCTGCACCTTCAACCgg                                                                                                                                                                                              >  2:332439/1‑300 (MQ=255)
                                                                                                         aaaTCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTg                                                                                                                                                                                          >  1:192624/1‑300 (MQ=255)
                                                                                                          aaTCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGt                                                                                                                                                                                         >  1:151819/1‑300 (MQ=255)
                                                                                                                 cTCATCGGTGTGCTCGACTTTCATACGCGGGATATCACCCGGCTTAGCGTTGGTTTCACCATTTAGATAAACCGCGCCGTCGATTTCCGGTACATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAACCACCAGAACTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTCCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTaaa                                                                                                                                                                                  <  2:104238/300‑1 (MQ=255)
                                                                                                                            gCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAAcgc                                                                                                                                                                       <  1:59247/300‑1 (MQ=255)
                                                                                                                                  cTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCa                                                                                                                                                                 >  1:536953/1‑300 (MQ=255)
                                                                                                                                      cacaCGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGAcg                                                                                                                                                             <  1:422158/300‑1 (MQ=255)
                                                                                                                                      cacaCGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGAcg                                                                                                                                                             <  2:337495/300‑1 (MQ=255)
                                                                                                                                              ggATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGTTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGAcgcatttctt                                                                                                                                                     >  2:284465/1‑293 (MQ=255)
                                                                                                                                              ggATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTNCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGcttctt                                                                                                                                                     >  1:502217/1‑300 (MQ=255)
                                                                                                                                              ggATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGcttctt                                                                                                                                                     >  1:275161/1‑300 (MQ=255)
                                                                                                                                              ggATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGcttctt                                                                                                                                                     >  2:19019/1‑300 (MQ=255)
                                                                                                                                              ggATATCAACCGGCTTAACGTTGGTTTCACCGTTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTT‑‑AGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCNCCTTCAAACGGGCTGTATTTAAAGCAGCAAACGCGATCCAGACGCGCNTCTTTc                                                                                                                                                   >  2:576649/1‑300 (MQ=255)
                                                                                                                                                          gCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGa                                                                                                                                         >  2:127444/1‑300 (MQ=255)
                                                                                                                                                          gCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCNTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGa                                                                                                                                         >  1:472413/1‑300 (MQ=255)
                                                                                                                                                              aaCGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAg                                                                                                                                     <  1:264335/300‑1 (MQ=255)
                                                                                                                                                                                              cGTCAATTTCCGGTGTATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCTCACTTTCTCTTGTAGGCGTTCTGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGa                                                                                                     <  2:98451/300‑1 (MQ=255)
                                                                                                                                                                                               gTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGaa                                                                                                    >  1:125420/1‑300 (MQ=255)
                                                                                                                                                                                                          gTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATaa                                                                                         <  1:599149/300‑1 (MQ=255)
                                                                                                                                                                                                             cATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAgg                                                                                      >  1:161266/1‑300 (MQ=255)
                                                                                                                                                                                                                  gCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAg                                                                                 >  1:288460/1‑300 (MQ=255)
                                                                                                                                                                                                                  gCCATGCTGCGACCAACCGCGCCTTCTTCGTCCACTTCGTAGATAATCACCCGAATTCCACGGCCCACTTTCTATTGCAGGCGTTCGGCGGAAATCTACTGCTGCAGCTGCATGAAACGGTTCCNGCGTTCTTCTTTAACTTCTTACGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCANCCGGGCTGTATTTAACGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCAGACAATAAAGGTTGAg                                                                                 >  2:545847/1‑300 (MQ=255)
                                                                                                                                                                                                                   ccATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGc                                                                                <  1:513907/300‑1 (MQ=255)
                                                                                                                                                                                                                          gcGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGt                                                                         <  1:413889/300‑1 (MQ=255)
                                                                                                                                                                                                                                aaTCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTc                                                                   <  1:513008/300‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTACGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTcc                                                                  >  2:85083/1‑300 (MQ=255)
                                                                                                                                                                                                                                                            aTAATCCCCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATNCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCANTTCCGGGCAGATTTCGCGCCACTGTTTGATg                                       <  2:492024/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gAATTTCACGGCCCACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCNAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGtt                             >  1:467276/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                gCCCACTCTCTATTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTAAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGGTTTCGCGCCACTGTTTGATGCGTGCCAGTTGGCGATCTAc                   >  2:343189/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                   cACTTTCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGa                <  2:6193/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     cTTTCTNTTGCATGCGTTCGTCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTNCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAAc              <  2:564959/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ttCTCTTGCAGGCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAAccc            <  1:71843/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 ggCGTTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTtca  >  2:618328/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                       
GAGTAAGAAAAACGGCAGGGCACAGATTAGCTGGCAGTAATCATTATTCGCAGGCCGGACAAATGATTTTGCCCGGCCTGAATTAATTAAACCCGGCTACCCCACAAATCGTACTCATCGGCGTGCTCGACTTTCACACGCAGGATATCACCCGGCTTAACGTTGGTTTCACCATTGAGATAAACCGCGCCGTCGATTTCCGGTGCATCTGCCATGCTGCGACCAATCGCGCCTTCTTCGTCCACTTCGTCGATAATCACCAGAATTTCACGGCCCACTTTCTCTTGCAGGCGCTCGGCGGAAATCTGCTGCTGCAGCTGCATGAAACGGTTCCAGCGTTCTTCTTTAACTTCTTCCGGAACCTGGTCAGGCAGGGCATTGGCGTCTGCACCTTCAACCGGGCTGTATTTAAAGCAGCCAACGCGATCCAGACGCGCTTCTTTCAGGAAGTCGAGTAGCATCTGGAAATCTTCTTCTGTCTCGCCAGGGAAGCCGACAATAAAGGTTGAGCGTAGGGTCAGTTCCGGGCAGATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCA  >  NC_000913/876624‑877212

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