Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,982,858 C→G 7.3% G1149G (GGC→GGG PP_5242 → GAF domain‑containing sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,982,8580CG 7.3% 170.8 / 4.9 55G1149G (GGC→GGGPP_5242GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (2/2);  total (29/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 = 7.08e-01

AACACCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGCTGCCGGCTTCATCCACGAGCAGGCGCTG  >  NC_002947/5982709‑5983005
                                                                                                                                                     |                                                                                                                                                   
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCTGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGc                                                                                                                                                     <  2:126450/150‑1 (MQ=255)
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGc                                                                                                                                                     <  1:482417/150‑1 (MQ=255)
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGc                                                                                                                                                     <  2:1277800/150‑1 (MQ=255)
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGc                                                                                                                                                     <  2:432774/150‑1 (MQ=255)
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGc                                                                                                                                                     >  1:432774/1‑150 (MQ=255)
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTGGTCTACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGc                                                                                                                                                     <  1:175602/150‑1 (MQ=255)
aacaCCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACAGCCGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCATGCGAGAGGGATTCAGCCGCTGGGCAGGc                                                                                                                                                     >  2:621743/1‑150 (MQ=255)
  cacCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCtt                                                                                                                                                   <  1:1209863/150‑1 (MQ=255)
  cacCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCtt                                                                                                                                                   <  2:285433/150‑1 (MQ=255)
   acCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTc                                                                                                                                                  >  1:483869/1‑150 (MQ=255)
             cACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCa                                                                                                                                        >  2:1201916/1‑150 (MQ=255)
                    gcACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGGTTGAAGCCGGcccagcg                                                                                                                                    >  1:809660/1‑142 (MQ=255)
                    gcACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGGTTGAAGCCGGcccagcg                                                                                                                                    <  2:809660/147‑6 (MQ=255)
                        aTCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCAACTTCCgcagc                                                                                                                             >  1:84276/1‑150 (MQ=255)
                        aTCAACCTGACCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCgcagc                                                                                                                             >  2:1137464/1‑150 (MQ=255)
                          cAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCagcaa                                                                                                                           <  2:1133241/150‑1 (MQ=255)
                              cTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAg                                                                                                                       >  1:15532/1‑150 (MQ=255)
                               tGGTCGACTATCACCTGTTAAGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCAACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGt                                                                                                                      >  2:213734/1‑150 (MQ=255)
                                         tCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGc                                                                                                                                          >  2:853003/1‑120 (MQ=255)
                                         tCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGc                                                                                                                                          <  1:853003/120‑1 (MQ=255)
                                         tCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGa                                                                                                            >  2:627140/1‑150 (MQ=255)
                                         tCACCTGTTCCGCGGCATGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGa                                                                                                            >  2:1417222/1‑150 (MQ=255)
                                            ccTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAAt                                                                                                               >  1:1277680/1‑144 (MQ=255)
                                            ccTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAAt                                                                                                               <  2:1277680/144‑1 (MQ=255)
                                            ccTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCAACTTCCGCAGCAAGCAGTTGCTAATCGAggg                                                                                                         <  2:84276/150‑1 (MQ=255)
                                             cTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGc                                                                                                        <  1:266526/150‑1 (MQ=255)
                                             cTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGc                                                                                                        <  1:1417222/150‑1 (MQ=255)
                                                ttCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTg                                                                                                     >  1:514988/1‑150 (MQ=255)
                                                            gCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTTATCGCCAAGCTc                                                                                         <  2:15532/150‑1 (MQ=255)
                                                               cTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGaa                                                                                      <  2:20805/150‑1 (MQ=255)
                                                                   gCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCg                                                                                  <  1:1201916/150‑1 (MQ=255)
                                                                      gCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACg                                                                               <  1:1137464/150‑1 (MQ=255)
                                                                       ccccGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACgg                                                                              >  2:39012/1‑150 (MQ=255)
                                                                        cccGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCCGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGc                                                                             <  2:1342975/150‑1 (MQ=255)
                                                                              gCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGc                                                                                                <  2:535269/125‑1 (MQ=255)
                                                                              gCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGc                                                                                                >  1:535269/1‑125 (MQ=255)
                                                                               cccACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCcg                                                                      >  2:774840/1‑150 (MQ=255)
                                                                                               cGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGGCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGa                                                      >  2:1023829/1‑150 (MQ=255)
                                                                                               cGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTTCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGGTCGAAGGCGACTGCCTGTCCGCGATCCGCTGCTTCGa                                                      <  2:234806/150‑1 (MQ=255)
                                                                                                   aaGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGaa                                                                                      <  1:1145721/114‑1 (MQ=255)
                                                                                                   aaGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGaa                                                                                      >  2:1145721/1‑114 (MQ=255)
                                                                                                     gCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGc                                                                                                                   >  2:337570/1‑83 (MQ=255)
                                                                                                     gCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGc                                                                                                                   <  1:337570/83‑1 (MQ=255)
                                                                                                            gAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTTATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATc                                         <  2:151451/150‑1 (MQ=255)
                                                                                                                tGCATGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATcgccg                                     >  2:401723/1‑150 (MQ=255)
                                                                                                                    ggCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGCGGCAAATctcctcgac                                 <  2:994740/150‑6 (MQ=255)
                                                                                                                      cccTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACg                                                                               <  2:1336468/102‑1 (MQ=255)
                                                                                                                      cccTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACg                                                                               >  1:1336468/1‑102 (MQ=255)
                                                                                                                       ccTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGc                              >  2:44204/1‑150 (MQ=255)
                                                                                                                           cgcgAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATc                                                                                                              >  1:551564/1‑66 (MQ=255)
                                                                                                                           cgcgAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATc                                                                                                              <  2:551564/66‑1 (MQ=255)
                                                                                                                               agcgaTTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGCTGCCGGCt                      <  1:39012/150‑1 (MQ=255)
                                                                                                                                    ttCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCt                                                                           <  2:45837/92‑1 (MQ=255)
                                                                                                                                    ttCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCt                                                                           >  1:45837/1‑92 (MQ=255)
                                                                                                                                         gaaggtggCCGGGTTGAAGCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATcgcc                                      <  1:1281355/118‑1 (MQ=255)
                                                                                                                                         gaaggtggCCGGGTTGAAGCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATcgcc                                      >  2:1281355/7‑124 (MQ=255)
                                                                                                                                            cTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGCTGCCGGCTTCATCCACGAGCa         >  2:701520/1‑150 (MQ=255)
                                                                                                                                               ggCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGCTGCCGGCTTCATCCACGAGCAGgc      >  2:753410/1‑150 (MQ=255)
                                                                                                                                                 ccGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCa                                                   <  2:482268/103‑1 (MQ=255)
                                                                                                                                                 ccGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCa                                                   >  1:482268/1‑103 (MQ=255)
                                                                                                                                                 ccGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGCTGCCGGCTTCATCCACGAGCAGgcgc    >  1:897691/1‑150 (MQ=255)
                                                                                                                                                   ggCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGCGCTCGGTAAGCTGCTTCGACCAGTCGCAAATCGCCGCACCCGCTGCCGGCTTCATCCACGAGCAGGCGCTg  >  1:738186/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACACCTGGGCCGCACCGGCGCACATCAACCTGGCCGACTATCACCTGTTCCGCGGCCTGGCCCTGGGCAGCCCCGAGGCCCACGGCAGCCTTGCCGACAAGCTGCGCGAGCTGCAGGCCCTGCGCGAGCGATTCAGCCGCTGGGCCGGCTTCAACCCGGCCACCTTCCGCAGCAAGCAGTTGCTAATCGAGGGCGTGATCGCCAAGCTCGAAGGCGACGGCCTGGCCGCGATCCGCTGCTTCGACCAGGCGCAAATCGCCGCCACCGCTGCCGGCTTCATCCACGAGCAGGCGCTG  >  NC_002947/5982709‑5983005

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: