Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 91,949 C→T I179I (ATC→ATT ftsI → peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091391,9490CT98.6% 232.6 / ‑4.2 70I179I (ATC→ATTftsIpeptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/35);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

CATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCT  >  NC_000913/91707‑92189
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
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                                                                                                                                                                                                       <  1:145685/251‑1 (MQ=255)
                                   gAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTgaga                                                                                                                                                                                                       <  2:404633/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)
                                                                                                                                                                                        catccgcgggaagaaactcGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTGTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCt                                                   <  2:118154/235‑1 (MQ=255)
                                                                                                                                                                                             gCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGGCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTTCCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGCCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGtt                                               <  2:234444/250‑1 (MQ=255)
                                                                                                                                                                                               gTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCAACTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAt                                            >  1:162506/1‑251 (MQ=255)
                                                                                                                                                                                                   agaatctcGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGCTTGAGGGCGTTGAGAAGAGTTTCGATACATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACAGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGAAAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGCTGAACGCCTGCCGGCGCTGGTTTATcgcg                                        >  1:394413/5‑251 (MQ=255)
                                                                                                                                                                                                          cGCCGTTATTATCCATCCGGCGAAGTGACTGCTCTCCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGATTTGCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACt                                    <  1:41876/248‑1 (MQ=255)
                                                                                                                                                                                                          cGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaa                                 <  1:472070/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGATAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaac                                >  1:493470/1‑250 (MQ=255)
                                                                                                                                                                                                            ccGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaac                                >  2:94608/1‑250 (MQ=255)
                                                                                                                                                                                                                    atccatccGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTg                       >  1:270520/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccatccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGc                     >  1:304079/4‑251 (MQ=255)
                                                                                                                                                                                                                        atccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCt                   >  1:6434/2‑251 (MQ=255)
                                                                                                                                                                                                                         tccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCtt                  >  1:292145/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGAGATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCttt                 <  1:76130/251‑1 (MQ=255)
                                                                                                                                                                                                                           cagCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCttt                 >  2:410580/3‑250 (MQ=255)
                                                                                                                                                                                                                            ggCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGGTATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGCACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTa                >  2:70131/1‑250 (MQ=255)
                                                                                                                                                                                                                              cGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAAc              >  1:292636/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gTGACTGCACACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGCTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGc         <  2:498136/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCt        >  2:176692/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGGCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCt        >  2:57807/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cACCTCATTGTCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAAGTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCt  <  1:313031/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
CATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCT  >  NC_000913/91707‑92189

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