Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,445,402 G→C 6.7% 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,4020GC 6.7% 186.8 / 4.7 60A51A (GCC→GCGlexA‑Itranscriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base C (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

TTGCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGTTTCAA  >  NC_002947/2445255‑2445551
                                                                                                                                                   |                                                                                                                                                     
ttGCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGc                                                                                                                                                     >  1:871152/1‑150 (MQ=255)
  gCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCaa                                                                                                                                                   >  1:1094597/1‑150 (MQ=255)
  gCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCaa                                                                                                                                                   <  1:58199/150‑1 (MQ=255)
   cTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAg                                                                                                                                                  >  2:455113/1‑150 (MQ=255)
   cTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAg                                                                                                                                                  <  2:1171524/150‑1 (MQ=255)
    tCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAgg                                                                                                                                                 >  1:195146/1‑150 (MQ=255)
    tCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAgg                                                                                                                                                 >  1:719411/1‑150 (MQ=255)
        tgtgCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCCCGCAAGGGcgc                                                                                                                                             <  1:494368/150‑3 (MQ=255)
          tgCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTg                                                                                                                                           >  1:211789/1‑150 (MQ=255)
          tgCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCCCGCAAGGGcgcga                                                                                                                                           <  1:527145/150‑5 (MQ=255)
           gCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGa                                                                                                                                          <  2:871152/150‑1 (MQ=255)
           gCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGa                                                                                                                                          <  2:719411/150‑1 (MQ=255)
            cTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAg                                                                                                                                         >  1:106492/1‑150 (MQ=255)
                  ggCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctc                                                                                                                                   >  1:82093/1‑150 (MQ=255)
                   gCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctcc                                                                                                                                  <  2:1094597/150‑1 (MQ=255)
                     aaGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCgg                                                                                                                                                       <  1:1198765/127‑1 (MQ=255)
                     aaGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCgg                                                                                                                                                       >  2:1198765/1‑127 (MQ=255)
                     aaGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctc                                                                                                                                   <  1:778632/147‑1 (MQ=255)
                     aaGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGctc                                                                                                                                   >  2:778632/1‑147 (MQ=255)
                              cgcACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCAtt                                                                                                                       >  1:269901/1‑150 (MQ=255)
                              cgcACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCAtt                                                                                                                       <  2:269901/150‑1 (MQ=255)
                                cACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTgg                                                                                                                     <  1:1259266/150‑1 (MQ=255)
                                            cccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAg                                                                                                                                         <  1:517106/118‑1 (MQ=255)
                                            cccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAg                                                                                                                                         >  2:517106/1‑118 (MQ=255)
                                             ccGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGcc                                                                                                        >  2:1328226/1‑150 (MQ=255)
                                               ggCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGc                                                                                                         >  2:954066/1‑147 (MQ=255)
                                               ggCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGc                                                                                                         <  1:954066/147‑1 (MQ=255)
                                               ggCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCa                                                                                                      >  1:1383559/1‑150 (MQ=255)
                                                       atgGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGa                                                                                              <  1:590072/150‑1 (MQ=255)
                                                       atgGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGa                                                                                              <  2:211789/150‑1 (MQ=255)
                                                           gCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCaa                                                                                          >  1:555860/1‑150 (MQ=255)
                                                             aGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAAtc                                                                                        >  2:538994/1‑150 (MQ=255)
                                                                ccGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGc                                                                                                                                               >  2:1304753/1‑92 (MQ=255)
                                                                ccGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGc                                                                                                                                               <  1:1304753/92‑1 (MQ=255)
                                                                   gCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAAt                                                                                         <  2:911424/143‑1 (MQ=255)
                                                                   gCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAAt                                                                                         >  1:911424/1‑143 (MQ=255)
                                                                      tcttcAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                               <  1:538994/150‑1 (MQ=255)
                                                                               ttGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCggg                                                                      >  2:1174537/1‑150 (MQ=255)
                                                                                  gCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCgg                                                                                    <  1:1109827/133‑1 (MQ=255)
                                                                                  gCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCgg                                                                                    >  2:1109827/1‑133 (MQ=255)
                                                                                  gCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAg                                                                   >  2:1348541/1‑150 (MQ=255)
                                                                                    ttCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGcc                                                                 <  2:1383559/150‑1 (MQ=255)
                                                                                       caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGCGCACTTGCGGCGCGGGTAGGCGGGAAGCCGtt                                                              <  1:1328226/150‑1 (MQ=255)
                                                                                       caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                              <  1:398704/150‑1 (MQ=255)
                                                                                       caggccaggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGtt                                                              <  1:1348541/150‑1 (MQ=255)
                                                                                             aggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGaa                                                                                                           <  1:545922/99‑1 (MQ=255)
                                                                                             aggGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGaa                                                                                                           >  2:545922/1‑99 (MQ=255)
                                                                                               ggatgcggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCa                                                      <  2:143054/150‑1 (MQ=255)
                                                                                                  tgcggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCa                                                                                                      <  1:982553/99‑1 (MQ=255)
                                                                                                  tgcggatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCa                                                                                                      >  2:982553/1‑99 (MQ=255)
                                                                                                      gatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGATTTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                               <  2:903981/118‑1 (MQ=255)
                                                                                                      gatgcCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCgg                                                                               >  1:903981/1‑118 (MQ=255)
                                                                                                           cgcgGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGc                                                                                                                                               <  2:530853/49‑1 (MQ=255)
                                                                                                           cgcgGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGc                                                                                                                                               >  1:530853/1‑49 (MQ=255)
                                                                                                           cgcgGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGCAGCCCAGTTTCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGTGCTtga                                          <  2:195146/150‑1 (MQ=255)
                                                                                                           cgcgGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTtga                                          <  2:106492/150‑1 (MQ=255)
                                                                                                                    cgcCCGGCGTCATTTCGCTCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGGGCCGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCga                                 <  2:555860/150‑1 (MQ=255)
                                                                                                                                tttCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGGAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCtt                     >  2:955290/1‑150 (MQ=255)
                                                                                                                                tttCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCtt                     >  2:364368/1‑150 (MQ=255)
                                                                                                                                    ctcaaggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTtggcg                 >  2:633720/7‑150 (MQ=255)
                                                                                                                                        aggcCCTTGCCCGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTtggcgtggc             >  1:1305813/3‑150 (MQ=255)
                                                                                                                                             cTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGATTGGCGTGGCGTCAg        >  2:992834/1‑150 (MQ=255)
                                                                                                                                                   ggCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGTTTCaa  <  2:530208/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGCTCGATGTGCTGTTCGGCAAGAATCGGCGCACCGGCAGCGACCCGGCCGATGATGGGCAGGCCGGCTTCTTCAGCCTTGGCTTCCAGGCCAGGGATGCGGATGCCGCGGGAAGCGCCCGGCGTCATTTCGATCGCGCCCTTGCGGGCAAGGGCCTTGAGGTGCTCCTCGGCGGCATTGGGCGACTTGAAGCCCAGCTCCTGAGCAATCTCGGCGCGGGTAGGCGGGAAGCCGTTGTCTTCCAGGCAGCGCTTGATGAACGCGAGAATTTCGGCTTGGCGTGGCGTCAGTTTCAA  >  NC_002947/2445255‑2445551

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