Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,252 C→A intergenic (‑9/+3) sgcC ← / ← sgcB putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,2520CA97.3% 107.2 / NA 37intergenic (‑9/+3)sgcC/sgcBputative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (11/25);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAG  >  NC_000913/4530006‑4530493
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
cTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAtttt                                                                                                                                                                                                                                               >  1:56722/1‑251 (MQ=255)
    gcacacggcAAAGACAGGCGGGGATAGGCCGAAACGCTAAAACATCCCTTTGATAGGGGGACTGAGGCTGGAGATGGCCAGGACGTACACAAGTCCCATACCGCCAAAACCATTCCCAACCGTTCCCACTGCCTGAATCGCCGGCAGAAAGGGGATGCGGAAGATAAGG‑CAAATTACGATAATGATGAAAGGGACGAAAAC‑TGTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTg                                                                                                                                                                                                                                           <  2:188772/243‑1 (MQ=255)
                      cAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGAt                                                                                                                                                                                                                         <  1:193150/251‑1 (MQ=255)
                      cAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGACGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGAt                                                                                                                                                                                                                         >  1:243042/1‑251 (MQ=255)
                        gggttAGGCCGAAACGCTCAATCATCACTTTGAAGGGGGGACTGAGGCTGCCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCGGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGAAAGGGACGAAAAC‑GGTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGAttt                                                                                                                                                                                                                       <  2:176800/247‑1 (MQ=255)
                              ggCCAAAACCCTCAATCACCCCTTGGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCATCCCCATCCCGACAAAACCAACCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGTGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACAAAAAC‑GTTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTgttgtt                                                                                                                                                                                                                 <  1:229082/251‑1 (MQ=255)
                                   aaaCGCACAATCACCACTTTAATAGGGGGACAGAGGCGTTCGATCGCCATGACGTTCACCAGTCCCATACCGACAAAACCAATCCCAACCGTGACCCCTGCCTTAATCGCCTGCAGCCGGGGGATGCGGAAGATAAGG‑CCAATCACGAGCATTATGATAGGGACGAAAAC‑GGTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTTTTGTTTTAAc                                                                                                                                                                                                            <  2:181089/251‑1 (MQ=255)
                                            accaccccttgtaaaGGGGGACGGAGGCTGTCGACCGCCATGACGACCACCAGCCCCATCCCGACAAAACAAATCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGAGGCGGAAGATAAGG‑CAAATCACGATCATGAGGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGAGGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTTTTGTTTTAACGcgtcgtcg                                                                                                                                                                                                   <  2:161041/237‑1 (MQ=255)
                                                       gATAGGGGAACTGAGGTTGCCGATCGCCAGGACGACCACCATCCCCAACCCGACAAAACCAATCCCAACCGTTCCCCCCGCCTTAATCGCCTGCAGCCAGGGGATCCGGAAGATAAGG‑CCAATCCCAATCATGATGATAGGGA‑GAAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTAGTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCtg                                                                                                                                                                                         <  1:156291/250‑1 (MQ=255)
                                                                   gAGGCTGTCGATCCCCATGACGATCCCCAGTCCCATACCGACAAAACCAATCCCACCCGTTCCCCCTGCTTAAATCGCCTGCAGCCAGGTGAGGCGGAATATAAGG‑CCAATCTCGATCATGATGATAGGGACGAAAAC‑GTTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                            <  1:247408/251‑1 (MQ=255)
                                                                     ggCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGgcgc                                                                                                                                                                          >  1:268074/1‑251 (MQ=255)
                                                                         gTCGATCGCCATGACGATCACCAGTCCCATACCGACACAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATACGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGTTCCCTGTGAGCAGGGCGGCGCCAtt                                                                                                                                                                      >  2:90622/1‑251 (MQ=255)
                                                                          tCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCACCCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTa                                                                                                                                                                                                              >  1:259028/1‑210 (MQ=255)
                                                                                   aggACGATCACCAGTCCCATACCGACAAACCCAACCCCAACCTTTCCCCCTGCTTAAACCGCCTGCAGCCAGGGGATGCGGAAGATAAGC‑CCAATCACGATAATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAggg                                                                                                                                                            <  2:56722/249‑1 (MQ=255)
                                                                                              aaGTCCCATACCGACAAACCAAAACCCACCCGTTCCCGCGGCCTTAACCGTCTCCGGCCAGGGTA‑GAGGAAAGTAAGGGCCATTCGCGATCAGGAGGATAGGGACGAAAGC‑GTTACCGCCGAGAGACAGGATGTAATCAAACATAAGTTCTCAGTTTATTGATGTAACAGAGCCGTGATGTGTTGGGTTAACGGGTCGGCGGTGATACCTGGGAGCAGGGCGGGGCCATTAAGCGTAGGGAGGCCGTAATc                                                                                                                                                 <  2:155223/250‑1 (MQ=255)
                                                                                                           acaaaacaaaGCCCAACAGTTCCCCCTCCCTAAATAGCCTGCAGCCAGGGGAGGCGGAAGATAAGG‑CCAATCACGATAATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                    <  2:62870/251‑1 (MQ=255)
                                                                                                                  caaccccAACCTTTCCCCCGGCCTTAATCGCCTGCAGCCAGGGGAGGGGGAAGATAAGG‑CCAATCACGAACATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAgg                                                                                                                             <  1:262504/247‑1 (MQ=255)
                                                                                                                     ccccAACCGTTACCCCTGCCTAAATCGCCGGCAGCCAGGGGATGCGGAAGAAAAGG‑CCAATCACGATCAGGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTa                                                                                                                          <  2:263264/250‑1 (MQ=255)
                                                                                                                           ccGTTACCCCCGCCTTAATCGCCGGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATGTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGGTACAATaa                                                                                                                    <  2:230832/251‑1 (MQ=255)
                                                                                                                                    cggcattAATCGCCTGCAGCGAGGGGGGGCGGAAGAGAAGG‑CCAATCACGAGAATGATGATAGGGCCGAAACC‑GGGACCGCCGAGAGACCGGATGTAATCAAACATAATGCCTCATTTTATTGAGTGAACAGAGCCTTGATTTGTTGTGTTAACGCGTCGGCGTTGATCCCTGTTAGCAGGGCGGCGCCATTAAGCGTAGGGGTGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATAc                                                                                                           <  1:27201/246‑1 (MQ=255)
                                                                                                                                                   gCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGCTGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCAGGGAGGTTACAATAAGGTCCATACCGTTACAGTTAag                                                                                              >  2:156482/1‑249 (MQ=255)
                                                                                                                                                   gCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAag                                                                                              <  1:190599/249‑1 (MQ=255)
                                                                                                                                                      gCCAGGGGATGCGGAAGATAAGG‑CCAATCACTATCATGATGAAAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAATCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGa                                                                                         <  2:243042/251‑1 (MQ=255)
                                                                                                                                                                   gAAGATAAGG‑CAAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACa                                                                            <  2:128190/251‑1 (MQ=255)
                                                                                                                                                                    aaGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAc                                                                             >  1:258401/1‑249 (MQ=255)
                                                                                                                                                                    aaGATAAGG‑CAAATCCCGATCAGGATGATAGGGACGAAAAC‑GGCACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTGTAGTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAAGGATGCCGTAAGCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAg                                                                           <  1:181042/251‑1 (MQ=255)
                                                                                                                                                                                aaccccGATAATGATGTAAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGGTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTATGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAATGTCCATACCGTTACAGTTAAGAGCGATCTCATTCAGACAGCATTGGGCGGt                                                                <  2:103801/246‑1 (MQ=255)
                                                                                                                                                                                        tCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTGCAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGGAGTGGCTg                                                        >  2:217686/1‑251 (MQ=255)
                                                                                                                                                                                            gatgatAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAAt                                                    >  1:206761/1‑251 (MQ=255)
                                                                                                                                                                                                 tAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCt                                               >  1:23575/1‑251 (MQ=255)
                                                                                                                                                                                                  aagggataaaactgtgactgccgAGAGACAGGATGTAATCAAACATATTTCCTCATTTTATTGAGGTAACAGAGCCTTGATTTGGTGTTGTATCGCGGCGTCGTTGCTCCCTGTGAGCAGGGCGGCGCTATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGTAGGTTACAATAATGTCCATACCGTTGCAGTTAAGAGGGATCTGATTCAAGCAGCACTGGTAGGTAGTGGCGTAAATGcgcg                                               <  2:258825/232‑4 (MQ=255)
                                                                                                                                                                                                        gcaaaC‑GGTACCGCCGAGAGACAGGATGTAATCCAACCTAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagt                                        >  2:200237/3‑251 (MQ=255)
                                                                                                                                                                                                           aaC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCGTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCTATTAAGCGTAGGGATGCCGGAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTTCATACCTTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtcag                                     <  1:129205/251‑1 (MQ=255)
                                                                                                                                                                                                                           agaCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCtt                      <  2:268074/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtgtg                 <  1:75717/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aTGTAATCAAACATTAAGCCTCATTCTATTGAGCTAACAGCGCCTTGATTTGTCGTTTTAACGCGTCGTCGTTGATCCATGCGACCAGGGCGGCCACCTTCAGCGTAGCGATGCCGTAATCACTGCCGGTCCGCATGGAGGTTACAATACGGTCCATACCGTTACAGTTAAGGGGGCTCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTAAGTAAGGAACTCTTGCAGCTTGTGTGc                >  2:122884/1‑250 (MQ=255)
                                                                                                                                                                                                                                                aaTGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGGCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAg  <  2:258401/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
CTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGG‑CCAATCACGATCATGATGATAGGGACGAAAAC‑GGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAG  >  NC_000913/4530006‑4530493

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