Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,340 C→G 5.3% P141R (CCG→CGG)  yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,3400CG 5.3% 237.8 / 4.4 75P141R (CCG→CGG) yrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base G (2/2);  total (35/40)
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.61e-01

AAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATG  >  NC_002947/4353192‑4353488
                                                                                                                                                    |                                                                                                                                                    
aaGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCg                                                                                                                                                     >  1:1272920/1‑150 (MQ=255)
   cAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAg                                                                                                                                                  <  2:1072977/150‑1 (MQ=255)
   cAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAg                                                                                                                                                  >  1:1449632/1‑150 (MQ=255)
    aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAgg                                                                                                                                                 >  1:1181837/1‑150 (MQ=255)
     ggCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGc                                                                                                                                                >  1:702108/1‑150 (MQ=255)
     ggCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGc                                                                                                                                                >  1:1227242/1‑150 (MQ=255)
      gCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGcc                                                                                                                                               >  2:474795/1‑150 (MQ=255)
          ttGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaa                                                                                                                                           <  2:1609008/150‑2 (MQ=255)
          ttGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAAc                                                                                                                                           <  1:1181021/150‑1 (MQ=255)
               ccTTATTACAGCGAAGTGGGCGCCTTTTTCCCGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAActgctg                                                                                                                                      <  1:1802519/150‑1 (MQ=255)
               ccTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaatggaa                                                                                                                                      <  1:1503387/150‑7 (MQ=255)
                cTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACtgctgc                                                                                                                                     >  2:725526/1‑150 (MQ=255)
                   attaCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCCGAGATGCCTGCTGGTCGAACCTGTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGc                                                                                                                                  <  1:1672252/150‑1 (MQ=255)
                     taCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                <  1:335260/150‑1 (MQ=255)
                     taCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                <  1:1813467/150‑1 (MQ=255)
                     taCAGCGAAGTGGGCGCCGATCTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                <  1:588793/150‑1 (MQ=255)
                        aGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTa                                                                                                                             >  2:430860/1‑150 (MQ=255)
                         gCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTaa                                                                                                                            >  2:721922/1‑150 (MQ=255)
                         gCGAAGGGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGACGCGTTGCGCCGGGAGATGGGAAGCTGCCATTTCGGCCTGCCGCAGGCCGAACCGCTGCAGCgagcaa                                                                                                                            >  2:940310/1‑145 (MQ=255)
                               tGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCa                                                                                                                      >  2:1671694/1‑150 (MQ=255)
                               ggggCGCGGTTTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCa                                                                                                                      <  1:1848408/149‑1 (MQ=255)
                                gggCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAg                                                                                                                     <  2:717788/150‑1 (MQ=255)
                                 ggCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGc                                                                                                                    <  1:1319976/150‑1 (MQ=255)
                                  gcgcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCg                                                                                                                   <  1:725526/150‑1 (MQ=255)
                                    gcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGt                                                                                                                 >  2:1420015/1‑150 (MQ=255)
                                          ttCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTGCGACGAGCAGAGCAGACGGCTGCTCTTAGCGTTTCCCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAg                                                                                                           <  1:1468552/150‑1 (MQ=255)
                                                  gcCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTc                                                                                                   <  2:899765/150‑1 (MQ=255)
                                                  ccccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTc                                                                                                   <  2:1227242/149‑1 (MQ=255)
                                                    ccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCNTGCTNGTCGAAGCGTTGCGCCCGGCAGTGGTGAGCTTCCATTTCGGCCTGCCACAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGAAATGTCCa                                                                                                 >  1:14075/1‑150 (MQ=255)
                                                                 cAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCgg                                                                                    >  2:1735257/1‑149 (MQ=255)
                                                                  aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGCGTAAAGGCCAGCGGTGCCAAGGTACTGTCCATCGCCACGACCGTg                                                                                   <  1:832322/150‑1 (MQ=255)
                                                                   aCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTgg                                                                                  <  2:1272920/150‑1 (MQ=255)
                                                                   aCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTgg                                                                                  <  1:940310/150‑1 (MQ=255)
                                                                    ccGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGCCCGTGga                                                                                 <  2:614525/150‑1 (MQ=255)
                                                                     cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGt                                                                                    <  2:1244831/146‑1 (MQ=255)
                                                                     cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGt                                                                                    >  1:1244831/1‑146 (MQ=255)
                                                                             tttCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                        >  1:1241374/1‑150 (MQ=255)
                                                                             tttCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                        >  1:884930/1‑150 (MQ=255)
                                                                                 gacgaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctgg                                                                    <  1:1915199/150‑1 (MQ=255)
                                                                                    gaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggctg                                                                 >  1:1854856/1‑150 (MQ=255)
                                                                                      gCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                               <  1:474795/150‑1 (MQ=255)
                                                                                        agagCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAgc                                                             >  2:1649018/1‑150 (MQ=255)
                                                                                             tgcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAAt                                                        <  1:1649018/150‑1 (MQ=255)
                                                                                              gcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATg                                                       <  1:1420015/150‑1 (MQ=255)
                                                                                               cctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATgg                                                      >  2:426709/1‑150 (MQ=255)
                                                                                               cctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACAGTTCCGAAAGCGGCCTGGCTGGAGCGCAATgg                                                      <  1:1483927/150‑1 (MQ=255)
                                                                                                ctgctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                        <  1:1890602/99‑1 (MQ=255)
                                                                                                ctgctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                        >  2:1890602/1‑99 (MQ=255)
                                                                                                  gctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtg                                                   >  2:303016/1‑150 (MQ=255)
                                                                                                   ctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                  <  1:1514210/150‑1 (MQ=255)
                                                                                                   ctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGACAAGGTCCTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                  >  2:1697553/1‑150 (MQ=255)
                                                                                                     ggTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcga                                                <  1:426709/150‑1 (MQ=255)
                                                                                                         gAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgcg                                            >  1:1556211/1‑150 (MQ=255)
                                                                                                          aaGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgcga                                           <  2:1449632/150‑1 (MQ=255)
                                                                                                           aGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCATGAACTGTCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTgcgatgcgat                                          <  1:1671694/150‑1 (MQ=255)
                                                                                                               ttGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGatcatc                                      <  2:702108/150‑1 (MQ=255)
                                                                                                                     ccGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATg                                >  1:1754958/1‑150 (MQ=255)
                                                                                                                      cGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATgg                               <  2:678936/150‑1 (MQ=255)
                                                                                                                         aGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGt                            <  1:1200716/150‑1 (MQ=255)
                                                                                                                          ggtggtGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTa                           >  1:841012/1‑150 (MQ=255)
                                                                                                                              gtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGaa                       >  2:291838/1‑150 (MQ=255)
                                                                                                                                gAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGAAACGATCATCGCCATGGGGTATGAAGc                     <  1:1697553/150‑1 (MQ=255)
                                                                                                                                    ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGt                                                                                    >  2:1774588/1‑83 (MQ=255)
                                                                                                                                    ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGt                                                                                    <  1:1774588/83‑1 (MQ=255)
                                                                                                                                     tCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                <  1:721922/150‑1 (MQ=255)
                                                                                                                                     tCCATTGCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGTCTGGCTGGAGCGCAATGGCGGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                <  2:732942/150‑1 (MQ=255)
                                                                                                                                        aTTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGcc             >  1:1726880/1‑150 (MQ=255)
                                                                                                                                         tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa            >  2:1923167/1‑150 (MQ=255)
                                                                                                                                         tttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa            >  1:253529/1‑150 (MQ=255)
                                                                                                                                          ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCAt           >  2:1624172/1‑150 (MQ=255)
                                                                                                                                          ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCAt           >  1:1650716/1‑150 (MQ=255)
                                                                                                                                            cGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATcg         <  1:430860/150‑1 (MQ=255)
                                                                                                                                                cTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgcg                                            >  1:599274/1‑111 (MQ=255)
                                                                                                                                                cTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgcg                                            <  2:599274/111‑1 (MQ=255)
                                                                                                                                                cTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGc     >  2:1274530/1‑150 (MQ=255)
                                                                                                                                                   ccGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcga                                                <  1:493551/104‑1 (MQ=255)
                                                                                                                                                   ccGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcga                                                >  2:493551/1‑104 (MQ=255)
                                                                                                                                                   ccGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATg  <  2:1241374/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATG  >  NC_002947/4353192‑4353488

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