Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,529,039 C→T G402D (GGT→GAT)  sgcC ← putative PTS enzyme IIC component SgcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,0390CT96.3% 151.6 / NA 55G402D (GGT→GAT) sgcCputative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (34/18);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCC  >  NC_000913/4528796‑4529281
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ttGTCTATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACat                                                                                                                                                                                                                                               >  1:198371/1‑251 (MQ=255)
               ttCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGCCAACGCCAATAATCCCAAGGTACATCAATTCAGAAATCGACCATTCGATCATATTGCCCACGGACa                                                                                                                                                                                                                                >  2:221308/1‑250 (MQ=255)
                 tcaacccAACCTGCGCACCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGTTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACaa                                                                                                                                                                                                                               <  1:291463/246‑1 (MQ=255)
                               gcATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTAATTGAAAGCACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGACACGACAACGCCAATAATCCCAAGCGACATAAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACACCGCTGAGCTTTgcgc                                                                                                                                                                                                                >  1:358576/1‑251 (MQ=255)
                                 aTCAACACTCGGGTCGCCGGGCAATGCACGTAGATGACACATGGCGATCACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTAACGTGAGAACAGCACGACGCAGGCGACAATACCGACCGCCACGACAACGCCAAAAATCACAAGAGACATCAATTCAGAAATCGACCAATCGACACTATTGCCCACGGACAACGCTTAGAttttggcgc                                                                                                                                                                                                              >  2:5567/1‑245 (MQ=255)
                                   cAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATTGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCAAGCCGCTGGCGCCAATACCGCCCGCCACGACACCGCCAATAATCCCAAGCGCCATCAATTCAGGAATCTACCAATCGAACATATTGCACCCGGCCAACGCTGAGTTGTGCGCACtt                                                                                                                                                                                                            >  2:552017/1‑251 (MQ=255)
                                                cccGGGCAATGCACGAAGATGGCACATGGCGATGACCGCTTTTCCCGTCCCAATACCCTCTTCCAGCCAACCCATGAAAAGTACCCCTTATCCCGATAATTCCCTTTAACGAAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAAc                                                                                                                                                                                               <  1:283376/250‑1 (MQ=255)
                                                     caaaTGCTCGTATATGGCACAGGGGTAGGACCGTTTTTCCCGCCCAAAAACCCCCTTTCAGCCAACACATGGAAATTACTCCTTATTCAAATAATTCACGTTTACGTAAGAACAGCACGCCGCTGGCGACAATACCTACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGaaa                                                                                                                                                                                          <  2:115823/249‑1 (MQ=255)
                                                                  aTGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGGACTCCTTATTCCGTTAATTCACGTTTACGTCAGAACAGCACGACGCTGGCGACAAAAACGACCACCACGACAACGCCAATAAGCCAAAGCGACCTCAATTCAGAACTCGACCAATCGAACATATTGCCCCAGGACAACGACGAGATTGGGGCACTTTCACTGACAAAACTGAAACCGCCGTTgga                                                                                                                                                                               >  2:92906/1‑247 (MQ=255)
                                                                                 gACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGt                                                                                                                                                              >  1:248259/1‑251 (MQ=255)
                                                                                         tttCCGTCCCAATAACCTCTTTCAGCCACCTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGAACAACGCTGAGATTTGCGCACTTTCACCGGCAAACCTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGcgcg                                                                                                                                                      >  1:129541/1‑251 (MQ=255)
                                                                                              gTCCCAATAAACTCTTTCAACCAACTCATTGAAAGTACTACTTATTCCGATAATTCACGTTTACGGAAGAACAGCACGACGCTGGCGACAATACAGACCGCCACGCCAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTTCGCACTTTCACCGGCAAAACTGAAACCGCCTTTAAGCGCCATTTCAGGGAAGTACGGCGCGAACTg                                                                                                                                                 >  2:213040/1‑251 (MQ=255)
                                                                                               tCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGAAAGGCGCGAACTgg                                                                                                                                                >  2:178649/1‑251 (MQ=255)
                                                                                                 ccAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTTTTCGGAAAATTCCCGTTAAGGAAGAACAGCCCCAACCCGGCGAAAATACCGACCGCCACGACAACGCCAATAATCCCCAGCGACATCAATTCAGAAATCGACCAAACGAACATATTGCCCACGGACAACGCTGAGATTTGGGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGGGGAGTAAGGCGCGAACTGGGTg                                                                                                                                             >  1:519268/1‑251 (MQ=255)
                                                                                                   ggtgtgctcttccgatcaaccTCATTGAAAGTCCTCCTTATCCCAAAAATCCCCGTTAACTAAAAAACACCCCGACGCGGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTg                                                                                                                                             <  1:31104/231‑1 (MQ=255)
                                                                                                         cccttctcaccaacagcatgaaaaGTACTCTTATTCCCAAAAATTCACTTTTAGATAAGAACACCACGCCGCTGGAGAAAATCCCGACCGCCACGACAACGCCAATAACCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                       <  1:50711/230‑1 (MQ=255)
                                                                                                           ctTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGCTAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATCCAGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                         >  2:554633/1‑182 (MQ=255)
                                                                                                           ctTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATcagc                                                                                                                                    >  2:223523/1‑251 (MQ=255)
                                                                                                           ctTTCAGCCAACTCATTGAAAGTACTCCTTATCCCGATAATTCACGTTTACGAAAGAACAGCACGACGCTGGCGACAAAACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                         <  1:554633/182‑1 (MQ=255)
                                                                                                                cAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAgg                                                                                                                                 >  1:517443/1‑250 (MQ=255)
                                                                                                                     aaCTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGCCAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTgagatcggaagagagtcgg                                                                                                                           >  2:31104/1‑233 (MQ=255)
                                                                                                                      ccTCATTAAAAGTCCTCCTTATCCCGAAAATCCCCGTTAACATAAAAACACCACGACCCGGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTaa                                                                                                                          <  1:413618/250‑1 (MQ=255)
                                                                                                                          aTTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAAAAGCACGACGCTGGCGACAATACCGACCGCCACCACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCACACGGACAACGCTGAGATTTGAGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGCAGTAAGGCGCGAACTGGGAGGCGATCAGCAGGACGGTAatcat                                                                                                                      >  1:31453/1‑251 (MQ=255)
                                                                                                                                  tACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACACTACCGACCGCCCCGCCAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATGATCACGATg                                                                                                               >  2:70543/1‑250 (MQ=255)
                                                                                                                                            tCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTAt                                                                                                    >  1:237761/1‑251 (MQ=255)
                                                                                                                                                 aTAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAAc                                                                                                >  1:171333/1‑251 (MQ=255)
                                                                                                                                                     ttCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGGAATCATCACGCTGACCCCGCTTATAAACGTcc                                                                                            >  2:194478/1‑251 (MQ=255)
                                                                                                                                                      tCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCt                                                                                           >  1:547421/1‑251 (MQ=255)
                                                                                                                                                             tACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAgg                                                                                    >  1:451260/1‑251 (MQ=255)
                                                                                                                                                                 tAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGtctc                                                                                >  1:473904/1‑251 (MQ=255)
                                                                                                                                                                     aaaaGCACGCCGCGGGCAACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCATGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCcgcg                                                                            <  1:194478/248‑1 (MQ=255)
                                                                                                                                                                       cAGCACGACGCTGGCGCCAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGCTTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAATAAGGCGCGAACTGGGTGGCGATAAGCAGGACGGGAATCATCACGATGACCCCGCTTATCAACGTCTTGCCCAGTTCTCCACgg                                                                           >  2:479193/1‑250 (MQ=255)
                                                                                                                                                                        aGCACGCCCCTGGCAACAATCCCCACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTg                                                                         <  1:433338/251‑1 (MQ=255)
                                                                                                                                                                           acgacgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAACTCGCCCAATCGAACATATTGCCCACGGACAACGCTGAGAGTTGCGCACTTTCACCGGCAAAACTGAACACGCCATTGAGCGCCATTTCAATGAAGTAAGGCGCGAACTGGGTGGCGATCATCAGGACGGTAATCAACACGCTGGCCCCGCTTATCACCGCCCTGCACAGGTCTACGAGGTgtat                                                                      >  2:136555/1‑248 (MQ=255)
                                                                                                                                                                              acgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAc                                                                   >  1:538192/1‑251 (MQ=255)
                                                                                                                                                                                         aaccccgaccgCCACGACAACGCAAATAATCCCAAGCGACATCAATACAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACAGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATAc                                                        <  1:419473/247‑1 (MQ=255)
                                                                                                                                                                                                       cgccAACGCCAATAATCCCAAGCGCCATCAATTGAGAAACCGACCAATCGAACAATGGGCATACGGATAACGCTGAGATTTGCGCACTTTAACCGGCAAACCCGAATCCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGACGACCCCGCTTATCCACGTCCTGACCAGGTCTCCTCGGTGGATAACCGTAGCCATACAGATAAAGAACGGt                                          <  1:541905/248‑1 (MQ=255)
                                                                                                                                                                                                       cgaaaaCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACAAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCTCACTTTCACCGGAAAAACTGAAACCGCCTTTGAGCGCCAGTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCATCGTCCTTACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGt                                          <  1:213040/247‑1 (MQ=255)
                                                                                                                                                                                                              gCCAATAACCCCAAGCGACATCAATTCAGAAAGCGACAAATCGAACATATGGCCCACGGACACCGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCg                                   <  2:171333/251‑1 (MQ=255)
                                                                                                                                                                                                               ccAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCgg                                  <  1:184979/251‑1 (MQ=255)
                                                                                                                                                                                                                  ataatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGc                                 >  1:316926/1‑249 (MQ=255)
                                                                                                                                                                                                                   taatCCCAAGAGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGATGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGCTCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACTTCCTAACCAGGTCTCCGCGGTGGATACCCGTCGCCATACAGATAAAGCACGGTGCTACCGGAAgg                              >  2:450563/1‑251 (MQ=255)
                                                                                                                                                                                                                     atCCCAAGCGACATCAATTCAGACATCGACCAATCGAACATATTGCCCACGGACAACGCTGCGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGCTCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGtg                            >  2:321553/1‑250 (MQ=255)
                                                                                                                                                                                                                     atCCCAAGCGACATCAATTCAGAACTCGCCCAATAGAACATATTGCCCACGGCCAACGCTGAGCTTTGCGCACTTTCACCGGCAAAACTGAAACAGCCTTTGAGCGCCATTTCAGTGAAGTACGGCGCGAACTGGGTGGCGATCAGCAGGCCGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACACGGTCTCCGCGGTGGATAACCGTAGCCAAACAGAAAAAGAACGGTGCTACCGGCAGGtt                            >  2:349229/1‑250 (MQ=255)
                                                                                                                                                                                                                       cccAAGCTACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACATCGCTGAGATTTGCGCACGTGCCCCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTg                          <  2:538192/251‑1 (MQ=255)
                                                                                                                                                                                                                         cAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGGCTGCg                        >  1:152735/1‑251 (MQ=255)
                                                                                                                                                                                                                           aGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCACCGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGGCTGCg                        >  1:275259/1‑249 (MQ=255)
                                                                                                                                                                                                                             cGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCCAggg                    >  1:206218/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGACCAACGCTGTGATTTGCGCACTTTCACCGGAAAAACTGAAACCGCCTTTGAGCGCCAGTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGt                  <  2:237761/251‑1 (MQ=255)
                                                                                                                                                                                                                                     atccAGAAACCGCCCAACCGAACATATTTCCCACGGACAACGCTGAGATTTGACCACTGCCACCGGCAAAAGTGAACCCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGTATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGTGGGGTAACACt           <  2:321452/248‑1 (MQ=255)
                                                                                                                                                                                                                                         aGAAATCGACCAATCGAACATATTGCCCACGGACAAAGCTGAGATTTGCGCACTTTCACCGGCAAAACTGACACCGCCTTTGCGCGCCATTTCAGTGAAGTAAGGCGCGCAATGGGGGGCGATCAGCAGGACGGGAATCATCACGCTGACCCCGCTTAGCAACGTCCTGACCAGGTATCCGCGGTGGATAACCGTAGCCATACAGATAAAGAGCGGTGCTACCGGCAGGTATGCGGGGGGAACAACTttgt       >  1:63675/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aaaaccGACCAATCGAAAATAGTCCCCACGGACAACTCTGAGATTTGCGCGTTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttgtt      <  2:473904/246‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaTCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttgtt      <  1:223523/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaTCGCCCAATCGAACATATTGCCCACGGACAACTCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTttgttg     <  2:519268/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGACTGCCGGGGGTAACCCTTTGGTGccc  >  1:139512/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
TTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCC  >  NC_000913/4528796‑4529281

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