Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 92,976 A→G 5.4% S522G (AGT→GGT)  ftsI → transpeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091392,9760AG 5.4% 299.0 / 4.5 93S522G (AGT→GGT) ftsItranspeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base A (45/43);  new base G (3/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACG  >  NC_000913/92828‑93122
                                                                                                                                                    |                                                                                                                                                  
aaGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                                                                    <  1:17572/149‑1 (MQ=255)
       gcgcTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGcc                                                                                                                                             <  1:254363/149‑1 (MQ=255)
         gcTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcgc                                                                                                                                           >  1:332803/1‑149 (MQ=255)
             ccAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTCCCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCtt                                                                                                                                       >  5:37099/1‑149 (MQ=255)
              ctggcggcggcggcGTGAAGGCGGCGATTAAAGGCTTTAGTATCGCCATTAAAACCGGTAACGCGATAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATAACGCAGGCGTTGCGCCTGAGAGTCAGCCGCGCTTc                                                                                                                                      <  6:377750/147‑1 (MQ=255)
               aGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAACACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATAATGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTcg                                                                                                                                     >  1:145354/1‑149 (MQ=255)
                     cgtcggctTGAAGGCGTCTATTAATGGCTATCGTATCGCAATTAAAACCGGTACCTCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTgg                                                                                                                               <  4:268640/146‑1 (MQ=255)
                       gcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgtt                                                                                                                             >  4:273872/1‑149 (MQ=255)
                          gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCTATTAAAACCGGTACCGCGAAAAAGGGCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                          <  8:184704/149‑1 (MQ=255)
                          gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                          <  8:114466/149‑1 (MQ=255)
                          gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                          <  5:285214/149‑1 (MQ=255)
                          gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                          <  3:309036/149‑1 (MQ=255)
                          gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                          <  5:63221/149‑1 (MQ=255)
                           cgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttg                                                                                                                         <  7:151277/149‑1 (MQ=255)
                             tGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgtt                                                                                                                       <  6:111494/149‑1 (MQ=255)
                               aaGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTAt                                                                                                                     >  7:106735/1‑149 (MQ=255)
                                  gcggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                                  >  2:209760/1‑149 (MQ=255)
                                  gcggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                                  >  2:87219/1‑149 (MQ=255)
                                  gcggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                                  >  3:22751/1‑149 (MQ=255)
                                  gcggcgATTAAAGGCAATCGTATCGCCATTAAAACCGGTACCGCGCAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATCTTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                                  <  6:200755/149‑1 (MQ=255)
                                   cggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAAc                                                                                                                 >  2:132936/1‑149 (MQ=255)
                                       gATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                                                                                             <  6:161745/149‑1 (MQ=255)
                                         ttAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGGTGGTTGTTGTTATCAACGATCCg                                                                                                           >  5:264409/1‑149 (MQ=255)
                                                 tatcgtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATc                                                                                                                    <  6:339292/132‑1 (MQ=255)
                                                 tatcgtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGGTGTTGTTATc                                                                                                                    >  5:339292/1‑132 (MQ=255)
                                                   tcgtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGt                                                                                                  >  2:16647/1‑148 (MQ=255)
                                                     gttgcgCCATTAAAACCGGTAGCGCGAAAAAGGTCGGGCCGGACGGACTCTACATCAATAAATATATTGCTTATACCGCAGGCGTGGCGCATGAGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaa                                                                                               <  1:171197/145‑1 (MQ=255)
                                                     gtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaa                                                                                               <  3:273872/149‑1 (MQ=255)
                                                      tatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaa                                                                                               >  1:253634/1‑148 (MQ=255)
                                                          gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCgg                                                                                                    >  6:234074/1‑139 (MQ=255)
                                                          gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCgg                                                                                                    <  5:234074/139‑1 (MQ=255)
                                                          gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtacta                                                                                          >  1:292443/1‑149 (MQ=255)
                                                           ccATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAt                                                                                              >  6:213938/1‑144 (MQ=255)
                                                           ccATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAt                                                                                              <  5:213938/144‑1 (MQ=255)
                                                                 aatccGGTACCGCGAAAAAGGTCGTGCCTGACGGTCGCTACATCAATAAATATATTGCTTATACCGCTGGCGTTGAGCCTTCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCATCGTTCCGCAGGCGGGTAAATACTACggcggc                                                                                   <  6:68091/146‑1 (MQ=255)
                                                                 aaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAta                                                                                             <  7:349818/139‑1 (MQ=255)
                                                                 aaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAta                                                                                             >  8:349818/1‑139 (MQ=255)
                                                                 aaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACAGCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc                                                                                   >  2:122983/1‑149 (MQ=255)
                                                                      ggTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGt                                                                              >  1:3435/1‑149 (MQ=255)
                                                                       gTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcgcg                                                                                                                                          >  2:190903/1‑88 (MQ=255)
                                                                       gTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcgcg                                                                                                                                          <  1:190903/88‑1 (MQ=255)
                                                                          ccGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTcc                                                                          >  7:283690/1‑149 (MQ=255)
                                                                           cgcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCgg                                                                                                    >  2:360049/1‑122 (MQ=255)
                                                                           cgcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCgg                                                                                                    <  1:360049/122‑1 (MQ=255)
                                                                             cgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAAGGCGGCGCCGTTTCcgcg                                                                       >  1:151225/1‑149 (MQ=255)
                                                                              gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                                      >  3:358397/1‑149 (MQ=255)
                                                                               aaaaaGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCGTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGcc                                                                     <  7:33077/149‑1 (MQ=255)
                                                                                aaaaGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCg                                                                    <  7:308790/149‑1 (MQ=255)
                                                                                       cGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTg                                                             <  6:37099/149‑1 (MQ=255)
                                                                                        ggggCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCCGCGCCGTTTCCGCGCCGGTCTTTg                                                             <  5:175459/148‑1 (MQ=255)
                                                                                        gggCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTgg                                                            <  5:179197/149‑1 (MQ=255)
                                                                                          gCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTg                                                          >  7:48320/1‑149 (MQ=255)
                                                                                          gCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTg                                                          >  7:30106/1‑149 (MQ=255)
                                                                                          gCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTg                                                          >  6:361912/1‑149 (MQ=255)
                                                                                              gACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCCGGTAAATACTACGGCGGCGCCGTTTCCGCGCTGGTCTTTGGTGCcat                                                      >  1:368417/1‑149 (MQ=255)
                                                                                                     gCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTAAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCTCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATggg                                                <  1:204967/148‑1 (MQ=255)
                                                                                                      cTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcg                                              <  8:283690/149‑1 (MQ=255)
                                                                                                         cATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGgcggcg                                           >  5:126705/1‑149 (MQ=255)
                                                                                                            cAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcggc                                            >  6:372232/1‑145 (MQ=255)
                                                                                                            cAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcggc                                            <  5:372232/145‑1 (MQ=255)
                                                                                                             aataaataTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTAtt                                       <  4:22751/149‑1 (MQ=255)
                                                                                                                  atatatTGCTTATACCGCAGGCGTTGCGCCGGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAAAGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTa                                  <  6:136987/149‑1 (MQ=255)
                                                                                                                   tatatTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGGCTTTGGTGCCATCATGGGCGGCGTATTGCGTAc                                 <  1:209760/149‑1 (MQ=255)
                                                                                                                    atatTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcgg                                                                                       >  8:211914/1‑94 (MQ=255)
                                                                                                                    atatTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcgg                                                                                       <  7:211914/94‑1 (MQ=255)
                                                                                                                      atTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCAt                              >  2:368575/1‑149 (MQ=255)
                                                                                                                       ttGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGAATTCCGTACCAt                              >  6:375113/1‑148 (MQ=255)
                                                                                                                        tGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATg                             >  4:337494/1‑148 (MQ=255)
                                                                                                                          cTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGaa                           >  6:250937/1‑148 (MQ=255)
                                                                                                                           ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                                      <  2:71516/104‑1 (MQ=255)
                                                                                                                           ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                                      >  1:71516/1‑104 (MQ=255)
                                                                                                                              taCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATc                       <  2:332803/148‑1 (MQ=255)
                                                                                                                                tcGCAGGCGCAACGCCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAg                    >  2:223919/2‑149 (MQ=255)
                                                                                                                                tcGCAGGCGCAACGCCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAg                    >  6:1788/2‑149 (MQ=255)
                                                                                                                                ccGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGa                     <  1:132936/148‑1 (MQ=255)
                                                                                                                                ccGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAg                    >  1:294132/1‑149 (MQ=255)
                                                                                                                                  gCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGcc                  <  4:358397/149‑1 (MQ=255)
                                                                                                                                  gCAGGCGCAACGCCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGcc                  >  8:149071/1‑149 (MQ=255)
                                                                                                                                   cAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTCCGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTCCCATGAACATCGAGCCg                 >  3:29778/1‑149 (MQ=255)
                                                                                                                                    aGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCATTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCgg                <  8:106735/149‑1 (MQ=255)
                                                                                                                                     ggCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGa               <  1:87219/149‑1 (MQ=255)
                                                                                                                                      gCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGAt              <  8:48320/149‑1 (MQ=255)
                                                                                                                                        gTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATgc            >  4:359718/1‑149 (MQ=255)
                                                                                                                                         ttGCGCGTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATgcg           >  8:140386/1‑149 (MQ=255)
                                                                                                                                           gcgcCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCTCAGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATgcgc          <  2:292443/148‑1 (MQ=255)
                                                                                                                                           acgcCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCt         <  5:1788/148‑1 (MQ=255)
                                                                                                                                           acgcCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCt         <  1:223919/148‑1 (MQ=255)
                                                                                                                                            cgcCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCt         >  4:145104/1‑148 (MQ=255)
                                                                                                                                             gcCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGa       >  6:159444/1‑149 (MQ=255)
                                                                                                                                             gcCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGa       >  2:301149/1‑149 (MQ=255)
                                                                                                                                                tGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACaa    >  1:83308/1‑149 (MQ=255)
                                                                                                                                                 gCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGc                                     <  4:106108/115‑1 (MQ=255)
                                                                                                                                                 gCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGc                                     >  3:106108/1‑115 (MQ=255)
                                                                                                                                                  cGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACg  <  1:277327/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACG  >  NC_000913/92828‑93122

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