Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,861 C→A 5.6% Q640K (CAA→AAA)  nagE → N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,8610CA 5.6% 276.9 / 2.1 79Q640K (CAA→AAA) nagEN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base C (41/33);  new base A (2/3);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGGCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGT  >  NC_000913/705722‑705997
                                                                                                                                           |                                                                                                                                        
cgcAGGGCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTc                                                                                                                                          >  1:398131/1‑140 (MQ=255)
   aGGGCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCaaa                                                                                                                                       <  1:455794/140‑1 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTTGCGGGTCAAaca                                                                                                                                     >  1:859479/1‑140 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAaca                                                                                                                                     >  1:1173109/1‑140 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAaca                                                                                                                                     >  1:306429/1‑140 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAaca                                                                                                                                     >  1:333722/1‑140 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAaca                                                                                                                                     >  1:333867/1‑140 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAaca                                                                                                                                     >  1:12682/1‑140 (MQ=255)
     ggCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCCGGTCAAaca                                                                                                                                     >  1:859252/1‑140 (MQ=255)
      gCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAacac                                                                                                                                    <  1:1536502/140‑1 (MQ=255)
               tCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACTACGCTGTATg                                                                                                                           >  2:1686038/1‑140 (MQ=255)
               tCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATg                                                                                                                           >  2:1686331/1‑140 (MQ=255)
               tCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATg                                                                                                                           <  2:882416/140‑1 (MQ=255)
                               gATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATc                                                                                                           <  2:1102529/140‑1 (MQ=255)
                               gATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATc                                                                                                           <  2:601885/140‑1 (MQ=255)
                                 ttACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTg                                                                                                         >  1:1401478/1‑140 (MQ=255)
                                    ccTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCtt                                                                                                      <  1:615608/140‑1 (MQ=255)
                                     cTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCttt                                                                                                     <  1:518278/140‑1 (MQ=255)
                                         aCGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGc                                                                                                 >  2:331745/1‑140 (MQ=255)
                                         aCGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCACTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGc                                                                                                 >  2:1588838/1‑140 (MQ=255)
                                              aaCGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGAt                                                                                            >  1:172637/1‑140 (MQ=255)
                                              aaCGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGAt                                                                                            >  1:1357718/1‑140 (MQ=255)
                                              aaCGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCAAATTGTGGCGGGTCAAACACCCCAGAAAGAAATCAAAAAGTAAACAGCTTTATGCCTGAt                                                                                            >  1:172629/1‑140 (MQ=255)
                                                cGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGc                                                                                          <  2:1640517/140‑1 (MQ=255)
                                                  cccGCTCGATGATTAGCCCGGTGGTTTGCAGCAAAATCGACGATTTCAGAGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGa                                                                                        >  1:543915/1‑140 (MQ=255)
                                                            tatTTGCCTGGTGGTTTGAAGCAATATCGACGATTTCAGTGGATTGATCATTAATGATCAGGGCAATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCg                                                                              <  1:1686038/139‑1 (MQ=255)
                                                            gatTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGATTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCg                                                                              <  1:1686331/140‑1 (MQ=255)
                                                              ttAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                            <  2:908021/140‑1 (MQ=255)
                                                                 gCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAACAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCa                                                                         >  2:1237332/1‑140 (MQ=255)
                                                                 gCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCa                                                                         <  1:695076/140‑1 (MQ=255)
                                                                 gCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCa                                                                         <  1:695137/140‑1 (MQ=255)
                                                                     ggtggtTTTCAGCAATATCGACGTTTTCTGTGGCTTGTTCATTAATGCTCAGGTCAATTTTGTGGCGGGTCAAACACCTCTGTTTGAATTCAAAAAGTTATCTGATTTATGCATGATGCGAGGCTTGAGCGTCgattcct                                                                     <  2:434639/140‑7 (MQ=255)
                                                                     ggtggtTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCca                                                                     <  2:641359/140‑1 (MQ=255)
                                                                     ggtggtTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCca                                                                     <  2:1011752/140‑1 (MQ=255)
                                                                        ggtTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCcaaca                                                                  <  1:834815/140‑1 (MQ=255)
                                                                            tGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGa                                                              >  2:1332645/1‑140 (MQ=255)
                                                                              cagaaATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCTGGGCCATATTGTGGAGGGTAAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACa                                                            <  2:1244103/140‑1 (MQ=255)
                                                                              aagcaATATCGACGATTTCATTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTTTGAAATCTAAAACTATTCTGCTTTATGCCTGATGCGACGAATGAGCGTCGCATCCAACAATGACa                                                            <  2:1038882/139‑1 (MQ=255)
                                                                                gcaATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAg                                                          <  2:262331/140‑1 (MQ=255)
                                                                                 caATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGc                                                         >  1:423932/1‑140 (MQ=255)
                                                                                        gacgaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGgag                                                  <  1:1185633/140‑1 (MQ=255)
                                                                                                cAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                          >  2:791016/1‑140 (MQ=255)
                                                                                                cAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                          >  2:791234/1‑140 (MQ=255)
                                                                                                       ttGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAATCGGTGGAGATCTTCTCTGCCGCt                                   >  1:806210/1‑140 (MQ=255)
                                                                                                       ttGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCt                                   >  1:47973/1‑140 (MQ=255)
                                                                                                               tAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAAACAAAATGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttttt                           >  2:1492765/1‑140 (MQ=255)
                                                                                                                      cAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCatcaat                    <  2:398131/140‑1 (MQ=255)
                                                                                                                      cAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCatcaat                    <  2:1038424/140‑1 (MQ=255)
                                                                                                                      cAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAAACAAAAAGTAATCTGCTTTTTGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTATTTTCatcaat                    <  2:132997/140‑1 (MQ=255)
                                                                                                                        gggCCATATGGTGGCGGGTCAAACACCGCTGTATTAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATATTCTCTGCCGCTTTTTTTTTCatcagtca                  <  2:231740/140‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAATAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAatcat                 <  1:607799/140‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTTGAGATCTTCTCTGCCGCTTTTTTTTTCATCAatcat                 <  2:217609/140‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAatcat                 <  2:630206/140‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAatcat                 <  1:1201393/140‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAatcat                 <  1:158385/140‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTAAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGa                                                                                        <  2:750602/69‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTAAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGa                                                                                        <  2:750085/69‑1 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTAAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGa                                                                                        >  1:750602/1‑69 (MQ=255)
                                                                                                                         ggCCATATTGTGGCGGGTAAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGa                                                                                        >  1:750085/1‑69 (MQ=255)
                                                                                                                            cATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATccc              <  2:1364091/140‑1 (MQ=255)
                                                                                                                              tatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCa            >  2:136018/1‑140 (MQ=255)
                                                                                                                                  gtgGCGGGTCTATCACCGCTGTTTGATATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAAt        <  1:1653086/140‑1 (MQ=255)
                                                                                                                                  gtgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAAt        <  1:1653114/140‑1 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCTCTGAATGAAATCAAAAAGTAATCAGCTATATGCCTGATTCGACCCTTGATCTTCGCATCCAACAATGACAAACTGAGGAAATCTAATCTGCCGCTATAATTTACAGCACACATCCCCaaaaac       >  2:1668687/1‑138 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATc       >  2:695901/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATc       >  2:1564553/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATc       >  2:367811/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATc       >  2:1566874/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTATTTTCATCAATCATCCCCATAATc       >  2:1668695/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTGTTTTTTCATCAATCAACCCCATAATc       >  2:1544359/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATcc      >  2:831571/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGAAAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATc       >  1:189416/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAAAGACAAGCGGTGGAGATCTTCTCTGCCGCTATTTTTTTCATCAATCATCCCCATAATc       >  2:1544448/1‑140 (MQ=255)
                                                                                                                                   tgGCGGGTCAAACACCGCAGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCATGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTACACTGCCGCTTTTTTTATCATCAATCAACCCCAAAATc       >  2:1544522/1‑140 (MQ=255)
                                                                                                                                     gCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCt     >  1:1234751/1‑140 (MQ=255)
                                                                                                                                     gCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCt     >  1:1234759/1‑140 (MQ=255)
                                                                                                                                       gggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTg   <  2:1247354/140‑1 (MQ=255)
                                                                                                                                        ggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCgtgt  >  1:1067127/1‑140 (MQ=255)
                                                                                                                                        ggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGt  >  1:1607580/1‑140 (MQ=255)
                                                                                                                                        ggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGt  >  1:499744/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                        
CGCAGGGCAACCGATTCTGGAAATGGATCTGGATTACCTGAACGCTAACGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGT  >  NC_000913/705722‑705997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑