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

GACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATA  >  NC_000913/701309‑701595
                                                                                                                                                    |                                                                                                                                            
gacgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCggg                                                                                                                                          >  1:817305/1‑151 (MQ=255)
  cgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCggggg                                                                                                                                        <  1:869782/151‑1 (MQ=255)
               cGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCttt                                                                                                                            <  2:483020/151‑1 (MQ=255)
               cGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCttt                                                                                                                            >  2:472747/1‑151 (MQ=255)
                    aTTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCg                                                                                                                       <  2:121660/151‑1 (MQ=255)
                           ggAAATTGCGGGTGTCGGTGGCGATGGGGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACAATGGCGAGGGGAAAGATGGCTTTGCTCGGCAGTc                                                                                                                 <  2:764019/150‑1 (MQ=255)
                               aTTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAAt                                                                                                            >  1:233199/1‑151 (MQ=255)
                                 tGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATc                                                                                                           <  1:363778/150‑1 (MQ=255)
                                      gTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTGTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCga                                                                                                      <  1:109180/150‑1 (MQ=255)
                                       tttCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcgcg                                                                                                       >  1:679302/1‑148 (MQ=255)
                                       tttCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcgcg                                                                                                       >  1:314540/1‑148 (MQ=255)
                                          cGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTg                                                                                                 <  1:579268/151‑1 (MQ=255)
                                               aCGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCg                                                                                            >  2:141933/1‑151 (MQ=255)
                                                 gATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTa                                                                                          <  2:775230/151‑1 (MQ=255)
                                                    ggAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACt                                                                                        <  2:350848/150‑1 (MQ=255)
                                                             gcggcgGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAg                                                                              <  2:691364/151‑1 (MQ=255)
                                                              cggcggCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAgg                                                                             <  1:128901/151‑1 (MQ=255)
                                                               ggcggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTTGTCAATCAGGc                                                                            >  1:893442/1‑151 (MQ=255)
                                                                 cggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGa                                                                          <  2:776137/151‑1 (MQ=255)
                                                                  ggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGAt                                                                         <  2:893371/151‑1 (MQ=255)
                                                                   gcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATa                                                                        <  1:64326/151‑1 (MQ=255)
                                                                     ccccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATaa                                                                       >  1:717256/1‑150 (MQ=255)
                                                                       ccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAAcc                                                                    <  2:794883/151‑1 (MQ=255)
                                                                           tGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAAccgcc                                                                 >  2:149342/1‑150 (MQ=255)
                                                                              aGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCcgcg                                                              >  1:232572/1‑150 (MQ=255)
                                                                                gCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCt                                                            >  1:202050/1‑150 (MQ=255)
                                                                                 ccTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGt                                                          >  1:584285/1‑151 (MQ=255)
                                                                                   tgctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGtt                                                         >  2:596446/1‑150 (MQ=255)
                                                                                    gctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGtt                                                         <  2:27635/149‑1 (MQ=255)
                                                                                    gctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCGGGCGGGGGGAAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTaa                                                       <  1:417559/151‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAg                                                      <  1:603818/149‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGc                                                     <  2:305676/150‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGc                                                     <  1:554741/150‑1 (MQ=255)
                                                                                       gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCt                                                    <  1:50026/151‑1 (MQ=255)
                                                                                        aTCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTg                                                   <  2:699018/151‑1 (MQ=255)
                                                                                          cAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGtt                                                 <  1:390591/151‑1 (MQ=255)
                                                                                          cAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGt                                                  <  2:782900/150‑1 (MQ=255)
                                                                                               tcttTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTa                                              <  2:442905/149‑1 (MQ=255)
                                                                                               tcttTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTAc                                             <  1:258338/150‑1 (MQ=255)
                                                                                                    gATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGAt                                        <  1:862627/150‑1 (MQ=255)
                                                                                                     aTCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCa                                      <  1:796976/151‑1 (MQ=255)
                                                                                                      tCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCa                                      >  1:848653/1‑150 (MQ=255)
                                                                                                        aGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACa                                   <  1:290606/151‑1 (MQ=255)
                                                                                                              cgTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCACGCGATAAACCGCCCCGCTGTTAAGCTGTTTTACGAGATCAACAtt                                 >  1:66082/1‑147 (MQ=255)
                                                                                                               gTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTa                            <  2:703340/151‑1 (MQ=255)
                                                                                                                 tCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTAt                           >  1:720711/1‑150 (MQ=255)
                                                                                                                 tCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTAt                           >  2:731061/1‑150 (MQ=255)
                                                                                                                  cGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATc                          <  2:717207/150‑1 (MQ=255)
                                                                                                                        gCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTGCTGGTCAATCAGGCGGTAAACCGCCGCGCTGTTAAGCTGTTTTACGGGATCAACATTACCTATCTGAGCtt                   >  2:70145/1‑151 (MQ=255)
                                                                                                                              gcgTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTcc               <  2:349382/149‑1 (MQ=255)
                                                                                                                               cgTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCATACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGcc            >  1:248160/1‑151 (MQ=255)
                                                                                                                                gTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCt           >  2:362296/1‑151 (MQ=255)
                                                                                                                                 tAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATTTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAAGCGATAAACCGCCGCGCTGTTAAGCTGTTTTTCGAGATCAACATTACCTATCTGAGCTTGTCcccc            >  2:569132/1‑149 (MQ=255)
                                                                                                                                  aaTTTTGGTTCCGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTgg         <  2:420810/151‑1 (MQ=255)
                                                                                                                                   aTTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGGACTGGGCAATCAGGCGAGAAACCGCCCCGCTGGTTAGCTTTTTTACGCGATCAAAATTACCTCTCTGAGCTTGGCCGCCTgg         >  2:270767/1‑150 (MQ=255)
                                                                                                                                    ttttGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGgtg       >  1:146070/1‑151 (MQ=255)
                                                                                                                                    ttttGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGgt        >  2:587608/1‑150 (MQ=255)
                                                                                                                                       tGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCGATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTc     <  2:198207/150‑1 (MQ=255)
                                                                                                                                       tGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTc     >  1:890441/1‑150 (MQ=255)
                                                                                                                                        ggTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCAt   <  1:197305/151‑1 (MQ=255)
                                                                                                                                          ttACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATa  >  2:476579/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATA  >  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: