Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,764,407 G→A 5.9% A213V (GCC→GTC)  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,4070GA 5.9% 205.5 / 2.5 68A213V (GCC→GTC) PP_2418cobalamin ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

ATGAAGTCGGCGACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCAC  >  NC_002947/2764263‑2764555
                                                                                                                                                |                                                                                                                                                    
aTGAAGTCGGCGACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTc                                                                                                                                                 <  2:1015549/150‑1 (MQ=255)
  gAAGTCGGCGACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGTCGCGGGAGAACCAGAACATGGCCGAAGTTGACCACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAg                                                                                                                                               >  1:1191436/1‑150 (MQ=255)
  gAAGTCGGCGACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCACGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGCGGCGCCCTTGCAATTCGGCGtccgg                                                                                                                                               >  1:657740/1‑146 (MQ=255)
        ggCGACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTgccgc                                                                                                                                         >  1:863585/1‑150 (MQ=255)
           gACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTgccgcgcc                                                                                                                                      >  2:936722/1‑150 (MQ=255)
               cccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCgg                                                                                                                                  >  2:1269728/1‑150 (MQ=255)
                ccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGt                                                                                                                                 >  2:702196/1‑150 (MQ=255)
                ccTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGt                                                                                                                                 >  2:1373405/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:542376/1‑150 (MQ=255)
                             cacaTAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGaa                                                                                                                    >  2:542239/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:542239/150‑1 (MQ=255)
                                                                              gAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAAc                                                                   <  1:542376/150‑1 (MQ=255)
                                                                                  ttGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGTTCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCtt                                                               <  1:936722/150‑1 (MQ=255)
                                                                                      cGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAg                                                           >  1:725712/1‑150 (MQ=255)
                                                                                      cGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAg                                                           >  2:120094/1‑150 (MQ=255)
                                                                                        aCAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGc                                                         <  1:1373405/150‑1 (MQ=255)
                                                                                          aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                       <  2:1287869/150‑1 (MQ=255)
                                                                                          aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                       <  1:1203961/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                            <  2:1018202/127‑1 (MQ=255)
                                                                                                                                            ttCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTcc                            >  1:1018202/1‑127 (MQ=255)
                                                                                                                                            ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGccc                      <  1:876733/133‑1 (MQ=37)
                                                                                                                                            ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGccc                      >  2:876733/1‑133 (MQ=37)
                                                                                                                                               ggCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCAc  <  1:653570/150‑1 (MQ=255)
                                                                                                                                                gCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGa           <  2:1165334/140‑1 (MQ=255)
                                                                                                                                                gCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGa           >  1:1165334/1‑140 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATGAAGTCGGCGACACCCTGGCGGCCGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCAC  >  NC_002947/2764263‑2764555

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: