Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,872,278 T→A 5.4% L174Q (CTG→CAG)  cobQ → adenosyl‑cobyric acid synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,872,2780TA 5.4% 293.8 / 5.9 93L174Q (CTG→CAG) cobQadenosyl‑cobyric acid synthase
Reads supporting (aligned to +/- strand):  ref base T (42/45);  new base A (2/3);  total (44/49)
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.35e-01

CCGGTGGTCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTACTGGGC  >  NC_002947/1872133‑1872420
                                                                                                                                                 |                                                                                                                                              
ccGGTGGTCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCtggtg                                                                                                                                            >  2:328377/1‑150 (MQ=255)
ccGGTGGTCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACACGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCtggtt                                                                                                                                            >  2:293376/1‑149 (MQ=255)
   gtggtCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGc                                                                                                                                         >  1:1112039/1‑150 (MQ=255)
    tggtCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCAGGTGGGcg                                                                                                                                        <  2:831100/150‑2 (MQ=255)
      gtCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACg                                                                                                                                      >  1:1490598/1‑150 (MQ=255)
      gtCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACg                                                                                                                                      >  2:1026559/1‑150 (MQ=255)
      gtCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACg                                                                                                                                      >  1:2031727/1‑150 (MQ=255)
       tCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGc                                                                                                                                     >  1:2397818/1‑150 (MQ=255)
       tCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCAGGTGGGcga                                                                                                                                       <  1:743084/148‑3 (MQ=255)
       tCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCAGGTGGGcga                                                                                                                                       >  2:743084/1‑146 (MQ=255)
         agggTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACACGGGCTTTGCCGAAGCGGCGGACTGTCCGGTGAGCCTGGTGGCCGACAGCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTg                                                                                                                                   <  2:274183/148‑1 (MQ=255)
           ggTGGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGCTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGa                                                                                                                                 <  2:403000/150‑1 (MQ=255)
             tCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAAt                                                                                                                               >  2:134926/1‑150 (MQ=255)
             tCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAAt                                                                                                                               >  2:2057539/1‑150 (MQ=255)
               gAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTg                                                                                                                             >  1:968250/1‑150 (MQ=255)
                  ggggCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTg                                                                                                                          >  1:679302/1‑150 (MQ=255)
                   gggCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCGCCTGGTGGGCACGCTGGAATTGCTGt                                                                                                                         <  2:2086234/150‑1 (MQ=255)
                   gggCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGt                                                                                                                         <  2:1112039/150‑1 (MQ=255)
                         ggTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCa                                                                                                                   <  2:2031727/150‑1 (MQ=255)
                           tCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACt                                                                                                                 <  2:1444893/150‑1 (MQ=255)
                           tCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGGGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACt                                                                                                                 >  2:853514/1‑150 (MQ=255)
                               cGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAAc                                                                                                             >  1:330255/1‑150 (MQ=255)
                                 gCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAg                                                                                                           <  1:1817921/150‑1 (MQ=255)
                                   agagATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGgc                                                                                                         >  2:186684/1‑150 (MQ=255)
                                      gaTCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCAGGTGGGcgaa                                                                                                                                      <  2:804407/118‑4 (MQ=255)
                                      gaTCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACAGCAACCGCGGTGGTGTGTTCGCCCACCAGGTGGGcgaa                                                                                                                                      >  1:804407/1‑115 (MQ=255)
                                       aTCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCgg                                                                                                     <  2:459827/150‑1 (MQ=255)
                                         cAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGt                                                                                                   <  2:679302/150‑1 (MQ=255)
                                         cAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGt                                                                                                   <  1:2323290/150‑1 (MQ=255)
                                          ccgtcgCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCgg                                                                                                     >  2:504360/6‑147 (MQ=255)
                                          ccgtcgCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCgg                                                                                                     <  1:504360/142‑1 (MQ=255)
                                           aCCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATGTCTTTCGACCACTGAACAGGCGCGGGTCa                                                                                                 >  2:266405/1‑150 (MQ=255)
                                             cTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCaaa                                                                                               >  1:1073731/1‑150 (MQ=255)
                                             cTGCGTGCAGGAGAAATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCACGGTGATCCTGGTGGCCGACATCACCCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCaaa                                                                                               <  2:236570/150‑1 (MQ=255)
                                                 gTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGCCGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGtt                                                                                           <  1:1230842/150‑1 (MQ=255)
                                                      ggCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGtc                                                                                      <  1:1961790/150‑1 (MQ=255)
                                                      ggCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGtc                                                                                      >  1:1702890/1‑150 (MQ=255)
                                                       gCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGtca                                                                                     <  1:813443/150‑1 (MQ=255)
                                                        cGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGtcat                                                                                    <  2:1490598/150‑1 (MQ=255)
                                                          atatCGCCAACATGGGCTTTGCCGAAGCGGTGGACTTCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGtcatca                                                                                  >  2:878943/1‑150 (MQ=255)
                                                          atatCGACAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTg                                                                                                                                   >  1:668856/1‑101 (MQ=255)
                                                          atatCGACAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTg                                                                                                                                   <  2:668856/101‑1 (MQ=255)
                                                               gCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCCACCGc                                                                             >  1:176840/1‑150 (MQ=255)
                                                                ccAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCaaa                                                                                               >  1:594579/1‑131 (MQ=255)
                                                                ccAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCaaa                                                                                               <  2:594579/131‑1 (MQ=255)
                                                                 cAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCtt                                                                           >  1:546881/1‑150 (MQ=255)
                                                                     aTGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGt                                                                       >  1:1839356/1‑150 (MQ=255)
                                                                       gggCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTgg                                                                     <  2:1839054/150‑1 (MQ=255)
                                                                        ggCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGc                                                                    <  2:330255/150‑1 (MQ=255)
                                                                             tGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACAt                                                               <  1:328377/150‑1 (MQ=255)
                                                                               ccGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTg                                                             >  1:200554/1‑150 (MQ=255)
                                                                                cGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGc                                                            >  1:82373/1‑150 (MQ=255)
                                                                                 gAAGCGGTGGCCTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCAACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGct                                                           >  2:499744/1‑150 (MQ=255)
                                                                                 gAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGct                                                           <  1:1880113/150‑1 (MQ=255)
                                                                                 gAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGct                                                           >  1:372516/1‑150 (MQ=255)
                                                                                  aaGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGctc                                                          <  1:1026559/150‑1 (MQ=255)
                                                                                  aaGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGctc                                                          <  1:809689/150‑1 (MQ=255)
                                                                                        tGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCtgctgc                                                    >  1:2379453/1‑150 (MQ=255)
                                                                                         ggACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCa                                                   <  2:546881/150‑1 (MQ=255)
                                                                                          gACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAg                                                  >  2:297101/1‑150 (MQ=255)
                                                                                             tGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCa                                               <  2:1073731/150‑1 (MQ=255)
                                                                                             tGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCa                                               >  2:717885/1‑150 (MQ=255)
                                                                                                  ggTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCt                                          >  2:721762/1‑150 (MQ=255)
                                                                                                      ttCCTGGTGGCGGACATCAACCGCGGTGGTGTGTTCGCCCCCCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGCCATTGCTCTGCTGCAGCCAGGGCTGGAc                                      <  2:1260506/149‑1 (MQ=255)
                                                                                                         cTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAg                                                                                              >  2:55444/1‑91 (MQ=255)
                                                                                                         cTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAg                                                                                              <  1:55444/91‑1 (MQ=255)
                                                                                                          tggtggCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGActggc                                  <  2:968250/150‑1 (MQ=255)
                                                                                                          tggtggCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTACGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGActggc                                  >  1:231383/1‑150 (MQ=255)
                                                                                                              ggCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAggg                                            <  1:864812/136‑1 (MQ=255)
                                                                                                              ggCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAggg                                            >  2:864812/1‑136 (MQ=255)
                                                                                                               gCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCgg                                                                                                     <  2:1776167/78‑1 (MQ=255)
                                                                                                               gCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCgg                                                                                                     >  1:1776167/1‑78 (MQ=255)
                                                                                                               gCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGaa                             <  2:1758170/150‑1 (MQ=255)
                                                                                                                ccGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGaac                            >  1:707192/1‑150 (MQ=255)
                                                                                                                   aCATCAACCGCGATGGGGTGTTCGCCCACTGGGTGTGCACGGTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGGCAGCAACCGCTTCCGTGTCGACATTTCTCTGCTGCAGGCAGGGCTGGACTGCCTGGaacaac                         <  2:231383/150‑1 (MQ=255)
                                                                                                                    cATCGACGGCGGGGGTGTAGTCGCCCGCGTGGTGGGCACGCTCGAATTGCTGTCGACCACTGAGCAGGCGCGGGTCGAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACg                        <  2:176840/150‑1 (MQ=255)
                                                                                                                         aCCGCGGTGGTGTGTTCGCCCACCTGGTGGGCATGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTg                                                                                        >  2:465952/1‑81 (MQ=255)
                                                                                                                         aCCGCGGTGGTGTGTTCGCCCACCTGGTGGGCATGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTg                                                                                        <  1:465952/81‑1 (MQ=255)
                                                                                                                          ccGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCgg                  <  2:372516/150‑1 (MQ=255)
                                                                                                                          ccGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCgg                  <  1:266405/150‑1 (MQ=255)
                                                                                                                           cgcgGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGc                 <  1:816878/150‑1 (MQ=255)
                                                                                                                              ggtggtGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCaga              >  2:2096868/1‑148 (MQ=255)
                                                                                                                              ggtggtGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCaaa              >  1:642679/1‑150 (MQ=255)
                                                                                                                               gtggtgTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGActggctg                                >  2:516933/1‑131 (MQ=255)
                                                                                                                               gtggtgTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGActggctg                                <  1:516933/131‑1 (MQ=255)
                                                                                                                               gtggtgTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAAc             <  2:2379453/150‑1 (MQ=255)
                                                                                                                                 ggtgTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCg           <  1:134926/150‑1 (MQ=255)
                                                                                                                                  gtgtgtTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCaa                                                                                                >  1:1704647/1‑64 (MQ=255)
                                                                                                                                  gtgtgtTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCaa                                                                                                <  2:1704647/64‑1 (MQ=255)
                                                                                                                                       ttCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTAc       >  2:181581/1‑148 (MQ=255)
                                                                                                                                       ttCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTAc       <  1:181581/148‑1 (MQ=255)
                                                                                                                                       ttCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTACTg     >  2:1865491/1‑150 (MQ=255)
                                                                                                                                        tCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCCAAAGCTTCCGTGGCGACATTGGTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTACTgg    <  1:297101/150‑1 (MQ=255)
                                                                                                                                        tCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTACTgg    <  1:721762/150‑1 (MQ=255)
                                                                                                                                          gCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTACTGGGc  <  2:2397818/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCGGTGGTCATGGTCGAGGGGGCGGGGTCGCCGGCAGAGATCAACCTGCGTGCCGGCGATATCGCCAACATGGGCTTTGCCGAAGCGGTGGACTGCCCGGTGATCCTGGTGGCCGACATCAACCGCGGTGGTGTGTTCGCCCACCTGGTGGGCACGCTGGAATTGCTGTCGACCACTGAACAGGCGCGGGTCAAAGGTTTTGTCATCAACCGCTTCCGTGGCGACATTGCTCTGCTGCAGCCAGGGCTGGACTGGCTGGAACAACGCACCGGCAAACCGGTACTGGGC  >  NC_002947/1872133‑1872420

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