Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,249 G→T 5.2% intergenic (+30/+27) PP_1425 → / ← nadB hypothetical protein/L‑aspartate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,625,2490GT 5.2% 301.8 / 3.1 96intergenic (+30/+27)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (39/52);  new base T (3/2);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GTTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCACTT  >  NC_002947/1625101‑1625397
                                                                                                                                                     |                                                                                                                                                    
gTTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGa                                                                                                                                                     <  1:1284928/150‑1 (MQ=255)
  tGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGc                                                                                                                                                   <  2:2277743/150‑1 (MQ=255)
      tGGGGGCAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                               >  1:1735240/1‑150 (MQ=255)
      tGGGGGAAAGCCAGGGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGGATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                               >  1:2135499/1‑150 (MQ=255)
         gggAAAGCCAGTCTATACCTCAGGAGGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGGGGTGGGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGaa                                                                                                                                            <  1:1001759/150‑1 (MQ=255)
                ccAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCacg                                                                                                                                     >  2:1796360/1‑150 (MQ=255)
                 cAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCacga                                                                                                                                    >  1:608453/1‑150 (MQ=255)
                 cAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCacga                                                                                                                                    >  1:2331249/1‑150 (MQ=255)
                    tgtATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGc                                                                                                                                 <  2:1211495/150‑1 (MQ=255)
                    tgtATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCACTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGc                                                                                                                                 <  2:647120/150‑1 (MQ=255)
                        taCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGat                                                                                                                             >  1:1708295/1‑150 (MQ=255)
                         aCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGCCGCCGtta                                                                                                                            <  2:2331249/150‑1 (MQ=255)
                          ccTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGatat                                                                                                                           >  1:1090598/1‑150 (MQ=255)
                           cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATc                                                                                                                          <  2:1708295/150‑1 (MQ=255)
                                  tGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCACAGATCGAGCCCGCGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCgg                                                                                                                   <  2:1557734/150‑1 (MQ=255)
                                       tAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTaa                                                                                                              >  2:834735/1‑150 (MQ=255)
                                       tAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTGGAAGTTCAGCCGGCGGCGCTGGGCTGTGGCGAGCTG‑TGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTaa                                                                                                              <  1:1570809/150‑1 (MQ=255)
                                         ggggCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAgg                                                                                                            >  1:1933720/1‑150 (MQ=255)
                                           ggCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGTTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGAt                                                                                                          <  1:341646/150‑1 (MQ=255)
                                            gCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                         >  2:802317/1‑150 (MQ=255)
                                            gCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGAGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                         <  1:540395/150‑1 (MQ=255)
                                              tGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGgt                                                                                                       <  2:1933720/150‑1 (MQ=255)
                                              tGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGAGGTGAGCTG‑TGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGgt                                                                                                       <  2:71100/150‑1 (MQ=255)
                                               gACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGgtg                                                                                                      <  1:662398/150‑1 (MQ=255)
                                                 cGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGT‑AGTCGTGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTc                                                                                                    <  1:1551929/150‑1 (MQ=255)
                                                 cGAGCAGCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTc                                                                                                    <  2:1420658/150‑1 (MQ=255)
                                                  gAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTcc                                                                                                   >  1:1958064/1‑150 (MQ=255)
                                                    gCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                         >  1:816836/1‑142 (MQ=255)
                                                    gCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                         <  2:816836/142‑1 (MQ=255)
                                                     cATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                                <  1:802317/150‑1 (MQ=255)
                                                     cATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                                <  1:72662/150‑1 (MQ=255)
                                                      aTCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGACGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTgg                                                                                               >  1:1255669/1‑150 (MQ=255)
                                                       tCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGc                                                                                              <  2:1260189/150‑1 (MQ=255)
                                                          gcgcTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTc                                                                                           >  2:780915/1‑150 (MQ=255)
                                                            gcTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGt                                                                                         >  2:931853/1‑150 (MQ=255)
                                                              ttGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTcc                                                                                       >  2:1374113/1‑150 (MQ=255)
                                                              ttGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTcc                                                                                       >  1:1973251/1‑150 (MQ=255)
                                                               tGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCa                                                                                      <  2:1735240/150‑1 (MQ=255)
                                                                gCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAg                                                                                     <  2:1973251/150‑1 (MQ=255)
                                                                        ggTTGAAGTTCAGCCGGCGGCGCTGGGCTGCTGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTg                                                                             <  2:548505/150‑1 (MQ=255)
                                                                         gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATc                                                                                                                          >  1:38792/1‑104 (MQ=255)
                                                                         gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATc                                                                                                                          <  2:38792/104‑1 (MQ=255)
                                                                         gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTgg                                                                            <  1:1374113/150‑1 (MQ=255)
                                                                         gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTgg                                                                            >  2:1831646/1‑150 (MQ=255)
                                                                          ttGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGa                                                                           <  2:1090598/150‑1 (MQ=255)
                                                                           tGAATTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGAt                                                                          <  1:931853/150‑1 (MQ=255)
                                                                             aaGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCAGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATaa                                                                        >  1:2185410/1‑150 (MQ=255)
                                                                               gTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                                      <  2:1748199/150‑1 (MQ=255)
                                                                               gTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTCCTTCGCGGGCTCGCCCGCGAAgaa                                                                                                                                         <  1:217582/83‑3 (MQ=38)
                                                                               gTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTCCTTCGCGGGCTCGCCCGCGAAgaa                                                                                                                                         >  2:217582/1‑81 (MQ=38)
                                                                                ttCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATcc                                                                     >  2:1033089/1‑150 (MQ=255)
                                                                                   aGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgt                                                                  <  1:2117474/150‑1 (MQ=255)
                                                                                   aGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgt                                                                  <  1:1070807/150‑1 (MQ=255)
                                                                                       ggcggcGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgtgt                                                              >  1:1201193/1‑150 (MQ=255)
                                                                                        gcggcgCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgtgt                                                              <  1:1647928/149‑1 (MQ=255)
                                                                                        gcggcgCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgtgt                                                              >  2:1647928/1‑149 (MQ=255)
                                                                                                   cTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGACCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTc                                                  >  1:948144/1‑150 (MQ=255)
                                                                                                   cTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAAAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCAAGTGTCTAATGCGAcccgc                                                  >  2:1661729/1‑148 (MQ=255)
                                                                                                     gCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCgg                                                <  2:1201193/150‑1 (MQ=255)
                                                                                                       ggTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgt                                                                <  1:2167541/132‑1 (MQ=255)
                                                                                                       ggTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgt                                                                >  2:2167541/1‑132 (MQ=255)
                                                                                                          gAGCCG‑GGCATAGCTCGAGACCGGGTTAGATTCTTTGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTTCAGCATCCCTGGATAATCCAGTGTGTAATGCAGCCCTCTGCTTTc                                           <  2:1255669/150‑1 (MQ=255)
                                                                                                          gAGCCG‑GGCATAGACCGAGCCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTc                                           <  1:696673/150‑1 (MQ=255)
                                                                                                            gCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCt                                         >  2:2117442/1‑150 (MQ=255)
                                                                                                               g‑ggCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTgc                                      <  1:780915/150‑1 (MQ=255)
                                                                                                                   cATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCTGCTTTCCTTGCGCt                                   >  1:1091867/1‑150 (MQ=255)
                                                                                                                      aGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCa                                <  1:1831646/150‑1 (MQ=255)
                                                                                                                       gATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATaa                                                                        <  2:2290693/109‑1 (MQ=255)
                                                                                                                       gATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATaa                                                                        >  1:2290693/1‑109 (MQ=255)
                                                                                                                        aTCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCAt                               <  1:2126913/149‑1 (MQ=255)
                                                                                                                        aTCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCAt                               >  2:2126913/1‑149 (MQ=255)
                                                                                                                         tCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                             <  1:1033089/150‑1 (MQ=255)
                                                                                                                         tCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCAGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                             <  2:2185410/150‑1 (MQ=255)
                                                                                                                         tCGAGACCGGGTTACCTTCTTAGCGAGCGGGCCCGGTAAGGCCACGACGCAGCTATCAGAGCGGGTTAGTTATTTTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                             <  1:1796360/150‑1 (MQ=255)
                                                                                                                         gcgggaCGGGGTTAGCTTCTTCGCGGGCGAGCACGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                             <  2:974326/146‑1 (MQ=255)
                                                                                                                          cGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGc                            <  1:2117442/150‑1 (MQ=255)
                                                                                                                                  ggTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGACACGACGCCGATATCAGAGCGGGTTAAGGATCGAGTAATTGAACTCGTCCAGCATCCCCCGATAATCGCGTGTGTTATGAAACCTTCGGCTTTCCTTGCGCGGAATGGCTGAAAGGa                    >  2:780293/1‑150 (MQ=255)
                                                                                                                                   gTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGat                   >  1:1741395/1‑150 (MQ=255)
                                                                                                                                      aGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGCTAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcat                <  2:1958064/150‑1 (MQ=255)
                                                                                                                                       gCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTTCGCTTCATGGCTGACAGGatcatc               >  2:2292376/1‑150 (MQ=255)
                                                                                                                                       gCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcatc               >  1:1359202/1‑150 (MQ=255)
                                                                                                                                        cttcttCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgt                                                                  <  2:920381/98‑1 (MQ=255)
                                                                                                                                        cttcttCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgt                                                                  >  1:920381/1‑98 (MQ=255)
                                                                                                                                         ttcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                         >  2:1359866/1‑58 (MQ=255)
                                                                                                                                         ttcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                         <  1:1359866/58‑1 (MQ=255)
                                                                                                                                         ttcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAg             >  1:1658675/1‑150 (MQ=255)
                                                                                                                                          tcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGc            <  2:169902/150‑1 (MQ=255)
                                                                                                                                           cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt           >  2:1867243/1‑150 (MQ=255)
                                                                                                                                           cttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                                      <  2:1637933/91‑1 (MQ=255)
                                                                                                                                           cttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                                      >  1:1637933/1‑91 (MQ=255)
                                                                                                                                             tcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTcc                                          >  2:1274573/1‑117 (MQ=255)
                                                                                                                                             tcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTcc                                          <  1:1274573/117‑1 (MQ=255)
                                                                                                                                               gcgGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGCATGGTGTCATTTGGCTCGTGCGGCAGCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGc       <  1:275709/150‑1 (MQ=255)
                                                                                                                                                   gCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAt                                                                       <  1:2134000/82‑1 (MQ=255)
                                                                                                                                                   gCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAt                                                                       >  2:2134000/1‑82 (MQ=255)
                                                                                                                                                    cGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGTCCTCGTCCACCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCTGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCACtt  <  2:608453/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCG‑GGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCACTT  >  NC_002947/1625101‑1625397

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