Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,546,8040TA79.5% 93.7 / 18.5 45D44E (GAT→GAA) murPN‑acetylmuramic acid‑specific PTS enzyme IICB component/anhydro‑N‑acetylmuramic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base A (15/20);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAACTGTAAAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTGATC  >  NC_000913/2546558‑2547047
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
cAACTGTAAAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCa                                                                                                                                                                                                                                                  <  1:102568/251‑1 (MQ=255)
 aaCTGTAAAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCa                                                                                                                                                                                                                                                  >  2:405309/1‑250 (MQ=255)
     gTAAAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGGGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAatat                                                                                                                                                                                                                                              <  2:27566/250‑1 (MQ=255)
               cATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTgg                                                                                                                                                                                                                                   >  2:224172/1‑251 (MQ=255)
                  ggtgATGGTGCTGAAAAATCGCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCTTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGTCAGTTCACTGGTTGATCCCAATATCAAAACGCGGGAAg                                                                                                                                                                                                                                <  2:318636/250‑1 (MQ=255)
                     gATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGc                                                                                                                                                                                                                              <  1:250698/250‑1 (MQ=255)
                     gATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGc                                                                                                                                                                                                                              <  1:224172/250‑1 (MQ=255)
                       tGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGAGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGGGTACATGACAGATCACTGGTTGATCCCAATATCACAACACTGGAAGGCGTg                                                                                                                                                                                                                           >  1:61011/1‑251 (MQ=255)
                                               cAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGACCACCATTCTTCCCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGCCAGTTCACTGGTTGACCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGAc                                                                                                                                                                                                   >  1:41565/1‑251 (MQ=255)
                                                                  ggATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTg                                                                                                                                                                                <  2:202716/251‑1 (MQ=255)
                                                                      cAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGcc                                                                                                                                                                            >  1:119598/1‑251 (MQ=255)
                                                                                  tttATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAgcgc                                                                                                                                                                 >  2:192223/1‑250 (MQ=255)
                                                                                     aTTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCg                                                                                                                                                             >  1:429126/1‑251 (MQ=255)
                                                                                                       aaGGAATAAACCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgag                                                                                                                                           <  2:179733/251‑1 (MQ=255)
                                                                                                             tAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCAGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGCGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCTCAAAAGCGATGAGCGAgct                                                                                                                                     >  1:342309/1‑251 (MQ=255)
                                                                                                                 ccATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTgg                                                                                                                                                                         >  2:274216/1‑211 (MQ=255)
                                                                                                                 ccATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTgg                                                                                                                                                                         <  1:274216/211‑1 (MQ=255)
                                                                                                                             agaTCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCg                                                                                                                    <  2:231863/251‑1 (MQ=255)
                                                                                                                             agaTCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGAGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCAc                                                                                                                      >  2:142730/1‑249 (MQ=255)
                                                                                                                             agaNCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCAGTTGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGATCAGCGACCAGGTGCAGGTCGTTTTTTGGCCTGGGAAAGCGCATCGTTCCGCAAAAGCGATGAGCGAGCTGTTGGG‑AGAAGCACCg                                                                                                                    <  2:440941/251‑1 (MQ=255)
                                                                                                                                        gagcTTCTGAACACCATTCTTACCCGCGTGGGGGGACCGGGAAATATCGCCAGTTGTGGTAACTGTAGGACGCGCCTGCGTATGGGTGTACATGACACTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGGCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCAGCGCGCCGCAAAAGCGATGAGCGAGCTGCGGTG‑AGAAGCACCGGTACAGGATg                                                                                                          <  2:63374/247‑1 (MQ=255)
                                                                                                                                         aaCTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGcc                                                                                                        >  1:197716/1‑251 (MQ=255)
                                                                                                                                          aCTTCTGAACCCCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCg                                                                                                       <  2:168463/251‑1 (MQ=255)
                                                                                                                                          aCTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTTCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCg                                                                                                       <  2:295229/251‑1 (MQ=255)
                                                                                                                                              cTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGc                                                                                                      >  2:116694/1‑248 (MQ=255)
                                                                                                                                                       aTTCTTATCCGTGTCGGCGAACCGGAAAATATCGCAAGTTGTGGTAAATGTCTGCCGCGCCTGCGTCTGGGTGTACATGACAGTTAACTGGATGAACCCCATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAACGCGATGAGCGAGCAGCTGGGCAGA‑GCACCGGTACAGGCTGCCGCAGAAATcgccgc                                                                                          >  2:489917/1‑251 (MQ=255)
                                                                                                                                                         tCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGccc                                                                                        >  1:442677/1‑251 (MQ=255)
                                                                                                                                                            tACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAg                                                                                      <  2:197716/250‑1 (MQ=255)
                                                                                                                                                                 gtgtCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGGCGTTTTTGGGCCTGGTAAAGAGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGAAGAAATCGCCGCCCAGAACaa                                                                                 >  1:214975/1‑250 (MQ=255)
                                                                                                                                                                   gtCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAAc                                                                               <  1:405309/250‑1 (MQ=255)
                                                                                                                                                                   gtCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGAGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCAAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACg                                                                              >  2:366338/1‑251 (MQ=255)
                                                                                                                                                                    gcggcggACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCAGTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCAGCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGt                                                                             <  2:41565/250‑1 (MQ=255)
                                                                                                                                                                      ggcggACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTc                                                                            >  1:321666/1‑250 (MQ=255)
                                                                                                                                                                                          cAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAAGCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCt                                                       <  2:206531/251‑1 (MQ=255)
                                                                                                                                                                                                    aaCTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGca                                             <  2:30966/251‑1 (MQ=255)
                                                                                                                                                                                                          aTGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAAtt                                       >  1:429465/1‑251 (MQ=255)
                                                                                                                                                                                                          aTGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAAtt                                       <  1:226212/251‑1 (MQ=255)
                                                                                                                                                                                                              cgcgCCTGCGTCTGGGTGTTCATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGTCGTTAACGTCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATCAGCCAGCTGCTGGG‑AGAATCACCGTTACAGGTTGCCGTTGAAATCGCCCCCCATAGCAACCGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTtctc                                   <  1:437516/251‑1 (MQ=255)
                                                                                                                                                                                                                 gcCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTTGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGcc                                >  1:432530/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAAt                            <  1:25666/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTAATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCCACAATTTCTCGCCaaa                             <  2:342309/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCtt               <  2:429126/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGTGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCGAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCTCAGAAATCGCCGCCCAGAACAGACGTCAGTTGATAGCCAAACAGACCTCCGGCGTGCAACAAGTTCTCGCCGAAATCGCGACCATCTTc              <  2:203532/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCAcgccg        >  2:259408/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ggTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTGATc  <  2:394862/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CAACTGTAAAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGG‑AGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTGATC  >  NC_000913/2546558‑2547047

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