Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 701,461 (G)5→6 coding (134/1221 nt) nagC ← N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913701,4561.G98.3% 204.3 / ‑4.2 59coding (139/1221 nt)nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (25/33);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑G‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATA  >  NC_000913/701309‑701595
                                                                                                                                                    |                                                                                                                                            
gacgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑g‑gg                                                                                                                                          >  1:817305/1‑151 (MQ=255)
  cgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑g‑gggg                                                                                                                                        <  1:869782/151‑1 (MQ=255)
               cGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCttt                                                                                                                            <  2:483020/151‑1 (MQ=255)
               cGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCttt                                                                                                                            >  2:472747/1‑151 (MQ=255)
                    aTTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCg                                                                                                                       <  2:121660/151‑1 (MQ=255)
                           ggAAATTGCGGGTGTCGGTGGCGATGGGGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACAATGGC‑GAGGGGAAAGATGGCTTTGCTCGGCAGTc                                                                                                                 <  2:764019/150‑1 (MQ=255)
                               aTTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAAt                                                                                                            >  1:233199/1‑151 (MQ=255)
                                 tGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATc                                                                                                           <  1:363778/150‑1 (MQ=255)
                                      gTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTGTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCga                                                                                                      <  1:109180/150‑1 (MQ=255)
                                       tttCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcgcg                                                                                                       >  1:679302/1‑148 (MQ=255)
                                       tttCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcgcg                                                                                                       >  1:314540/1‑148 (MQ=255)
                                          cGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTg                                                                                                 <  1:579268/151‑1 (MQ=255)
                                               aCGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCg                                                                                            >  2:141933/1‑151 (MQ=255)
                                                 gATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTa                                                                                          <  2:775230/151‑1 (MQ=255)
                                                    ggAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACt                                                                                        <  2:350848/150‑1 (MQ=255)
                                                             gcggcgGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAg                                                                              <  2:691364/151‑1 (MQ=255)
                                                              cggcggCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAgg                                                                             <  1:128901/151‑1 (MQ=255)
                                                               ggcggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTTGTCAATCAGGc                                                                            >  1:893442/1‑151 (MQ=255)
                                                                 cggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGa                                                                          <  2:776137/151‑1 (MQ=255)
                                                                  ggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGAt                                                                         <  2:893371/151‑1 (MQ=255)
                                                                   gcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATa                                                                        <  1:64326/151‑1 (MQ=255)
                                                                     ccccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATaa                                                                       >  1:717256/1‑150 (MQ=255)
                                                                       ccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAAcc                                                                    <  2:794883/151‑1 (MQ=255)
                                                                           tGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAAccgcc                                                                 >  2:149342/1‑150 (MQ=255)
                                                                              aGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCcgcg                                                              >  1:232572/1‑150 (MQ=255)
                                                                                gCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCt                                                            >  1:202050/1‑150 (MQ=255)
                                                                                 ccTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGt                                                          >  1:584285/1‑151 (MQ=255)
                                                                                   tgctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGtt                                                         >  2:596446/1‑150 (MQ=255)
                                                                                    gctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGtt                                                         <  2:27635/149‑1 (MQ=255)
                                                                                    gctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCGGGCGG‑GGGGAAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTaa                                                       <  1:417559/151‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAg                                                      <  1:603818/149‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGc                                                     <  2:305676/150‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGc                                                     <  1:554741/150‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCt                                                    <  1:50026/151‑1 (MQ=255)
                                                                                        aTCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTg                                                   <  2:699018/151‑1 (MQ=255)
                                                                                          cAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGtt                                                 <  1:390591/151‑1 (MQ=255)
                                                                                          cAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGt                                                  <  2:782900/150‑1 (MQ=255)
                                                                                               tcttTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTa                                              <  2:442905/149‑1 (MQ=255)
                                                                                               tcttTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTAc                                             <  1:258338/150‑1 (MQ=255)
                                                                                                    gATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGAt                                        <  1:862627/150‑1 (MQ=255)
                                                                                                     aTCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCa                                      <  1:796976/151‑1 (MQ=255)
                                                                                                      tCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCa                                      >  1:848653/1‑150 (MQ=255)
                                                                                                        aGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACa                                   <  1:290606/151‑1 (MQ=255)
                                                                                                              cgTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCACGCGATAAACCGCCCCGCTGTTAAGCTGTTTTACGAGATCAACAtt                                 >  1:66082/1‑147 (MQ=255)
                                                                                                               gTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTa                            <  2:703340/151‑1 (MQ=255)
                                                                                                                 tCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTAt                           >  1:720711/1‑150 (MQ=255)
                                                                                                                 tCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTAt                           >  2:731061/1‑150 (MQ=255)
                                                                                                                  cGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATc                          <  2:717207/150‑1 (MQ=255)
                                                                                                                        gCTGACGCGTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTGCTGGTCAATCAGGCGGTAAACCGCCGCGCTGTTAAGCTGTTTTACGGGATCAACATTACCTATCTGAGCtt                   >  2:70145/1‑151 (MQ=255)
                                                                                                                              gcgTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTcc               <  2:349382/149‑1 (MQ=255)
                                                                                                                               cgTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCATACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGcc            >  1:248160/1‑151 (MQ=255)
                                                                                                                                gTAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCt           >  2:362296/1‑151 (MQ=255)
                                                                                                                                 tAATTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATTTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAAGCGATAAACCGCCGCGCTGTTAAGCTGTTTTTCGAGATCAACATTACCTATCTGAGCTTGTCcccc            >  2:569132/1‑149 (MQ=255)
                                                                                                                                  aaTTTTGGTTCCGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTgg         <  2:420810/151‑1 (MQ=255)
                                                                                                                                   aTTTTGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGGACTGGGCAATCAGGCGAGAAACCGCCCCGCTGGTTAGCTTTTTTACGCGATCAAAATTACCTCTCTGAGCTTGGCCGCCTgg         >  2:270767/1‑150 (MQ=255)
                                                                                                                                    ttttGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGgtg       >  1:146070/1‑151 (MQ=255)
                                                                                                                                    ttttGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGgt        >  2:587608/1‑150 (MQ=255)
                                                                                                                                       tGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCGATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTc     <  2:198207/150‑1 (MQ=255)
                                                                                                                                       tGGTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTc     >  1:890441/1‑150 (MQ=255)
                                                                                                                                        ggTTACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCAt   <  1:197305/151‑1 (MQ=255)
                                                                                                                                          ttACGCTGGCGG‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATa  >  2:476579/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑G‑GGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATA  >  NC_000913/701309‑701595

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