Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,445,359 G→C 8.5% A51A (GCC→GCG)  lexA‑I ← transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,445,3590GC 8.5% 182.0 / 5.6 59A51A (GCC→GCG) lexA‑Itranscriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base C (3/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TTGCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGTTTC  >  NC_002947/2445212‑2445506
                                                                                                                                                   |                                                                                                                                                   
ttGCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCggg                                                                                                                                                    <  1:157433/149‑1 (MQ=255)
  gCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGc                                                                                                                                                   <  2:80693/148‑1 (MQ=255)
     cGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAgg                                                                                                                                               >  5:190580/1‑149 (MQ=255)
      gATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGAAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAggg                                                                                                                                              >  4:159775/1‑149 (MQ=255)
            cTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCtt                                                                                                                                          >  1:282323/1‑147 (MQ=255)
            cTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGAAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTg                                                                                                                                         <  3:159775/148‑1 (MQ=255)
               ttCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAgg                                                                                                                                      <  6:190580/148‑1 (MQ=255)
                   gCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctc                                                                                                                                 >  3:250104/1‑149 (MQ=255)
                    cAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctcc                                                                                                                                >  1:248294/1‑149 (MQ=255)
                     aaGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctcct                                                                                                                               >  5:299534/1‑149 (MQ=255)
                      aGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctcctc                                                                                                                              >  3:123619/1‑149 (MQ=255)
                        aaTCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTcgg                                                                                                                            <  5:174394/149‑1 (MQ=255)
                          tCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTcggcg                                                                                                                          <  7:36673/149‑1 (MQ=255)
                               gcACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCAtt                                                                                                                     >  8:66264/1‑149 (MQ=255)
                                   cGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGc                                                                                                                 <  4:66784/149‑1 (MQ=255)
                                      cAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGAc                                                                                                              >  2:161477/1‑149 (MQ=255)
                                      cAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGAc                                                                                                              >  2:301958/1‑149 (MQ=255)
                                         cGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACtt                                                                                                            >  5:162314/1‑148 (MQ=255)
                                            cccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAg                                                                                                        <  6:299534/149‑1 (MQ=255)
                                              cGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGcc                                                                                                      >  7:235788/1‑149 (MQ=255)
                                               ggCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGaa                                                                                                         <  1:300074/145‑1 (MQ=255)
                                               ggCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGaa                                                                                                         >  2:300074/1‑145 (MQ=255)
                                                  cGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGc                                                                                                  >  4:287325/1‑149 (MQ=255)
                                                  cGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGc                                                                                                  <  4:250104/149‑1 (MQ=255)
                                                  cGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGc                                                                                                  >  8:265473/1‑149 (MQ=255)
                                                         gggCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGa                                                                                            <  2:282323/148‑1 (MQ=255)
                                                         gggCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAg                                                                                           >  8:96953/1‑149 (MQ=255)
                                                          ggCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAg                                                                                           <  8:103718/148‑1 (MQ=255)
                                                                 cGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCg                                                                                   >  4:168011/1‑149 (MQ=255)
                                                                        ttcAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCggg                                                                            >  8:40829/1‑149 (MQ=255)
                                                                             ccTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGc                                                                       <  1:301958/149‑1 (MQ=255)
                                                                             ccTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCATTGAGGTGATCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGc                                                                       >  6:115698/1‑149 (MQ=255)
                                                                              cTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCg                                                                      <  3:168011/149‑1 (MQ=255)
                                                                               ttGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCgg                                                                     >  1:246931/1‑149 (MQ=255)
                                                                                    ttCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGc                                                                >  7:69301/1‑149 (MQ=255)
                                                                                    ttCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGc                                                                >  4:12477/1‑149 (MQ=255)
                                                                                        aggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                            <  3:44823/149‑1 (MQ=255)
                                                                                        aggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                            <  2:246931/149‑1 (MQ=255)
                                                                                          gccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGt                                                          >  4:244624/1‑149 (MQ=255)
                                                                                          gccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGt                                                          >  6:254909/1‑149 (MQ=255)
                                                                                            caggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCt                                                        <  2:248294/149‑1 (MQ=255)
                                                                                             aggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCtt                                                       <  2:45567/149‑1 (MQ=255)
                                                                                              gggATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTc                                                      <  6:162314/149‑1 (MQ=255)
                                                                                                gatgcggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCa                                                    <  7:40829/149‑1 (MQ=255)
                                                                                                 atgcggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAg                                                                                                   <  1:11793/101‑1 (MQ=255)
                                                                                                 atgcggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAg                                                                                                   >  2:11793/1‑101 (MQ=255)
                                                                                                    cggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCTGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGc                                                 <  4:123619/148‑1 (MQ=255)
                                                                                                          ccGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCtt                                          <  5:254909/149‑1 (MQ=255)
                                                                                                                 aaGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGACAGCGCTTGATGAAc                                   <  7:66264/149‑1 (MQ=255)
                                                                                                                   gcgcCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAAcgc                                 <  3:12477/149‑1 (MQ=255)
                                                                                                                     gcCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTTTTCCAGGCAGCGCTTGATGAACGCga                               <  1:161477/149‑1 (MQ=255)
                                                                                                                                tttCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCt                    >  8:237314/1‑149 (MQ=255)
                                                                                                                                    ctcaaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCg                       >  2:300391/7‑142 (MQ=255)
                                                                                                                                    ctcaaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCg                       >  2:295728/7‑142 (MQ=255)
                                                                                                                                    ctcaaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCg                       <  1:300391/136‑1 (MQ=255)
                                                                                                                                    ctcaaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCg                       <  1:295728/136‑1 (MQ=255)
                                                                                                                                     tcaaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTtggcg               >  1:121352/6‑149 (MQ=255)
                                                                                                                                     aTCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTtggcg               >  8:197164/1‑149 (MQ=255)
                                                                                                                                         cgcCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTtggcgtggc           >  1:295903/1‑149 (MQ=255)
                                                                                                                                            ccTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTc        <  4:10640/149‑1 (MQ=255)
                                                                                                                                                  gggCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGttt   >  6:264501/1‑148 (MQ=255)
                                                                                                                                                  gggCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGTTTc  <  7:265473/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGTTTC  >  NC_002947/2445212‑2445506

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