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,0390CT98.6% 218.8 / NA 71G402D (GGT→GAT) sgcCputative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (36/33);  minor base G (0/1);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTA  >  NC_000913/4528796‑4529285
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
ttggcGATAACCCAGTTCATCCCCAACTGCGCATCAAAACGCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGTTTTTTCCGTCCCAATACCCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCGCGTTTACGTAAGAACAGCACGACGCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAAc                                                                                                                                                                                                                                                   <  1:720528/245‑1 (MQ=255)
 tGTCGATAACCCATATCATCCCCACCTCGGCATAAAAACGCGGGTCTCCGGGCATTTCCCGTATATGGCCCATGGCGATGCCCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTTATTGAAAGTACTCCTGACACCGATAATCCTCTTTTTCGTAAGAACATCACGATTCGTGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGCCATCAATTCATAAATCGACCAATCGAACat                                                                                                                                                                                                                                                 <  1:680534/250‑1 (MQ=255)
  gTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCAGGTCGCAGGGCAATGCACGAAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATCATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACatat                                                                                                                                                                                                                                               >  2:731296/1‑251 (MQ=255)
    cGATAACCCAGTTCCTCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTg                                                                                                                                                                                                                                             >  1:101628/1‑251 (MQ=255)
        aaCCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCCTTGACAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATAGACCAATCGAACATATTGccc                                                                                                                                                                                                                                          >  1:364249/1‑250 (MQ=255)
           ccAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGACCAGCACGACGCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGACATCACCCAATCGAACCTATTGCCCAccg                                                                                                                                                                                                                                      >  2:450237/1‑249 (MQ=255)
             aGTTCATCCCAAACGGCGCATCAAACCTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTCCTCCTTATTCCGAAAATCCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATAGTGCCCACGGAc                                                                                                                                                                                                                                    <  1:543868/251‑1 (MQ=255)
              gTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGAc                                                                                                                                                                                                                                    >  1:448148/1‑250 (MQ=255)
                    ccccAACTGCGCATCAAACCTCGGGCCGCCGGGCAATGCCCGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATACCCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATCCCGATAATTCACGTTTACGTAAGAACAGCCCGACGCGGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTg                                                                                                                                                                                                                             <  1:399429/251‑1 (MQ=255)
                      ccAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTga                                                                                                                                                                                                                            >  2:399845/1‑250 (MQ=255)
                       cAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTgcg                                                                                                                                                                                                                           >  2:405598/1‑248 (MQ=255)
                               cagcAAACCTCGGGTCGCGGGGCATTACCGGTAGATGGCACTGGTCCATGCCCGCTTTTCCCTTCCCAATAACCTTTGTCATCCACCGCATAGAAAGTACTCCTATTCCCTAAAATTCACTTTTACGAAAGAACACAACAACTCTGGCTACAATCCCGGCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCATAACTCGACCAATCGAACATATTGCCCACGGACTACGCTGAGATTTGCGCa                                                                                                                                                                                                                  <  2:4675/248‑1 (MQ=255)
                                    aaccccGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGCCCGCTTTTCCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATCCCGATAATTCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTc                                                                                                                                                                                                             <  1:310289/246‑1 (MQ=255)
                                     acaccGGGGCGCCGCGCAATGCACGTAGATGACACATGGCAATGACCGCTGTCTCCGTCCCAATAACCGCTGTCGGCCAACTCATTGAAAGTACTCCTTAGCCCTATAATTCGCGTTTCCGTAAGAACAGCGCGACTCGGGCGACAATACCTACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                            <  1:67958/247‑1 (MQ=255)
                                     aaCTCGGGTCGCCGGGCAATGCACGTAGAGGGCACAGGGCGATGACCGCTTTGTCCGGCCCAATACCCGTTTTCAGCCAACTAATTAAAAGTCCTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCGGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                            <  1:472591/251‑1 (MQ=255)
                                       cTCGGGTCGCCGGGCACTGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCACGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCAcc                                                                                                                                                                                                          >  2:688445/1‑251 (MQ=255)
                                             ccGCCGGGCAATGCACGTAGAGGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTCCTCCTTATCCCGATAATTCACGTTTACGTAAGAACAGCACGACGCGGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCaaa                                                                                                                                                                                                    <  2:172976/250‑1 (MQ=255)
                                                         gCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGcc                                                                                                                                                                                        >  2:409906/1‑251 (MQ=255)
                                                                   ggcccaTGGCAACGACCGCTTTTCCCGCCCCAAAACCCCTTTGCAGCCAACGCTTTGAAAGTACCCCTTATCCCGAAAAGTCCCGTTTACGAAAAAGCAGCACGACGCTGGCAACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAAGTCAGAAATCGACCAATCGAACAGATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGcc                                                                                                                                                                              <  1:447616/247‑1 (MQ=255)
                                                                    gCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCACTCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGcc                                                                                                                                                                              >  2:424773/1‑250 (MQ=255)
                                                                             cgattcccGCTTTCTCCCTCCCATTACCCTCTTCCAACCAACACATGGAATGTACACCTAATCCCAAAAATTCACGTTCACGCAAAAACAGCCCGAGGCGGGCGAAACTCCCGCCCGCCACGACAACGCCAATAACCCCAAGCGACATCAATTCAGAAATCGACCAATCGAAAATATTGCCCACGGACAACGCTGAGATTTGCGCACTTCCACCGGCAAAACTGAACCCGCCTTTGAGCGCCATATCAGTg                                                                                                                                                                     <  2:419016/245‑1 (MQ=255)
                                                                              caggACCGCTTTTTCCGTCCCAAGACCCACTTCCACCCAACTCATTAAAAGTGCTCCTTAATCCGATAATTCACTTTTACTTAAGAATAGCACGTCGCTGGCAAAATATCCGACCGCCACTACAACGCCAATAAGCCCAGGCGACATCAATTCAGAAAGCGACCAATGGAACATATTGCCGCCGGACAACGGTGAGATTTGCGGTCTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGa                                                                                                                                                                    <  1:78717/248‑1 (MQ=255)
                                                                                   ccGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTa                                                                                                                                                                >  1:98890/1‑250 (MQ=255)
                                                                                       tttttCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGgcg                                                                                                                                                           >  2:380781/1‑251 (MQ=255)
                                                                                              gTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTg                                                                                                                                                    >  2:150395/1‑251 (MQ=255)
                                                                                               tCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTgg                                                                                                                                                   >  2:252461/1‑251 (MQ=255)
                                                                                                cccAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTggg                                                                                                                                                  >  1:687489/1‑251 (MQ=255)
                                                                                                        ccTCTTCCAGCCAACTCATTAAAAGTACTCCTTATTCCAATAATTCACGTTTACGTAAAAACAGCACGACGCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                          <  1:601387/251‑1 (MQ=255)
                                                                                                        ccTCTTCCAGCCAACTCATGAAAATTCCCCCTTACTCCAATAATACCGGTTTACGAAAGAACACCCCGCGGCTGGCGAAAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAACCTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                          <  2:457437/251‑1 (MQ=255)
                                                                                                               ccacccAACTCATTGAAATACCTCTTAATTCGGATAATTCCCTTTTACGAAAGAACAGCACGCCGCTGGGGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAgg                                                                                                                                    <  1:399845/247‑1 (MQ=255)
                                                                                                                             gAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCAc                                                                                                                      >  2:340235/1‑251 (MQ=255)
                                                                                                                                  tACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATg                                                                                                                  >  2:205221/1‑250 (MQ=255)
                                                                                                                                   aCTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGAc                                                                                                                >  1:308878/1‑251 (MQ=255)
                                                                                                                                       ccTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATTAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGAcccc                                                                                                             >  1:231840/1‑251 (MQ=255)
                                                                                                                                          tATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGAAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATAACGATGACCCCGCt                                                                                                          >  1:712382/1‑251 (MQ=255)
                                                                                                                                          tATTCCGAAAATTCGCGTTTACGAAAGAACGGCCCGCCGCGGGCGACAACCCCGCCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTGGCGCACTTTCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCGTCACGATGACCCCGc                                                                                                           <  2:687489/250‑1 (MQ=255)
                                                                                                                                                         aCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCTGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTg                                                                                             >  1:276656/1‑250 (MQ=255)
                                                                                                                                                         aCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGACGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTg                                                                                             >  2:262998/1‑250 (MQ=255)
                                                                                                                                                            tttACGTAAGACCAGCCCGACGCTGGCGACAATACCGACCGCCAAGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGCGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAcc                                                                                          >  2:725754/1‑250 (MQ=255)
                                                                                                                                                             ttACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAcc                                                                                          >  1:275747/1‑249 (MQ=255)
                                                                                                                                                               aCGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGt                                                                                      >  1:49187/1‑251 (MQ=255)
                                                                                                                                                                    aGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCcg                                                                                 >  1:73676/1‑251 (MQ=255)
                                                                                                                                                                      aaaaGCACGACGCGGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCcgcg                                                                               <  1:28955/248‑1 (MQ=255)
                                                                                                                                                                        cAGCACGACGCGGGCCAAAATCCCGCCCGCCACGACAACGCCAATAATCCAAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGt                                                                             <  1:673746/251‑1 (MQ=255)
                                                                                                                                                                           caggacgCTGCCGACAATCCCGACCGCCACGACAGCGCCAATAATCCCAAGCGACATCAATACAGAAATCAACCAATCGAACATATTGCCCACGGAAAACGCTGAGATTTGCGCACTTGCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGTTGGCGATCAGCAGGACGGGAATCAGCACGATGACCCCGCTTATGAACGGCTTGACCAGGTCTCCGCTGTGGa                                                                          <  2:364249/248‑1 (MQ=255)
                                                                                                                                                                              gacgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATa                                                                        >  1:11408/1‑250 (MQ=255)
                                                                                                                                                                                cgcgggCGAAAATCCCGCCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAcc                                                                     <  2:98890/247‑1 (MQ=255)
                                                                                                                                                                                cgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCACGCGCCATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACACCGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAcc                                                                     >  2:56147/1‑251 (MQ=255)
                                                                                                                                                                                 gCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCg                                                                    >  1:230471/1‑251 (MQ=255)
                                                                                                                                                                                  gctggcaaaaatcccaaccgCCACGACAACGCCAATAATCCCAAGCGCCATCAAGTAAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGTCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTACGCGGTGGATAACCg                                                                    <  1:192567/235‑1 (MQ=255)
                                                                                                                                                                                       gtcAATACCGACCGCCACGACAACGCCAATAATCCCTAGCGACATCAATTCAGAAATCGACCAATCGCACATATTGCCCACGGACAACGCTGCGATTTGCTCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGCGGCGATCAGCAGGACGGTATTCATCCCGATGACCCCGCTTATCAACGTCCTGACCAGGACTCCGCGGTGGATAACCGTAGcc                                                               >  2:733914/3‑250 (MQ=255)
                                                                                                                                                                                          caaccccaccagCCACGACAACGCCAATAATCCCAACCGACAACAATTCAGAAAGCGCCCAAGCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCATTTTCACCGTCAAAACTGAAACCGCCTTTGAGCGCCTGTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGAGGACCCCGCTTATCAGCGTCCTGACCAGGTCGCCGCTGTGGATTACCGTAGCCATa                                                            <  2:375409/240‑1 (MQ=255)
                                                                                                                                                                                          aaacaccaaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGCCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATa                                                            <  2:135079/243‑1 (MQ=255)
                                                                                                                                                                                               ccacccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGAt                                                       <  2:475521/247‑1 (MQ=255)
                                                                                                                                                                                                 gtccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGTTGAGATTTGCGCACTTGCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGTCGCGAACTGGGTGGCGAGCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATaa                                                     <  1:584378/249‑1 (MQ=255)
                                                                                                                                                                                                 gaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGCCCCCGCTTATCAACGTCATGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAgcta                                                      >  2:417312/1‑247 (MQ=255)
                                                                                                                                                                                                  accgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATaa                                                     >  2:121146/1‑250 (MQ=255)
                                                                                                                                                                                                     gcaactaCAACGCCAATAATCCCAAGCAACATCAATTCAGAAATCGACCAATCGAACATATGGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGaa                                                 <  1:13440/245‑1 (MQ=255)
                                                                                                                                                                                                            cAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCt                                          <  1:409906/251‑1 (MQ=255)
                                                                                                                                                                                                             aaCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTa                                         <  1:380781/251‑1 (MQ=255)
                                                                                                                                                                                                                   aagaatCCCAAGCGACATCAATGAAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCa                                   <  2:73676/248‑1 (MQ=255)
                                                                                                                                                                                                                                 aCATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGt                     <  1:150395/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAac                  <  2:11408/250‑1 (MQ=255)
                                                                                                                                                                                                                                         tCAGAAATCGACCAATCGAACATATTGCCCCCGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACtt              <  2:275747/250‑1 (MQ=255)
                                                                                                                                                                                                                                           agacaTCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGGAACAATttgt           >  2:545405/5‑251 (MQ=255)
                                                                                                                                                                                                                                             aaaTCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGCAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCCTACAGATCAAGAACGTTGCTACCGGCAGGTTTGCTCGGGGTAACACTtttttg         >  2:399617/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGGAACACTTTGTTGCCCg     >  1:448115/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCgg    <  1:340235/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCgg    <  1:725754/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGACTGAGGGGGGGAACACTTTTTTGCCCgg  >  2:612301/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    ccAACCGAGCATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAACCTGAAACCGCCTTTGAGCGCCATTGCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGTTAATCATCACGATGACCCCACTTATCAACGTCCTGACCAGGTCTCCGCTTTTGATAACCGTAGCCATACAGATAAATATCGTTGCTTCCGGCAGGTCTGCGAGGGGTAACTCTTTGTTGCCCGGTa  <  1:354462/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
TTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTA  >  NC_000913/4528796‑4529285

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