Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,398 T→G A228A (GCT→GCG nhaR → DNA‑binding transcriptional activator NhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091319,3980TG98.5% 218.3 / NA 68A228A (GCT→GCGnhaRDNA‑binding transcriptional activator NhaR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (23/44);  total (23/45)
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.79e-01

ATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTA  >  NC_000913/19153‑19641
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
atGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATTATGCCGCGTtgat                                                                                                                                                                                                                                                 <  1:535492/251‑1 (MQ=255)
atGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTtgat                                                                                                                                                                                                                                                 <  2:184256/251‑1 (MQ=255)
    tCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGaa                                                                                                                                                                                                                                              >  1:224140/1‑250 (MQ=255)
        ttCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCTAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCt                                                                                                                                                                                                                                          <  2:333971/250‑1 (MQ=255)
          cTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAAGGATTAAACGTAGAAATCCTCGGCGAGTTTGATTATGCCCCGTTGATATAAGCtttt                                                                                                                                                                                                                                       >  2:262493/1‑251 (MQ=255)
               tGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTACCGGCTTGTCTGGAAGAACGGCGACGTTTGATTACTGGGCGACGTTCAATGTTAGGGCGCCAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGGAGCCGCGTTGATGAAAGCTTTTGGttg                                                                                                                                                                                                                                  >  2:514340/1‑249 (MQ=255)
                   cgagaGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTTGTGTACAAGTCCACCACCAGAAAATCCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCGCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGAGATg                                                                                                                                                                                                                              <  2:212593/247‑1 (MQ=255)
                                     agGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACAGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGtt                                                                                                                                                                                                            >  1:332191/1‑251 (MQ=255)
                                     agGAAAGCCTGTTCTCCGTGAGAATTGGCGAATGTTGCGTGAGTGTCTGGTGTACAAATCCACCACCAGAAAAACCGTTTCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGGTTAACTCCCAGGGATTAAACGTAGAAAGCCTCGGCGAGTTTGAGGATGCCTCGTTGATGAAAGCGTTTGGTGCGATGCACAATGCAATCTTCGtt                                                                                                                                                                                                            <  2:53568/251‑1 (MQ=255)
                                           gCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCa                                                                                                                                                                                                      <  1:94412/251‑1 (MQ=255)
                                            ccTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCCATGTTAGGGCTCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCAGCGAGATTGATAATGACCCGTTGATGAAAGCTTTTGGTGCGTTGCACAATTCAATCTTCGGTGCCCCa                                                                                                                                                                                                      >  2:313335/1‑250 (MQ=255)
                                               gTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACg                                                                                                                                                                                                   <  1:330213/250‑1 (MQ=255)
                                               gTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACg                                                                                                                                                                                                   <  1:445981/250‑1 (MQ=255)
                                                                gCGAATGTGGCGTTATGTTCTGGTGTCCAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAATGTCGACTTTTGATTCCTTGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTTACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCTCAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACttt                                                                                                                                                                                 <  1:43072/251‑1 (MQ=255)
                                                                    aggttgCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATg                                                                                                                                                                             <  1:656115/246‑1 (MQ=255)
                                                                          cGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACCATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCggt                                                                                                                                                                        >  2:119309/1‑248 (MQ=255)
                                                                           gTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATa                                                                                                                                                                       >  1:527244/1‑250 (MQ=255)
                                                                            tGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCTAGGGATTAAACGTAGAAATCCTCGGCGTGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTAGTTGCCCCAAGGCTTTATGCATATGACTTTTATGACGATaa                                                                                                                                                                      >  1:617396/1‑250 (MQ=255)
                                                                            tGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATaa                                                                                                                                                                      >  1:617413/1‑250 (MQ=255)
                                                                             gAGTTGCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATaaa                                                                                                                                                                     <  1:22936/250‑1 (MQ=255)
                                                                              aGTTTCTGGTGTACATATACACCTACAAAAAAATCGTTCCCTTCTGTTTGGGAAGAACGGCGACTTTTGATTCCTGGGCGTCGTTCAATGTTAGGGCGCAAATTTCTTAACTGGTTTAACTCCCAGGGAGTAAACGTAGAAAGCCTCTGCGAGTCTGAGGATGCCGCGTTGATTAAATCTTTTTGTGTGTTGCATAAGGCAATCTTCGTTGGCCCAACGCTTTATACATATGACTTTTATGACGATAAAAc                                                                                                                                                                   <  2:670512/251‑1 (MQ=255)
                                                                                 ttCTGGTGTACACATCCACCACCAGAAAAAACGTTCACGGCTTGTCTGGAAGAACGGCCACTTTTGCTTCCTGGGAGACGGTCAATGCTAGGGCGCAAATCGCTTAACTGGTTTAACTCCCAGGGACTAAACGTAGAAATCCTCGGAGAGTTTGATGATGCCCCGTTGGTGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGccagttctgaagg                                                                                                                                                                >  2:56277/1‑240 (MQ=255)
                                                                                 ttCTGGTGTACAAATCCTCCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACGCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGcc                                                                                                                                                                           <  1:56277/240‑1 (MQ=255)
                                                                                 ttCTGGTGCACAAATCCACCACCAGAGAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGTCTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTg                                                                                                                                                                 <  2:451673/250‑1 (MQ=255)
                                                                                     ggTGTACAAATCCTCCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTGCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTc                                                                                                                                                               <  1:232593/248‑1 (MQ=255)
                                                                                                          aaaaaCCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGTGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAAt                                                                                                                                       <  1:437740/251‑1 (MQ=255)
                                                                                                            aaaCCGTTCCCGGCTTGTCGGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAAtgt                                                                                                                                     <  2:258984/251‑1 (MQ=255)
                                                                                                                cGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATg                                                                                                                                 <  1:469441/251‑1 (MQ=255)
                                                                                                                    cccGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGaa                                                                                                                              <  1:158967/250‑1 (MQ=255)
                                                                                                                        gCTTGGCTGGAAGAACGGCGACTTTTGATTCCTGTGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGGTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACGGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTa                                                                                                                         <  1:675133/251‑1 (MQ=255)
                                                                                                                         cTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTAc                                                                                                                        >  1:612947/1‑251 (MQ=255)
                                                                                                                               gggAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACGATGCt                                                                                                                  <  1:228011/250‑1 (MQ=255)
                                                                                                                                 gaagaaCGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAATCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCGCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTAt                                                                                                                <  2:55690/251‑1 (MQ=255)
                                                                                                                                   aagaaCGTCTACTTTTTATTCCTGGGATACGTGATATGTTAGGGCGCAAATTGCGTAACTGGTTTAACTCCCAGTGATTAAACGTAGAAATCCGCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTGCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCGAtt                                                                                                               <  2:474093/251‑1 (MQ=255)
                                                                                                                                   aagaaCGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTAtt                                                                                                               >  1:13536/1‑251 (MQ=255)
                                                                                                                                    agaaCGGCGACTGTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTAttt                                                                                                              <  2:641019/251‑1 (MQ=255)
                                                                                                                                           cGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGa                                                                                                       <  1:8291/251‑1 (MQ=255)
                                                                                                                                           cGACTTTTGATTCCTGGGCGACGGTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTACCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGAAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGCTGGAAGAGTACCATGCTATTTTTGCTgg                                                                                                       >  2:607267/1‑250 (MQ=255)
                                                                                                                                             ccTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGc                                                                                                     <  1:485571/250‑1 (MQ=255)
                                                                                                                                                 ttGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAAc                                                                                                                                                                   <  2:370903/185‑1 (MQ=255)
                                                                                                                                                 ttGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAAc                                                                                                                                                                   >  1:370903/1‑185 (MQ=255)
                                                                                                                                                  tGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTGTGCTGAGCGgat                                                                                                 <  2:445248/250‑1 (MQ=255)
                                                                                                                                                    atgccTGGGCGACGTTCAATGTTAGGGCTCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGgatgat                                                                                              <  2:671879/248‑1 (MQ=255)
                                                                                                                                                    aTTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCGTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGgatgat                                                                                              <  1:474530/251‑1 (MQ=255)
                                                                                                                                                       ccTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGAGATTAAACGGAGAAATCCTCGGCGATGTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTTAGCGGATGATTc                                                                                            <  2:639336/250‑1 (MQ=255)
                                                                                                                                                              gACGTTCAATGTTAGGGCGCAACTTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAACTTGGTCGCGTCgag                                                                                                                                          >  2:539815/1‑197 (MQ=255)
                                                                                                                                                              gACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCgag                                                                                                                                          <  1:539815/197‑1 (MQ=255)
                                                                                                                                                                cGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGCTTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTc                                                                                                                                                               >  2:255165/1‑174 (MQ=255)
                                                                                                                                                                cGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTc                                                                                                                                                               <  1:255165/174‑1 (MQ=255)
                                                                                                                                                                  ttCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGAt                                                                                                                                                                        >  1:18052/1‑163 (MQ=255)
                                                                                                                                                                  ttCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGAt                                                                                                                                                                        <  2:18052/163‑1 (MQ=255)
                                                                                                                                                                        gTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGc                                                                          <  1:121436/251‑1 (MQ=255)
                                                                                                                                                                                   aaTTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTAt                                                                                                                                                                                             <  1:614456/125‑1 (MQ=255)
                                                                                                                                                                                   aaTTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATACTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTAt                                                                                                                                                                                             >  2:614456/1‑125 (MQ=255)
                                                                                                                                                                                        cTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACACTGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGGAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCGGCAATCCggg                                                          >  2:582941/1‑250 (MQ=255)
                                                                                                                                                                                           aaCTGGTTAACCACCCAGGGATTAACCGTACAAATCCCCCGCGAGTTTGATGATGCCGCGCTGATGAAAGCTTTTGGTGCGATGCACTATGCAATCTCCGTTGCTCCAACGCTTTATGCATATGACTTTTATGCCGAGAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGCTGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAACGAATCTGCAATACGGatta                                                       <  2:1922/251‑1 (MQ=255)
                                                                                                                                                                                               gggttAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGattagt                                                     <  2:91865/246‑3 (MQ=255)
                                                                                                                                                                                               ggTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTc                                                    <  1:462684/250‑1 (MQ=255)
                                                                                                                                                                                               ggTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTc                                                    <  1:119309/250‑1 (MQ=255)
                                                                                                                                                                                                   tAACCCCCAGGGATTAAACGTAGAAATCCTCGGCTAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTgcg                                                <  2:391796/250‑1 (MQ=255)
                                                                                                                                                                                                        cccAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCttt                                            <  1:311824/249‑1 (MQ=255)
                                                                                                                                                                                                                 tAAACGAAGCAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGAACCATGCTATTTTTGCTGAGCGGTTGATTCAGAACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCg                                 >  2:330983/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaTCCTCGTCGAGGTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGCTGCAAAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCTAGAATGTGATGGAAGAGTACCATGTTAGTTTTGCTGAGCGGATGATTCAGTACCCGGCGGTACAGCTAATCTGCAATACGGTTTATTCTGTGCTTGTTAGGCTAGCGTCGCGTTAAt                       <  2:630658/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGGACCATGCTATTTTTGCTGAGCGGCTGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCgg                    >  2:522855/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCggc                               <  2:556730/233‑2 (MQ=255)
                                                                                                                                                                                                                                     gAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCggc                               >  1:556730/1‑232 (MQ=255)
                                                                                                                                                                                                                                         ttGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTacacggcaccttccccc         >  2:646022/1‑235 (MQ=255)
                                                                                                                                                                                                                                                 gCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTa  <  1:646022/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTA  >  NC_000913/19153‑19641

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