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

TTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTA  >  NC_000913/4528796‑4529285
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
ttggcGATAACCCAGTTCATCCCCAACTGCGCATCAAAACGCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGTTTTTTCCGTCCCAATACCCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCGCGTTTACGTAAGAACAGCACGACGCTGGCGACAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAAc                                                                                                                                                                                                                                                   <  1:720528/245‑1 (MQ=255)
 tGTCGATAACCCATATCATCCCCACCTCGGCATAAAAACGCGGGTCTCCGGGCATTTCCCGTATATGGCCCATGG‑CGATGCCCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTTATTGAAAGTAC‑TCCTGACACC‑GATAATCCTCTTTTTCGTAAGAACATCACGATTCGTGCGACAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGCCATCAATTCATAAATCGACCAATCGAACat                                                                                                                                                                                                                                                 <  1:680534/250‑1 (MQ=255)
  gTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCAGGTCGCAGGGCAATGCACGAAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATCATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACatat                                                                                                                                                                                                                                               >  2:731296/1‑251 (MQ=255)
    cGATAACCCAGTTCCTCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTg                                                                                                                                                                                                                                             >  1:101628/1‑251 (MQ=255)
        aaCCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCCTTGACAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATAGACCAATCGAACATATTGccc                                                                                                                                                                                                                                          >  1:364249/1‑250 (MQ=255)
           ccAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGACCAGCACGACGCTGGCGACAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGACATCACCCAATCGAACCTATTGCCCAccg                                                                                                                                                                                                                                      >  2:450237/1‑249 (MQ=255)
             aGTTCATCCCAAACGGCGCATCAAACCTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTCC‑TCCTTATTCC‑GAAAATCCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATAGTGCCCACGGAc                                                                                                                                                                                                                                    <  1:543868/251‑1 (MQ=255)
              gTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGAc                                                                                                                                                                                                                                    >  1:448148/1‑250 (MQ=255)
                    ccccAACTGCGCATCAAACCTCGGGCCGCCGGGCAATGCCCGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATACCCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATCCC‑GATAATTCACGTTTACGTAAGAACAGCCCGACGCGGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTg                                                                                                                                                                                                                             <  1:399429/251‑1 (MQ=255)
                      ccAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTga                                                                                                                                                                                                                            >  2:399845/1‑250 (MQ=255)
                       cAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTgcg                                                                                                                                                                                                                           >  2:405598/1‑248 (MQ=255)
                               cagcAAACCTCGGGTCGCGGGGCATTACCGGTAGATGGCAC‑TGGTCCATGCCCGCTTTTCCCTTCCCAATAACCTTTGT‑CATCCACCGCATAGAAAGTAC‑TCCT‑ATTCCCTAAAATTCACTTTTACGAAAGAACACAACAACTCTGGCTACAATC‑CCGGCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCATAACTCGACCAATCGAACATATTGCCCACGGACTACGCTGAGATTTGCGCa                                                                                                                                                                                                                  <  2:4675/248‑1 (MQ=255)
                                    aaccccGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGCCCGCTTTTCCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATCCC‑GATAATTCCCGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTc                                                                                                                                                                                                             <  1:310289/246‑1 (MQ=255)
                                     acaccGGGGCGCCGCGCAATGCACGTAGATGACACATGG‑CAATGACCGCTGTCTCCGTCCCAATAACCGCTGT‑CGGCCAACTCATTGAAAGTAC‑TCCTTAGCCC‑TATAATTCGCGTTTCCGTAAGAACAGCGCGACTCGGGCGACAATA‑CCTACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                            <  1:67958/247‑1 (MQ=255)
                                     aaCTCGGGTCGCCGGGCAATGCACGTAGAGGGCACAGGG‑CGATGACCGCTTTGTCCGGCCCAATACCCGTTTT‑CAGCCAACTAATTAAAAGTCC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCGGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCa                                                                                                                                                                                                            <  1:472591/251‑1 (MQ=255)
                                       cTCGGGTCGCCGGGCACTGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCACGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCAcc                                                                                                                                                                                                          >  2:688445/1‑251 (MQ=255)
                                             ccGCCGGGCAATGCACGTAGAGGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTCC‑TCCTTATCCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCGGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCaaa                                                                                                                                                                                                    <  2:172976/250‑1 (MQ=255)
                                                         gCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGcc                                                                                                                                                                                        >  2:409906/1‑251 (MQ=255)
                                                                   ggcccaTGG‑CAACGACCGCTTTTCCCGCCCCAA‑AACCCCTTTGCAGCCAACGCTTTGAAAGTAC‑CCCTTATCCC‑GAAAAGTCCCGTTTACGAAAAAGCAGCACGACGCTGGCAACAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAAGTCAGAAATCGACCAATCGAACAGATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGcc                                                                                                                                                                              <  1:447616/247‑1 (MQ=255)
                                                                    gCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCACTCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGcc                                                                                                                                                                              >  2:424773/1‑250 (MQ=255)
                                                                             cgattcccGCTTTCTCCCTCCCATTACCCTCTTC‑CAACCAACACATGGAATGTAC‑ACCTAATCCC‑AAAAATTCACGTTCACGCAAAAACAGCCCGAGGCGGGCGAAACTC‑CCGCCCGCCACGACAACGCCAATAACCCCAAGCGACATCAATTCAGAAATCGACCAATCGAAAATATTGCCCACGGACAACGCTGAGATTTGCGCACTTCCACCGGCAAAACTGAACCCGCCTTTGAGCGCCATATCAGTg                                                                                                                                                                     <  2:419016/245‑1 (MQ=255)
                                                                              caggACCGCTTTTTCCGTCCCAAGACCCACTTC‑CACCCAACTCATTAAAAGTGC‑TCCTTAATCC‑GATAATTCACTTTTACTTAAGAATAGCACGTCGCTGGC‑AAAATATCCGACCGCCACTACAACGCCAATAAGCCCAGGCGACATCAATTCAGAAAGCGACCAATGGAACATATTGCCGCCGGACAACGGTGAGATTTGCGGTCTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGa                                                                                                                                                                    <  1:78717/248‑1 (MQ=255)
                                                                                   ccGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTa                                                                                                                                                                >  1:98890/1‑250 (MQ=255)
                                                                                       tttttCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGgcg                                                                                                                                                           >  2:380781/1‑251 (MQ=255)
                                                                                              gTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTg                                                                                                                                                    >  2:150395/1‑251 (MQ=255)
                                                                                               tCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTgg                                                                                                                                                   >  2:252461/1‑251 (MQ=255)
                                                                                                cccAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTggg                                                                                                                                                  >  1:687489/1‑251 (MQ=255)
                                                                                                        ccTCTTC‑CAGCCAACTCATTAAAAGTAC‑TCCTTATTCC‑AATAATTCACGTTTACGTAAAAACAGCACGACGCTGGCGACAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                          <  1:601387/251‑1 (MQ=255)
                                                                                                        ccTCTTC‑CAGCCAACTCATGAAAATTCC‑CCCTTACTCC‑AATAATACCGGTTTACGAAAGAACACCCCGCGGCTGGCGAAAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAACCTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                          <  2:457437/251‑1 (MQ=255)
                                                                                                               ccacccAACTCATTGAAA‑TACCTCTTAATTCG‑GATAATTCCCTTTTACGAAAGAACAGCACGCCGCTGGGGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAgg                                                                                                                                    <  1:399845/247‑1 (MQ=255)
                                                                                                                             gAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCAc                                                                                                                      >  2:340235/1‑251 (MQ=255)
                                                                                                                                  tAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATg                                                                                                                  >  2:205221/1‑250 (MQ=255)
                                                                                                                                   aC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGAc                                                                                                                >  1:308878/1‑251 (MQ=255)
                                                                                                                                       ccTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATTAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGAcccc                                                                                                             >  1:231840/1‑251 (MQ=255)
                                                                                                                                          tATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGAAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATAACGATGACCCCGCt                                                                                                          >  1:712382/1‑251 (MQ=255)
                                                                                                                                          tATTCC‑GAAAATTCGCGTTTACGAAAGAACGGCCCGCCGCGGGCGACAACC‑CCGCCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTGGCGCACTTTCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCGTCACGATGACCCCGc                                                                                                           <  2:687489/250‑1 (MQ=255)
                                                                                                                                                         aCGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCTGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTg                                                                                             >  1:276656/1‑250 (MQ=255)
                                                                                                                                                         aCGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGACGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTg                                                                                             >  2:262998/1‑250 (MQ=255)
                                                                                                                                                            tttACGTAAGACCAGCCCGACGCTGGCGACAATA‑CCGACCGCCAAGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGCGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAcc                                                                                          >  2:725754/1‑250 (MQ=255)
                                                                                                                                                             ttACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGAcc                                                                                          >  1:275747/1‑249 (MQ=255)
                                                                                                                                                               aCGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGt                                                                                      >  1:49187/1‑251 (MQ=255)
                                                                                                                                                                    aGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCcg                                                                                 >  1:73676/1‑251 (MQ=255)
                                                                                                                                                                      aaaaGCACGACGCGGGCGACAATC‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCcgcg                                                                               <  1:28955/248‑1 (MQ=255)
                                                                                                                                                                        cAGCACGACGCGGGCCAAAATC‑CCGCCCGCCACGACAACGCCAATAATCCAAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGt                                                                             <  1:673746/251‑1 (MQ=255)
                                                                                                                                                                           caggacgCTGCCGACAATC‑CCGACCGCCACGACAGCGCCAATAATCCCAAGCGACATCAATACAGAAATCAACCAATCGAACATATTGCCCACGGAAAACGCTGAGATTTGCGCACTTGCACCGGCAAAACGGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGTTGGCGATCAGCAGGACGGGAATCAGCACGATGACCCCGCTTATGAACGGCTTGACCAGGTCTCCGCTGTGGa                                                                          <  2:364249/248‑1 (MQ=255)
                                                                                                                                                                              gacgCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATa                                                                        >  1:11408/1‑250 (MQ=255)
                                                                                                                                                                                cgcgggCGAAAATC‑CCGCCCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAcc                                                                     <  2:98890/247‑1 (MQ=255)
                                                                                                                                                                                cgCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCACGCGCCATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACACCGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAcc                                                                     >  2:56147/1‑251 (MQ=255)
                                                                                                                                                                                 gCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCg                                                                    >  1:230471/1‑251 (MQ=255)
                                                                                                                                                                                  gctggcaaaaatcccaaccgCCACGACAACGCCAATAATCCCAAGCGCCATCAAGTAAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGTCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTACGCGGTGGATAACCg                                                                    <  1:192567/235‑1 (MQ=255)
                                                                                                                                                                                       gtcAATA‑CCGACCGCCACGACAACGCCAATAATCCCTAGCGACATCAATTCAGAAATCGACCAATCGCACATATTGCCCACGGACAACGCTGCGATTTGCTCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGCGGCGATCAGCAGGACGGTATTCATCCCGATGACCCCGCTTATCAACGTCCTGACCAGGACTCCGCGGTGGATAACCGTAGcc                                                               >  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)
                                                                                                                                                                                                                                                   aCCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGACTGAGGGGGGGAACACTTTTTTGCCCggg   >  2:612301/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    ccAACCGAGCATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAACCTGAAACCGCCTTTGAGCGCCATTGCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGTTAATCATCACGATGACCCCACTTATCAACGTCCTGACCAGGTCTCCGCTTTTGATAACCGTAGCCATACAGATAAATATCGTTGCTTCCGGCAGGTCTGCGAGGGGTAACTCTTTGTTGCCCGGTa  <  1:354462/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
TTGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGG‑CGATGACCGCTTTTTCCGTCCCAATAACCTCTTT‑CAGCCAACTCATTGAAAGTAC‑TCCTTATTCC‑GATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATA‑CCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGTA  >  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: ‑