Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 700,928 C→A G223* (GGA→TGA)  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_000913700,9280CA98.6% 263.9 / ‑3.4 69G223* (GGA→TGA) nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (42/26);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCC  >  NC_000913/700782‑701068
                                                                                                                                                  |                                                                                                                                            
gTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCagcg                                                                                                                                           >  2:220012/1‑150 (MQ=255)
gTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGa                                                                                                                                          <  2:589192/151‑1 (MQ=255)
  tGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATg                                                                                                                                        >  1:103584/1‑151 (MQ=255)
  tGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATg                                                                                                                                        >  2:896173/1‑151 (MQ=255)
    gCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGa                                                                                                                                      <  1:1030232/151‑1 (MQ=255)
      aGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGac                                                                                                                                     >  1:68291/1‑150 (MQ=255)
         gATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGc                                                                                                                                 >  1:181657/1‑151 (MQ=255)
         gATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGc                                                                                                                                 >  2:346812/1‑151 (MQ=255)
           tAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCAc                                                                                                                               >  2:800520/1‑151 (MQ=255)
           tAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCAc                                                                                                                               >  1:630683/1‑151 (MQ=255)
                tCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGa                                                                                                                           <  2:250353/150‑1 (MQ=255)
                tCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGa                                                                                                                           <  2:751370/150‑1 (MQ=255)
                     gCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATgg                                                                                                                       >  2:655571/1‑149 (MQ=255)
                      cAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGt                                                                                                                    >  2:766425/1‑151 (MQ=255)
                       aGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTc                                                                                                                   >  1:447372/1‑151 (MQ=255)
                       aGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTc                                                                                                                   >  1:852861/1‑151 (MQ=255)
                        gCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCt                                                                                                                  <  1:242901/151‑1 (MQ=255)
                            aaGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGc                                                                                                              <  1:214953/151‑1 (MQ=255)
                            aaGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGc                                                                                                              <  1:1009340/151‑1 (MQ=255)
                                 gCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATcc                                                                                                         >  2:882522/1‑151 (MQ=255)
                                 gCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATc                                                                                                          <  2:934371/150‑1 (MQ=255)
                                          gCAGCGTTCACCCAGCGGTTCGTCCTGAATATGGCCAATTTCCCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTg                                                                                                  <  1:711141/149‑1 (MQ=255)
                                               gTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGa                                                                                            >  2:114187/1‑150 (MQ=255)
                                                ttCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCg                                                                                             >  1:726431/1‑148 (MQ=255)
                                                ttCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGa                                                                                            >  2:508823/1‑149 (MQ=255)
                                                  cACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                        <  1:246324/151‑1 (MQ=255)
                                                    cccAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGCTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCGCGGCCCCGGTTCAGCGATGGAGACGCACCTGAATGGAGTCTTCGCAATCCTGACTTGCCCCGAag                                                                                          >  1:414778/1‑147 (MQ=255)
                                                       aGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCt                                                                                   <  1:782964/151‑1 (MQ=255)
                                                        gCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTc                                                                                  <  1:745966/151‑1 (MQ=255)
                                                              cGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGCTCAGCGATGGACACGCACCAGAATGGAGTCTTCGAAATCCTGACTTGCACCGAAGTAGTGCTCcg                                                                                >  2:779974/1‑147 (MQ=255)
                                                                 ccTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgc                                                                         >  1:83054/1‑151 (MQ=255)
                                                                    gAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgcc                                                                        >  2:877282/1‑149 (MQ=255)
                                                                    gAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgcc                                                                        <  2:977047/149‑1 (MQ=255)
                                                                    gAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgcc                                                                        <  2:833637/149‑1 (MQ=255)
                                                                      atatGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGAc                                                                    >  1:340969/1‑151 (MQ=255)
                                                                               aTTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCTCCCGGTAGAGATAATCCCGCCCCCGGTTCAGCGATGGACACGCACCAGAACGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGat                                                             >  1:1004738/1‑149 (MQ=255)
                                                                                 ttCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATc                                                          <  2:388961/150‑1 (MQ=255)
                                                                                  tCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCg                                                         <  2:361197/150‑1 (MQ=255)
                                                                                    aCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGa                                                      <  1:844236/151‑1 (MQ=255)
                                                                                     ccGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAc                                                     >  1:860806/1‑151 (MQ=255)
                                                                                       gACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAcc                                                    >  1:976892/1‑150 (MQ=255)
                                                                                        aCGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCa                                                   >  2:541893/1‑150 (MQ=255)
                                                                                           ttgccgttgcGGCCAATAAAAATGCGTCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAcc                                                    <  2:682425/146‑1 (MQ=255)
                                                                                            tgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGa                                               >  2:247909/1‑150 (MQ=255)
                                                                                              ccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGaaa                                             >  1:696259/1‑150 (MQ=255)
                                                                                              ccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAAc                                            >  2:996520/1‑151 (MQ=255)
                                                                                                gttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAAc                                            >  2:644517/1‑149 (MQ=255)
                                                                                                ggtgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTGCAGCGATGGACACCCACCAGAATGGAGTCTTCGCAAGCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAg                                          <  1:1010829/149‑1 (MQ=255)
                                                                                                    cGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTc                                       >  2:279032/1‑150 (MQ=255)
                                                                                                     ggCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCAc                                     >  2:515100/1‑151 (MQ=255)
                                                                                                        cAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACttt                                  <  2:402479/151‑1 (MQ=255)
                                                                                                        cAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACttt                                  >  2:938575/1‑151 (MQ=255)
                                                                                                          aTAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTa                                 <  1:33515/150‑1 (MQ=255)
                                                                                                             aaaaTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAAc                             >  1:929397/1‑151 (MQ=255)
                                                                                                              aaaTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACg                            >  2:537488/1‑151 (MQ=255)
                                                                                                                   cgcCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttctt                       <  1:468766/151‑1 (MQ=255)
                                                                                                                    gcCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttcttc                      >  1:1060484/1‑151 (MQ=255)
                                                                                                                      ccGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCa                    <  2:572923/151‑1 (MQ=255)
                                                                                                                         ttAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCag                   >  2:825959/1‑149 (MQ=255)
                                                                                                                         ttAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCag                   >  2:718583/1‑149 (MQ=255)
                                                                                                                          tAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGc                >  1:532944/1‑151 (MQ=255)
                                                                                                                              gaTAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCt            <  2:202859/151‑1 (MQ=255)
                                                                                                                               aTAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTa           <  1:360678/151‑1 (MQ=255)
                                                                                                                                tAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAc          >  1:499190/1‑151 (MQ=255)
                                                                                                                                 aaTCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAc          <  2:668643/150‑1 (MQ=255)
                                                                                                                                  aTCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAcc         >  2:816921/1‑150 (MQ=255)
                                                                                                                                   tCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAcc         >  1:5833/1‑149 (MQ=255)
                                                                                                                                     ccGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGcc     >  1:336977/1‑151 (MQ=255)
                                                                                                                                     ccGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGcc     >  2:977544/1‑151 (MQ=255)
                                                                                                                                          cccGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGccccc  >  1:826814/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCC  >  NC_000913/700782‑701068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: