Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,338 T→G 16.2% intergenic (+55/‑88) PP_3338 → / → PP_3339 cytochrome bd‑type quinol oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,778,3380TG16.2% 168.2 / 4.0 90intergenic (+55/‑88)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base G (14/3);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 6.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACC  >  NC_002947/3778191‑3778486
                                                                                                                                                   |                                                                                                                                                     
tCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgc                                                                                                                                                     <  1:529642/150‑1 (MQ=255)
  gATCAGCATCTGGGAAGCCGCCGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGgcgc                                                                                                                                                   >  2:1724196/1‑150 (MQ=255)
     cAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTg                                                                                                                                                <  2:5255321/150‑1 (MQ=255)
      aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                               <  2:3367202/150‑1 (MQ=255)
       gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCATCGCCGTGGGGGGCCGGCg                                                                                                                                              >  1:4683930/1‑150 (MQ=255)
        cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGa                                                                                                                                             <  1:4523250/150‑1 (MQ=255)
         aTCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                            <  2:175704/150‑1 (MQ=255)
         aTCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGAc                                                                                                                                            >  1:5385602/1‑150 (MQ=255)
            tGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACaaa                                                                                                                                         >  1:393102/1‑150 (MQ=255)
             gggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCCGCaaaa                                                                                                                                        >  1:3971754/1‑150 (MQ=255)
             gggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGACaaaa                                                                                                                                        >  1:3971540/1‑150 (MQ=255)
              ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaaa                                                                                                                                       >  1:4005523/1‑150 (MQ=255)
              ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaaa                                                                                                                                       <  2:4153846/150‑1 (MQ=255)
               gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAg                                                                                                                                      <  1:1798683/150‑1 (MQ=255)
               gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAg                                                                                                                                      <  1:2244030/150‑1 (MQ=255)
               gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCGGGTAGGCACACTGTGCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAg                                                                                                                                      <  1:2491502/150‑1 (MQ=255)
                aaGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAgg                                                                                                                                     <  2:2246157/150‑1 (MQ=255)
                 aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                    <  2:273709/150‑1 (MQ=255)
                  gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                   <  1:3476342/150‑1 (MQ=255)
                   ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCAAAAACgcgc                                                                                                                                  >  2:1062234/1‑150 (MQ=255)
                    cGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGca                                                                                                                                 <  1:804699/150‑1 (MQ=255)
                     gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCAAAAAGGCGcac                                                                                                                                >  2:2094154/1‑150 (MQ=255)
                     gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATgcgccc                                                                                                                                >  2:4749036/1‑148 (MQ=255)
                     gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                >  2:4749035/1‑150 (MQ=255)
                     gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                <  1:91782/150‑1 (MQ=255)
                          gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCAGCGACAAAAAAGCGccccaccc                                                                                                                           >  1:5098587/1‑148 (MQ=255)
                               cccAGTGCAATGGTTTCAGGGTGGTTGGCACTCTGTTCAGCGTGCCCTGTTCGCTTGGCTGATCCGAGTTTTCATGTTCGCCGCCGGCTGGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACtt                                                                                                                      <  2:3431679/150‑1 (MQ=255)
                                ttagagCCATGGTTTATTGCTGTTTGGCACAGTGTTCTGCCTGCCCTGTTCGCTTGGGTTAGCGGAATTTTCATTTTGGCCGCGGGCTTGCCCAGCAATTGAGCCCTGGCGGTGGTGCGCGTGCCAAAAAAAGGCGCACCACGGAACTTg                                                                                                                     <  2:4731385/148‑1 (MQ=255)
                                ccAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAAGGCGCACCCCCGAACATg                                                                                                                     >  2:3395040/1‑150 (MQ=255)
                                  agagCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTa                                                                                                                   <  2:4525097/150‑1 (MQ=255)
                                   gagCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGGGGCAAAAAAGCGCACCACCGCACGTGTaa                                                                                                                  >  1:3673022/1‑149 (MQ=255)
                                    agCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCAGGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGGCACAAAGGCGCCCCACGGAACGTGTACa                                                                                                                 >  2:440329/1‑150 (MQ=255)
                                        aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAAGGCGCACCACCCAACATGTAAAG‑ccc                                                                                                            >  1:5360332/1‑150 (MQ=255)
                                        aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAACGTGTACAG‑ccc                                                                                                            >  2:420977/1‑150 (MQ=255)
                                                      gtcggCACACTGTACAGCCTGCCCCGTTCGGTGGGCTGATCCGCCTGTTCATGAGGGCCGCCGGCAGGCGCCGCCATTGCGCCCCCGCGGGGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGc                                                                                              <  2:4100898/147‑1 (MQ=255)
                                                        aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAGGCGCACCCCCGCCCATGTACAG‑CCCGCCACCGATGCGGcgt                                                                                            >  2:1610419/1‑148 (MQ=255)
                                                         ggCACACTTTTCAGCGTGCCCTGTTCGTTTGGCTTATCCGACTGTTCATGTTGGCCGCGGGATTGCGCAGCAATTGACCCCTCGCCGTGGTGCGCCTGCCAAAAACAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTa                                                                                           <  2:689224/150‑1 (MQ=255)
                                                         ggCACACGGTTCAGCCTGCCCTGTTCGCTTGGCTTATCCGAGTTTTCATTTTGGCCGCGGGCTTGCGCAGCATTTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTa                                                                                           <  2:689280/150‑1 (MQ=255)
                                                            acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACGTGTTAAGCCCCCCAACGG‑TGCGGGTTAATg                                                                                        >  1:3123925/1‑150 (MQ=255)
                                                            acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACATGTACAG‑CCCCCCCCCGATGCCGCTGAATg                                                                                        >  2:5207281/1‑150 (MQ=255)
                                                            acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCCAAAAGGCGCCCCACCCGACGTGTTAAG‑CCCGCAACCGCTGCGGCTTAATg                                                                                        >  1:295609/1‑150 (MQ=255)
                                                            acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAAGCGCCCCACGGCGCGTGTAAAG‑CCCCCCACCGCTGCGGCGGAATg                                                                                        >  1:295481/1‑150 (MQ=255)
                                                             cacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCAAAAAGGCGCCCCACGGACCGTGTACAG‑CCCGCAACGGAAGCGGCTTAAttc                                                                                       >  2:4165831/1‑148 (MQ=255)
                                                                tGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGGCAAAAAGGCGCCCCACCGAACATGTACAG‑CCCCCAACCGATGCGGCGTAATgcgcc                                                                                    >  2:2641421/1‑149 (MQ=255)
                                                                 gTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAAGCGCACCACCGCACATTTAAAG‑CCCCCAACCGGTGGGGGATAAAgcgcgc                                                                                   >  2:2427176/1‑150 (MQ=255)
                                                                 gTTCACCCGGCCCAGGGCGCTGGGCCGATCCGACGGTTCCTGATGGCCGCCGGCATGCGCAGCCATTGCGCCCGCGCCTTGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATgcgcgc                                                                                   <  2:4114303/150‑1 (MQ=255)
                                                                     aTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAATGCGCACCACGGAAAATTTAAAG‑CCCGCCACGGAAGAGGGTTAATGCGCGCCata                                                                               >  2:2874043/1‑150 (MQ=255)
                                                                       cgtGCCCCGTGCGCTTGGCTGATCCGAGTGTTCATGTTGGCCGCCGGCATGCGCAGCCTTTGAGCCCGCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATAtt                                                                             <  2:1144361/148‑1 (MQ=255)
                                                                       ccTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCCCCACCGCACTTGTACAG‑CCCGCAACGGATGCGGCTTAATTCGCGCCAAAtt                                                                             >  2:5239782/1‑150 (MQ=255)
                                                                       ccTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAAGCGCACCACGGAACATGTACAC‑CCCGCCACGGAGTCGGCATAATGCGCGCCAAAtt                                                                             >  2:5239801/1‑150 (MQ=255)
                                                                         tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCCCACCACCGAACTTGTAAAA‑CCCCCAACGGAAGCGGCTTAATGCGCGCCctaatag                                                                           >  2:2244691/1‑146 (MQ=255)
                                                                          gCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGGGTCAAAAAGGCGCACGACGGTACTTGTTGAT‑CCCGCTACGGATGCGGCTTAATTCGCGGCCTAttggg                                                                          >  2:4082518/1‑147 (MQ=255)
                                                                             ctgTTCGGTTGGCTTAGCCGACTGTTCATGTTGGCCGCCGGCTTGCGCAGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCg                                                                       <  2:4621023/148‑1 (MQ=255)
                                                                             ccgTTCGCTTGGCTGATCCGCGTGTTCTTGTGGGCCGCCGGCTTGCGCCGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCg                                                                       <  2:4621025/148‑1 (MQ=255)
                                                                              aGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACCAAAAGGCGCACCCCCGGCCGTGTAAAC‑CCCGCCACCGATGCGGCTTAAATCGCGCCCAATTAGCGCgg                                                                      >  1:5323207/1‑150 (MQ=255)
                                                                                ttCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCAAAAAAGCGCCCCCCCGCCCATGTACAG‑CCCCCCACCCCGGCGGCGTTATGCGCGCCCAAGTCGCGCgggg                                                                    >  1:3454693/1‑149 (MQ=255)
                                                                                ttCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCAATGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACGGCACAAGTACAG‑ACCCCAACGGAAGACGCATAAAGCCCGACAAATTCGCTCGGGt                                                                    >  2:1727134/1‑150 (MQ=255)
                                                                                     ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGGCTGGGTGAAAAATTGGCGCCACGGGAGTTTTTAAT‑CCCGCGACGGGTGGGGGTTTATGGGGGCCATTTTAGGTGGGGTgtgtg                                                               >  2:511560/1‑146 (MQ=255)
                                                                                      tGGCTTATCCGACTTTTAATGTTGGCCGCGGGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCCTGCCACAAAATGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCa                                                              <  1:5296077/150‑1 (MQ=255)
                                                                                      tGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTTCGACAAAAATGCGCACCACGGAACTTGTACAG‑GCCGCAACGGATGGGGCTTTATGGGCGCCATATTGGGTCGGGTgtgtga                                                              >  2:2663031/1‑145 (MQ=255)
                                                                                         cttaTCGGACTGTTCATGTTGGCCGCCGGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGc                                                           <  2:421388/147‑1 (MQ=255)
                                                                                          tGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGCCAAAAAGGCGCACCACGGAAATATTACAG‑CCCCCAACGGATGCGGCTTAATTGGCGCCAAATTAGCGCGGGGGAATCATGCt                                                          >  2:702456/1‑150 (MQ=255)
                                                                                           gATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAAGGCGCACCCCGGCAAATGTACAG‑CCCGCCACCGATGCGGCGTAATTCGCGCCAGATTAGCGCGGGTGGCGCATGcgt                                                         >  2:2015287/1‑148 (MQ=255)
                                                                                            ttcgttgtGTTCTTGCGGGCCGCGGGCTTGCCCAGCCTTTGCCCCCTGGCGGTGGTGCCCCTGCAAAAAAAAGCCCCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCttt                                                        <  1:5035771/143‑1 (MQ=255)
                                                                                                   tGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAATGCGCACCACCGAACGTGTAAAG‑CCCCCCACCGATGCGGCGTAATGGGCGCCCTATTAGCGCGGGTGGCGCATGCTTTcacgtgc                                                 >  2:1742420/1‑150 (MQ=255)
                                                                                                    gTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCCACAAAAAGGCGCCCCACCGCACGTGTACAG‑CCCCCCACCGCTGCGGCGTAAATCGCGCCAAATTAGCGCGGGTGGCGAATGCTTTgacgtgcc                                                >  2:2475753/1‑149 (MQ=255)
                                                                                                     ttCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTaaagtgcac                                               >  2:3641033/1‑150 (MQ=255)
                                                                                                        atgaggGCCGCCGGCATGCGCAGCCATTGAGCCCCCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                            <  1:4030551/150‑1 (MQ=255)
                                                                                                         tgttgGCCGGCGGCTTGCGCAGCCTTTGCGCCCGCGGGGGGGGGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                           <  1:4745637/147‑1 (MQ=255)
                                                                                                         tgtgggCCGCCGGCTTGCGCCGCATTTGAGCCCTCGCGGTGGTGCGCCTGCAAAAAAATGGCGAACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                           <  1:4745742/146‑1 (MQ=255)
                                                                                                         tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAAAGTGTACAG‑CCCCCCACGGCTGCGGCGTAATGCGCGCCCTATTCGCTCGGGTGGCTCAAGCTGTCacgggcacctgc                                           >  2:2758286/1‑150 (MQ=255)
                                                                                                         ggcgggCCGCCGGCTTGCGCAGCCTTTGTGCCCTCGCCGTGGTGCGCCTCCGAAAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                           <  1:3056706/146‑1 (MQ=255)
                                                                                                             ggCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCCAAAAGGCGCACCACGGCACTTGTACAG‑CCCGCAACGGATGCGGCGTAATGCGCGCCCAATTCGCTCGGGGGGCGCCAGCGTTCACGTGCCCGGGcccaa                                       >  2:1706675/1‑149 (MQ=255)
                                                                                                             ggCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACAt                                       <  2:4541479/150‑1 (MQ=255)
                                                                                                               ccgcggGCTTGCGCCGCCTTTGCGCCGTCGGGGTGTTGCGCCTGCCAAAAAACGGCGCCCCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCAcgtgcacgtgcaagcat                                     <  1:1353625/150‑6 (MQ=255)
                                                                                                               ccgccgGCATGCGCAGCCTTTGAGCCATCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCt                                     <  2:4683930/150‑1 (MQ=255)
                                                                                                                 gccgGCTTGCGCAGCATTTGAGCCCTCGCGGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                   <  1:4306539/150‑1 (MQ=255)
                                                                                                                 gccgGCATGCGAAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                   <  2:2492028/150‑1 (MQ=255)
                                                                                                                  ccgccctGCGCCCCCATTGCCCCCGCGCCGTGGTCCCCCTGCTCAAAAACGCCGCCCCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTg                                  <  1:2012121/144‑1 (MQ=255)
                                                                                                                  ccgGCATGCGCAGCCATTGCGCCCGCGCGGGGGTGCGCCTGCGACAAAAAGGCGCCCCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTg                                  <  1:5327465/150‑1 (MQ=255)
                                                                                                                   gggaaTGCGCAGCTTTTGCGCCCTCGCGGTGGTGCGCGTGCGAAAAAATGGCCCCCCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                 <  1:2094154/146‑1 (MQ=255)
                                                                                                                     gCATGCGCAGCCTTTGAGCCCTCGGGGTGGGGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                               <  1:2803297/150‑1 (MQ=255)
                                                                                                                           gcAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACGGAACTTGTAAAG‑CCCGCAACGGATTCGGCTTAATTGGCGCCCTATTAGCTCGGGTGTCTCATGCGTTCACGTGCACGGTCACAACATTTTGTGTTCgc                         >  2:1750787/1‑150 (MQ=255)
                                                                                                                                 aTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACCGAAAATGTAAAG‑CCCGCAACCGATGCGGCTTAAATCGCGCCAAATTAGCTCGGGTGACGCATGCTTTAACGTGCAAGTGCACATATTTTTATGTTAgccttcca                   >  2:4676297/1‑149 (MQ=255)
                                                                                                                                     cgCCCTCGCCGTGGGGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc               <  2:3325777/149‑1 (MQ=255)
                                                                                                                                     aGCCCTCGCCGTGTTGCGCCTGCTAAAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc               <  1:2663031/150‑1 (MQ=255)
                                                                                                                                      gCCCTGGCGGTGGTGCGCCTGCGACAAACTGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCg              <  1:4033975/150‑1 (MQ=255)
                                                                                                                                       cccTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAAAATTTAAAG‑CCCCCAACGGATGCGGCGTAATGCGCGCCAAATTAGCTCGGGTGACTCAAGCGTTAACCTGCACGTGCCCAAAATTGTCTGTTCGCCAGCCGCGcccg             >  2:2266435/1‑148 (MQ=255)
                                                                                                                                       cccGCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt             <  2:2835019/150‑1 (MQ=255)
                                                                                                                                         ctcggggtggtgCGCCTCCAAAAAAATGCCCCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc           <  1:4676297/144‑1 (MQ=255)
                                                                                                                                             ccGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGtt       >  2:1026249/1‑150 (MQ=255)
                                                                                                                                               gtggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa     >  1:1639609/1‑150 (MQ=255)
                                                                                                                                                tggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCaa    >  1:4899216/1‑150 (MQ=255)
                                                                                                                                                 ggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc   <  2:4358598/150‑1 (MQ=255)
                                                                                                                                                 ggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc   <  2:328829/150‑1 (MQ=255)
                                                                                                                                                 ggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc   >  2:861182/1‑150 (MQ=255)
                                                                                                                                                  gtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAcc  >  2:2155012/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAG‑CCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACC  >  NC_002947/3778191‑3778486

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