Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,764,422 C→G 7.2% G208A (GGC→GCC)  PP_2418 ← cobalamin ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,764,4220CG 7.2% 218.8 / 7.0 69G208A (GGC→GCC) PP_2418cobalamin ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base G (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTT  >  NC_002947/2764274‑2764570
                                                                                                                                                    |                                                                                                                                                    
gACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTgccgcgcc                                                                                                                                                     >  2:936722/1‑150 (MQ=255)
    cccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCgg                                                                                                                                                 >  2:1269728/1‑150 (MQ=255)
     ccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGt                                                                                                                                                >  2:1373405/1‑150 (MQ=255)
     ccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGt                                                                                                                                                >  2:702196/1‑150 (MQ=255)
     ccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGt                                                                                                                                                <  1:589597/150‑1 (MQ=255)
           ggccggccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGaa                                                                                                                                          >  1:1287869/1‑150 (MQ=255)
                 ccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTc                                                                                                                                               >  2:1190231/1‑139 (MQ=255)
                 ccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTc                                                                                                                                               <  1:1190231/139‑1 (MQ=255)
                  cacaTAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGaa                                                                                                                                   >  2:542239/1‑150 (MQ=255)
                  cacaTAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGaa                                                                                                                                   >  2:542376/1‑150 (MQ=255)
                     aTAGGGGTCGATGTGCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGAt                                                                                                                                <  1:567758/150‑1 (MQ=255)
                     aTAGGGGTCGATGACCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGAt                                                                                                                                >  1:291800/1‑150 (MQ=255)
                         gggTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCg                                                                                                                            <  2:1018350/150‑1 (MQ=255)
                          ggTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCgg                                                                                                                           <  2:78565/150‑1 (MQ=255)
                           gTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGAtc                                                                                                                               >  1:659419/1‑145 (MQ=255)
                           gTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGAtc                                                                                                                               <  2:659419/145‑1 (MQ=255)
                           gTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGc                                                                                                                          >  2:293554/1‑150 (MQ=255)
                              gATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAg                                                                                                                       >  1:378189/1‑150 (MQ=255)
                                   ccAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGc                                                                                                                  <  2:72205/150‑1 (MQ=255)
                                            cgcTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCt                                                                                                                 >  1:1484277/1‑142 (MQ=255)
                                            cgcTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCt                                                                                                                 <  2:1484277/142‑1 (MQ=255)
                                                 gagaACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAg                                                                                                    >  2:1203961/1‑150 (MQ=255)
                                                  agaACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                       <  1:1502863/114‑1 (MQ=255)
                                                  agaACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                       >  2:1502863/1‑114 (MQ=255)
                                                   gaACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGAt                                                                                                  <  2:288841/150‑1 (MQ=255)
                                                      ccAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCt                                                                                               >  1:126234/1‑150 (MQ=255)
                                                       cAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTg                                                                                              >  1:1093751/1‑150 (MQ=255)
                                                        aGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGt                                                                                             <  2:1093751/150‑1 (MQ=255)
                                                                   gAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAAc                                                                                  <  1:542376/150‑1 (MQ=255)
                                                                   gAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAAc                                                                                  <  1:542239/150‑1 (MQ=255)
                                                                       ttGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGTTCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCtt                                                                              <  1:936722/150‑1 (MQ=255)
                                                                           cGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAg                                                                          >  2:120094/1‑150 (MQ=255)
                                                                           cGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAg                                                                          >  1:725712/1‑150 (MQ=255)
                                                                             aCAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGc                                                                        <  1:1373405/150‑1 (MQ=255)
                                                                               aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                      <  1:1203961/150‑1 (MQ=255)
                                                                               aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                      <  2:1287869/150‑1 (MQ=255)
                                                                                 gTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGAt                                                                                                  >  2:282863/1‑120 (MQ=255)
                                                                                 gTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGAt                                                                                                  <  1:282863/120‑1 (MQ=255)
                                                                                    cTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGcc                                                                 >  2:2631/1‑150 (MQ=255)
                                                                                      tGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGt                                                               <  2:1191436/150‑1 (MQ=255)
                                                                                         cGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTg                                                            >  1:1358278/1‑150 (MQ=255)
                                                                                            ccAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCAc                                                         <  1:2631/150‑1 (MQ=255)
                                                                                              aGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCg                                                       >  2:1280508/1‑150 (MQ=255)
                                                                                                  gtgCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCa                                                   <  1:293554/150‑1 (MQ=255)
                                                                                                      ccTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGAt                                                                                                  <  2:364734/99‑1 (MQ=255)
                                                                                                      ccTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGAt                                                                                                  >  1:364734/1‑99 (MQ=255)
                                                                                                      ccTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGtt                                               >  1:109082/1‑150 (MQ=255)
                                                                                                      ccTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGtt                                               <  2:126234/150‑1 (MQ=255)
                                                                                                      ccTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGtt                                               >  2:1080217/1‑150 (MQ=255)
                                                                                                       cTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTg                                              >  2:1476656/1‑150 (MQ=255)
                                                                                                         ggCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGTACCGACCAGGTTGTc                                            <  1:874681/150‑1 (MQ=255)
                                                                                                          gcgcTGTCGAGGCGCCCTTGCAATGCGGCGTTCAGTGCCGCGCCCCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTcc                                           <  1:758433/150‑1 (MQ=255)
                                                                                                           cgcTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCt                                          <  1:1269728/150‑1 (MQ=255)
                                                                                                               ccgatctcgcCCTTGCAATTCGGCGTTCAGTGCCGCGCAGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGcc                                      <  1:1250676/143‑1 (MQ=255)
                                                                                                                      cgcCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAg                               <  1:702196/150‑1 (MQ=255)
                                                                                                                       gcCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTAACTGGAACGCTTGGAGGCCATTACCTACATCGCCGTCCTGGTTGTCCTTGcccagatcgg                              >  2:1250676/1‑143 (MQ=255)
                                                                                                                              cAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATgg                       <  1:120094/150‑1 (MQ=255)
                                                                                                                               tcttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCgcc                                                              >  2:1123192/3‑108 (MQ=21)
                                                                                                                               tcttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTc                                                              <  1:1123192/108‑1 (MQ=21)
                                                                                                                                aTTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGgca                     >  1:574680/1‑150 (MQ=255)
                                                                                                                                 ttCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTcc                                           >  1:1018202/1‑127 (MQ=255)
                                                                                                                                 ttCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTcc                                           <  2:1018202/127‑1 (MQ=255)
                                                                                                                                 ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGccc                                     <  1:876733/133‑1 (MQ=37)
                                                                                                                                 ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGccc                                     >  2:876733/1‑133 (MQ=37)
                                                                                                                                    ggCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCAc                 <  1:653570/150‑1 (MQ=255)
                                                                                                                                     gCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGa                          >  1:1165334/1‑140 (MQ=255)
                                                                                                                                     gCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGa                          <  2:1165334/140‑1 (MQ=255)
                                                                                                                                     acGTTCAGGACCGCGGCGCGGCACTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACg                >  1:877518/2‑150 (MQ=255)
                                                                                                                                            gTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTgg         <  2:291800/150‑1 (MQ=255)
                                                                                                                                            gTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTgg         <  2:1358278/150‑1 (MQ=255)
                                                                                                                                               ccgcgccgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGa      <  2:435030/150‑1 (MQ=255)
                                                                                                                                               ccgcgccgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGa      >  2:1359687/1‑150 (MQ=255)
                                                                                                                                                   gccgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCtt  <  2:863585/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTT  >  NC_002947/2764274‑2764570

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: