Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 91,904 G→A E164E (GAG→GAA ftsI → peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091391,9040GA95.3% 196.0 / 1.0 65E164E (GAG→GAAftsIpeptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (29/32);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

TAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAT  >  NC_000913/91657‑92147
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
gaGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAAtct                                                                                                                                                                                                                                                  <  1:562177/250‑1 (MQ=255)
               cGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTAtcc                                                                                                                                                                                                                                    >  1:481644/1‑250 (MQ=255)
                   aaaaGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTTTTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCCTCTGCGTGAAGAATCTCGCCGTTAATAtccatccg                                                                                                                                                                                                                               >  1:146220/1‑251 (MQ=255)
                    aaaGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGCATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTCTTTATCTGGCGCGTCAGGTGAACCCTGACCTGGCCGAATACATTCAAAAACTGAAACTGCCGGTGATTCCTCTGACTGAAGAATCTCGCCGGTATTATCCATCCgg                                                                                                                                                                                                                              >  2:189311/1‑251 (MQ=255)
                    aaaGCGATTTAGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGGCGGTGGCCGCTGGGTGGCGCTGGCTAACGCGCCCCATATTCCGCTGGAACAAGCTTCAGCCCGCATTAACTCCAACCCCGAAGGGCGCTTTATTTCTCTGGCGCGGCAGGGGAACCCTGGCCTGGCGGACTGCATACAAAAACTGAAACTGGCGGGGGTTCATCTGCGGGAAAAATCTCCCCGTTCTTATCCATCCgg                                                                                                                                                                                                                              >  2:297977/1‑251 (MQ=255)
                    aaaGCGATTGGGGCGGACCGGAAAGAAGGGCATGAGGCTGGCGGTATCAGCGGCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCgg                                                                                                                                                                                                                              <  2:84945/251‑1 (MQ=255)
                      aGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGCATCCGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGGCTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTAGCCGTTTTTATCCATCCGGCg                                                                                                                                                                                                                            >  1:234444/1‑251 (MQ=255)
                         gATTTGGGCTGACCCGAAAGAAGTGCATGCCGCTGGCGGTACCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGTATCAGCTTTCAGCCCGAATTAACGCCAACCCGAAAGGGCGCTTTATTTATCGGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAg                                                                                                                                                                                                                         <  2:353624/251‑1 (MQ=255)
                         gATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGc                                                                                                                                                                                                                                               >  2:152545/1‑228 (MQ=255)
                         gATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGc                                                                                                                                                                                                                                               <  1:152545/228‑1 (MQ=255)
                          aTTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAAGCGCTGGCTCACGCGCTCAATATTCCGCTGGATCCGCTTTCAGCCCGCATTAACGCCAAACCGACAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTattttcaatccggcgacgg                                                                                                                                                                                                                        >  1:106266/1‑233 (MQ=255)
                                   gACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAGACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTc                                                                                                                                                                                                                <  2:230161/250‑1 (MQ=255)
                                      ccgaaagaattgcatatcgcgttcGGTACCAGCGTAGGTGCCCCCTGGAAGGCGCTGGCTAACTCGCTCGAGATTCCGCTGGATCAGCTTTCACCCCCCATTAACGCCAACCCGAAAGGGCCCTTTATTTATCTGGCGCGTCAGGTTAACCCTGACATGGCGGACTACAGCAAAATACTGAAACTGCCGGGGATACATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCt                                                                                                                                                                                                            <  2:565596/228‑1 (MQ=255)
                                                  cATGACGCTGGCGGTATCAGCGTAGGTGACCGCTGGAAGGCGCTGGCTAAAGCGCTCAATATTCCGCTGGATCCACTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTTTCTGGCGCGTCAGGTGAACCCTGACATGGCGGAATACATCAAAAAACTGAAACGGCCGGGGATTCATCTGCGTGTAGAATCTCGCCGTTATTATCCATCCGGCGAAGGGACTGCTCCCCCCATTGGCTTTa                                                                                                                                                                                                 >  2:562122/1‑250 (MQ=255)
                                                      acggtggtggTAACAGCGGCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAAc                                                                                                                                                                                            <  2:63005/243‑1 (MQ=255)
                                                      aaggtGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCTCTGGCTAACGCTCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAAc                                                                                                                                                                                            <  2:397174/247‑1 (MQ=255)
                                                           ggCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCg                                                                                                                                                                                        <  1:220349/250‑1 (MQ=255)
                                                                 aTCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTc                                                                                                                                                                                  <  1:312145/250‑1 (MQ=255)
                                                                          ggTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTgg                                                                                                                                                                        >  1:218932/1‑250 (MQ=255)
                                                                              aCCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGc                                                                                                                                                                    >  1:221321/1‑251 (MQ=255)
                                                                                     gAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTgaga                                                                                                                                                             <  2:404633/251‑1 (MQ=255)
                                                                                     gAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTgaga                                                                                                                                                             <  1:145685/251‑1 (MQ=255)
                                                                                           gcTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGtt                                                                                                                                                       >  1:258757/1‑251 (MQ=255)
                                                                                              ggCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGAGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCg                                                                                                                                                    >  1:531397/1‑251 (MQ=255)
                                                                                                  aaCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATa                                                                                                                                                 <  2:221321/250‑1 (MQ=255)
                                                                                                    cgcgCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAAt                                                                                                                                              <  1:53678/251‑1 (MQ=255)
                                                                                                       gcTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACAGATCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGc                                                                                                                                           >  1:498136/1‑251 (MQ=255)
                                                                                                        cTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCt                                                                                                                                          >  1:421208/1‑251 (MQ=255)
                                                                                                         tCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCtt                                                                                                                                         >  2:471981/1‑251 (MQ=255)
                                                                                                          cAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTAAGGTGAACCCTGACAGGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCGGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCAGTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCtt                                                                                                                                         <  2:211395/250‑1 (MQ=255)
                                                                                                                  cGCTGGAGCAGCTTGCAGCCCGCATTAACGCCAACCAGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGGGAACCCTGACATGGCGTACGACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAg                                                                                                                                <  2:218932/251‑1 (MQ=255)
                                                                                                                   gCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGc                                                                                                                               <  1:521073/251‑1 (MQ=255)
                                                                                                                     tGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCg                                                                                                                             >  1:184360/1‑251 (MQ=255)
                                                                                                                             gCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAgc                                                                                                                      <  1:401958/251‑1 (MQ=255)
                                                                                                                                  cAGCCCGCATTAACGCCAACCCGAAAGGCCGTTTTATTTATCTGGCGCGTCAGGGGAACCCTGACATGTCGGACTACATCAAAAAACTGAATCTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCTTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAgcgcttt                                                                                                                 <  2:360016/251‑4 (MQ=255)
                                                                                                                                   aGCCCGCATTAACGCCATCCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATtg                                                                                                                <  2:473114/251‑1 (MQ=255)
                                                                                                                                                gcaaccccGAAAGGGCCCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAACCTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGATGAGTTTCGAGAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTTTGCGTAAAGCCCg                                                                                                   <  2:124161/246‑1 (MQ=255)
                                                                                                                                                gCCATCCCGAAAGGGCTCTTTATTTAGCTGGTGCGTCAGTTGACCCCGGACATGGCGGACTACAGCAAAAAACTGAAACTGCCTGGGATTCATCTGCGTGAATAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCg                                                                                                   <  2:151940/251‑1 (MQ=255)
                                                                                                                                                       cGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGt                                                                                            >  1:552613/1‑251 (MQ=255)
                                                                                                                                                         aaaGGCCGCTTTATTTACCTGGCTCGTCAGGGGAACCCTGACATGGCGGACTACATCAAAAAACTGAACCTTCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTcg                                                                                          <  2:106108/251‑1 (MQ=255)
                                                                                                                                                           aGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTcgcg                                                                                        <  1:183132/251‑1 (MQ=255)
                                                                                                                                                           aGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTcgcg                                                                                        >  2:274620/1‑251 (MQ=255)
                                                                                                                                                             ggCGCTTTATTTATCGGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTa                                                                                      <  2:258757/251‑1 (MQ=255)
                                                                                                                                                               cgcTTGATTTATCTGGCGCGTGAGGTCAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTCTCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGAAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAAt                                                                                    <  2:46733/251‑1 (MQ=255)
                                                                                                                                                                      tttatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACACCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGa                                                                              <  2:184360/250‑1 (MQ=255)
                                                                                                                                                                       ttatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGat                                                                             >  1:501778/1‑250 (MQ=255)
                                                                                                                                                                        tatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGatat                                                                           >  1:319532/1‑251 (MQ=255)
                                                                                                                                                                          tcagtggggtatggTGACCCCTGACATGCGGGACTCCAACAAAAACTGTAACCGTCCGGGTATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCTCCTCAATGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAACCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATAtt                                                                          <  2:115636/238‑1 (MQ=255)
                                                                                                                                                                           cTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATttc                                                                        >  1:139255/1‑251 (MQ=255)
                                                                                                                                                                                      gttcagcccTGACAGGGCGGATTACACTAAAAAACGCCAACTGCCGTGGTTTCTTCGACGTGAAGATTGGCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGATAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACa                                                             <  2:305676/245‑1 (MQ=255)
                                                                                                                                                                                           aCCCTGACATGGCGGACTACATAAAAAAACGGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAg                                                        <  2:146220/251‑1 (MQ=255)
                                                                                                                                                                                                 aCATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTgcg                                                  >  2:206906/1‑251 (MQ=255)
                                                                                                                                                                                                   aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGCGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAGGATATTTCTTCTACTGACAGCCAGGCTgcgcc                                                >  2:185969/1‑250 (MQ=255)
                                                                                                                                                                                                       cGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAAc                                            >  1:436966/1‑251 (MQ=255)
                                                                                                                                                                                                        ggACTACAGCAAAAAACTGAAACTTCCGGGTATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCGCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAAcc                                           <  1:376435/251‑1 (MQ=255)
                                                                                                                                                                                                              cATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGgcgc                                     >  2:24728/1‑251 (MQ=255)
                                                                                                                                                                                                                       actgaaactgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATtg                            >  1:561786/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggaaacggCCGGGGATTCATCTGCTGGAAGAATCTCGCCGTTATTATCCAGCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATtgat                          <  2:112821/250‑1 (MQ=255)
                                                                                                                                                                                                                           aaactgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATtgatga                        <  2:207250/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGCTTAACGCCTg                 >  2:44450/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ggggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGc                >  1:173583/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggATTCATCTGCGTGAAGAAGCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAg              <  2:499346/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGgc            <  1:297977/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGGCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTTCCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGCCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGtt     <  2:234444/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCAACTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAt  >  1:162506/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
TAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAT  >  NC_000913/91657‑92147

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