Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,772 T→C F236F (TTT→TTC amiB → N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,7720TC98.3% 190.8 / ‑6.9 58F236F (TTT→TTCamiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/34);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGAT  >  NC_000913/4396541‑4397006
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          
tCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGGACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTg                                                                                                                                                                                                                         >  1:190048/1‑251 (MQ=255)
   gattAAAACGGAAAGGATCCGCCCTACGGGTGTTATCAGCAGGAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACg                                                                                                                                                                                                                      <  2:298535/249‑1 (MQ=255)
     ttAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGGTTTAACCCGTGACggg                                                                                                                                                                                                                    >  1:199146/1‑251 (MQ=255)
        aaaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAAGACGGTAACGCGGCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTTGTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCTGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGAt                                                                                                                                                                                                                 <  1:448534/251‑1 (MQ=255)
            ggAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGGAAAATGTCACCATCGCCATCGCGCGTAAATTTCGTACTTTGCTCCATGAAGATCCGATGTTCCAAGGGGTTTTTACCCGTGACGGGGGTTACt                                                                                                                                                                                                             >  2:223963/1‑251 (MQ=255)
            ggAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGGACGCGGGAGAAAAATTTCACCATCGCCATCGCGCGTAAATTGCGTACTTTTCTCAATGTCGATCCGATGTTCAAAGGCGTTTTACACCGTGACGGGGAATACt                                                                                                                                                                                                             >  1:115689/1‑251 (MQ=255)
               aaGTAACCGCACGCCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGAGCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACGTTTCTTAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTa                                                                                                                                                                                                          <  1:331925/251‑1 (MQ=255)
                    aCCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCCCGCGACGGCTAACACTGGCGATAACATTATCATCGCAATTGATGCCGGACACGGCGGGCAGGAACCTAGCGCTATCGGCACCCGGGGGACGCGGGGGAAAAATGACACCCTCGCCATCGGGCGCAAACTGAGTACATTGCTCAAAGCTGCTCCGAAGTGCGAAGGGGGTTTAACCCAGGACGGGGATTACTTTAtcc                                                                                                                                                                                                       >  2:456435/1‑248 (MQ=255)
                    aCCGCACTACGGGTGTTAGCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGAGAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCg                                                                                                                                                                                                     <  1:304115/251‑1 (MQ=255)
                         aCTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGAt                                                                                                                                                                                                <  1:481024/251‑1 (MQ=255)
                            aCGGGTGTTATCAGCAGTAATACGGTAACGCGGCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTAGTGATGCCGGACACGGCGGTCAGTACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGGTGTAACCCGTGACGGGGATTACTTTATCTCGGTGATgg                                                                                                                                                                                              <  2:203468/250‑1 (MQ=255)
                                      tCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGAGACGGCTAAAACTGACGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGGACCCGGGAGAAAAATTTCACCAACGCCCTCCCGCCGAAATTGAGGACCTTTCTCCATTACGAACCGATGTTTAAAAGCGGTTTTACCCGTGGAGGGGTTTTCGTTATCTCGGGGAAgggggacaggg                                                                                                                                                                                     >  2:337523/1‑242 (MQ=255)
                                                 aCGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTaa                                                                                                                                                                        <  1:223963/251‑1 (MQ=255)
                                                 aCGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTaa                                                                                                                                                                        <  1:425309/251‑1 (MQ=255)
                                                  cGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGGTTACTTTATCTCGGGGATGGGGCGCAGTGATGTGGCACGTa                                                                                                                                                                         >  1:494000/1‑249 (MQ=255)
                                                           gTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTAGTACGCGGGAGGAAAAAGTCACCATCGCCATCGCGCGTAAATTGCGTAATTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGGGATGGGGCGCAGAGATGTGGCACGTAAGCAAAACGcc                                                                                                                                                              >  1:398762/1‑251 (MQ=255)
                                                             ccGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCaa                                                                                                                                                            <  1:187582/251‑1 (MQ=255)
                                                                 cAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTc                                                                                                                                                        <  1:503659/251‑1 (MQ=255)
                                                                       gcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgtg                                                                                                                                                  <  1:523048/251‑1 (MQ=255)
                                                                       gcgcgACGGCTAACACGGGCGATAAAATTATCAGCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgtt                                                                                                                                                  <  1:523060/251‑2 (MQ=255)
                                                                        cgcgactgtgaaCACTGGCTATAAAATAACCATCGCGATTGATGCCGGACACGGCGTTCAGGACCCTGGCGCTATCGGCCCTGGTGGTACGCGGAAGAAAAAGGGCACCAGCGCCATCGCGCGGAAATTGCGTGCTTTGCTCAATGACGAGCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgtgt                                                                                                                                                 <  2:4228/241‑1 (MQ=255)
                                                                         gcgaggtatAACACTGGCGATAAAATTATCATCGCTATTGACGCCGGACACGGGGGTCAGGCCCCTGGCGCTATCGGCCCTGGTGGTACGCGGGAGGAATATGTCACCGTCGCCATCGCGCGTAAATTGCGTGCTTTGCTCAATGACGATCCGATGTTCAAAGGCGTGTTAACCCTTGACGGGGATTACTTTATCTCGGTGATGGGGCGGAGTGAAGTGGCACGGAAGGAAAACGCCAAtgtcaccgtgtc                                                                                                                                                <  2:68323/243‑12 (MQ=255)
                                                                         gcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTc                                                                                                                                                <  1:456089/251‑1 (MQ=255)
                                                                          cgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCg                                                                                                                                               <  1:64465/251‑1 (MQ=255)
                                                                           gACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGGACGCGGGAGAAAAATGACACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATTTTCAAAGGCGTTTTAACCCTTGACGGGGATTACTTTATCTCGGTGATGGGGGGCAGTGATGTGGCACGGTAGCAACACGCCAATTTCCTCGTGTCg                                                                                                                                               >  2:396592/1‑250 (MQ=255)
                                                                                        ggCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGAt                                                                                                                                   <  1:311838/249‑1 (MQ=255)
                                                                                        ggCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGc                                                                                                                                 <  1:478107/251‑1 (MQ=255)
                                                                                                aaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCATCGCCCTCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGGTGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGGGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCg                                                                                                                          >  1:343214/1‑250 (MQ=255)
                                                                                                   gcacagcgctagtgatgctgggatcggcggGCAGGACCCTGGCGCTTACGGCCCCGGGGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGAGCCGATGGTCAAATGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAAc                                                                                                                       <  2:157625/226‑1 (MQ=255)
                                                                                                      catcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCa                                                                                                                   <  1:239150/251‑1 (MQ=255)
                                                                                                            tATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCTGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCg                                                                                                              >  1:144361/1‑250 (MQ=255)
                                                                                                             actgctgcgggACACGGCGGTCGCGTTCCTGGCTCTACCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAAAGTCGATCCGATGTTCAAAGGCGTTTTAGCCCGTGACGGGGATTACTTTAGCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCg                                                                                                              <  2:180787/240‑1 (MQ=255)
                                                                                                                                gTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGAGGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTg                                                                                         <  2:494000/251‑1 (MQ=255)
                                                                                                                                        ccgggaccTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCCCCATCGCCATCGCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCGCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGTGTGCTCTCTaa                                                                                 <  2:287550/244‑1 (MQ=255)
                                                                                                                                        ccTGGCGCTATCGGCACCGGGGGTACGCGGCAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGAACCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTTCGACTGGCGGTTCCGTATGGGTGcatgcta                                                                                  >  2:295840/1‑244 (MQ=255)
                                                                                                                                           tgcgcTAGCGGCCCCTGTGGCACGAGGGAGAAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTCGCTCAATGACGAGCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAGCGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCGCCGACCCGCAGTTCGACTGGCGCTTTCGTATGGGTGCTCTCTAAcc                                                                               <  2:229287/249‑1 (MQ=255)
                                                                                                                                                tttacggctccggcagtacctgggAGAAAAATGTCCCCATCGCCATCGACCGTAAAGGTCGTACTTTGCTCAATGACGATCCGATGGTCAAAGGCGTTTTAACCCGTGTCGGGGATTACGTTATCTCGGTGATGGGGCGCATTTATGGGGCACGTAAGCAAAACGCCAATTTCCGCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCTACTGGCGCTTCCGTATGTGGGTGCTCTAACCGTc                                                                            <  2:327608/228‑1 (MQ=255)
                                                                                                                                                   cGGCCCCGGTGGTACGCGGGAGAAAAATGACACCATCGCCATCGCGCGTAAATTGCGTACTTTGCGCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCa                                                                      <  1:518840/251‑1 (MQ=255)
                                                                                                                                                     ggccccgtaggcaCGCGGGAGAAAAATGGCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGTGGCGCAGTGATGTGGCACGTAAGCAAAACGCCCATTTCATCGTATCGGTTGACGCTGTTGGCGCACCGATCCGACGTGCGCGAGACGCTTCCTTATGGACGCTCTCTCACCGTCGCGCCac                                                                     >  2:540357/13‑250 (MQ=255)
                                                                                                                                                                   gcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCa                                                       <  1:502064/251‑1 (MQ=255)
                                                                                                                                                                           aaaaTGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAACCGCCAATTTCCTCGTGTCGATTCACGCTGCTGCCGCCCCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACCGTGTAATGGCCAgttggctg                                               >  2:76639/1‑251 (MQ=255)
                                                                                                                                                                                      atcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACg                                     >  1:544599/1‑250 (MQ=255)
                                                                                                                                                                                          ccatcgtGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGt                                                                             >  1:139028/1‑206 (MQ=255)
                                                                                                                                                                                          ccatcgtGCGTAAATTGCGTACTTTGCGCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGt                                                                             <  2:139028/206‑1 (MQ=255)
                                                                                                                                                                                              cgcgcgTAAAGTGCGTACTTGGCTCAATGACGCGCCGATGTTCAAAGTCGGTATATCCCGGGACGGGGATTACTGTATCTCGGTGATGGGGCGCAGTGATGTGGCCCGGAAGCAAAACGCCAATTTCCTCGGGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACTAGAAACAg                             <  2:545072/250‑1 (MQ=255)
                                                                                                                                                                                                       aTTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgctgc                   <  2:199146/251‑1 (MQ=255)
                                                                                                                                                                                                        ttGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgctgct                  >  1:196572/1‑251 (MQ=255)
                                                                                                                                                                                                              aCTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACAGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggtgg            <  1:104532/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCgg       >  1:425975/1‑251 (MQ=255)
                                                                                                                                                                                                                   gCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCgg       >  1:318682/1‑251 (MQ=255)
                                                                                                                                                                                                                   gCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCg        <  1:414218/250‑1 (MQ=255)
                                                                                                                                                                                                                     tCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGGTTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAAACGCAGTGCGACTGGCGCTTCCTTATGGGTGCTCTCTAACCGGCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGGACTGCTGTGTGGGGCgggg     >  1:284506/1‑250 (MQ=255)
                                                                                                                                                                                                                      cAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTg    <  1:344403/251‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGCTGGGGCGCAGTGATGTGGCACGTAAGCCAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGACGCACCGAACAGAAGTGCGCCTGGCGCTTCCGTATGGGTGATCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTAGAGCTCAACGAGAAACAGGCGGAGCTGATGGGTgggggggggtg   >  2:23672/1‑248 (MQ=255)
                                                                                                                                                                                                                       aaTGACGATCCGATGTTCAAAGGCGTTTTAACCCGTAACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGCCTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCCGGCAGGCGCACCACGACAAACAATCGGAGCTTGTGAGTGGGGCggggt    >  2:120060/1‑248 (MQ=255)
                                                                                                                                                                                                                        aTGACGATCCGATGTTCAAAGGCGTTTTAACCCGAGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCAAGAGAAACAGTCGGAGCTGCTGGGTGGGGCggggggt  >  1:156815/1‑247 (MQ=255)
                                                                                                                                                                                                                        aTGACGATCCGATGTTCAAAGGCGGTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGATGCTGGGTGGGGCGGGTGAt  <  1:490719/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGCGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTAACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGCGCTCTCTAAcc                                                                               >  1:299281/1‑168 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          
TCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGAT  >  NC_000913/4396541‑4397006

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