Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,105 T→G 11.0% intergenic (+69/+49) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,1050TG11.0% 384.4 / 6.7 175intergenic (+69/+49)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base T (84/64);  new base G (10/15);  total (96/79)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

GCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352956‑353254
                                                                                                                                                      |                                                                                                                                                     
  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                                                                     <  2:4633642/150‑1 (MQ=255)
   aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                                                                    <  1:48385/150‑1 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCAGTgg                                                                                                                                                   >  1:1000135/1‑150 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCGGCGATCTCTGTgg                                                                                                                                                   >  1:1000204/1‑150 (MQ=255)
      aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGa                                                                                                                                                 <  2:3189190/150‑1 (MQ=255)
         ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGa                                                                                                                                                 >  1:2906181/1‑147 (MQ=255)
         ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCg                                                                                                                                              >  2:561262/1‑150 (MQ=255)
            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                           <  2:5571105/150‑1 (MQ=255)
            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                           <  1:976789/150‑1 (MQ=255)
              cTTTGGTGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                                                         <  2:1439301/150‑1 (MQ=255)
               tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTg                                                                                                                                        >  2:6820838/1‑150 (MQ=255)
                ttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                                                                       >  2:3007345/1‑150 (MQ=255)
                 gggCGAGTACACAGTGCCGGAGGCCAAGGCCGACAGGGCTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCGTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGCGATCTCAGCAACCCTGTGGGAGCGGGCGTGcc                                                                                                                                      <  1:5512271/149‑1 (MQ=21)
                    cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                                                                   <  2:2932017/150‑1 (MQ=255)
                        tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCGCCCCCAGGGTTGGGGGGGTGTGGGGGGGAGCGGGCGGGGCCGCGAAg                                                                                                                               >  1:4736594/1‑150 (MQ=37)
                        tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                               <  2:1651681/150‑1 (MQ=255)
                         acacTGCGCCGTCGGCCACGGCATACTGGGGTTGGGGGGGGATGGGGGTTGGGTGATCCATTTACCGGGGCTGCTTTGCAGCCATTCGGGGGCACGCCGGCTCCCACAGGTTTGATGATTTCTCGGTGGGAGGGGGCACGCCcgctccca                                                                                                                              <  2:3320640/150‑6 (MQ=255)
                         acacAGGGCCGGTGGCCAAGGCATACTGGGCTTCGCGCGGGATCGTCGTTCGCTGAGCCATTGACCGGGGCTGGTTTGCAGCCCTTCGGGGGCACGCCGGCTCCCAGAGGTTTGCTGATATCTCGGGGGGAGCGGGCACGCCcgctccca                                                                                                                              <  2:6790414/150‑6 (MQ=255)
                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGc                                                                                                                            >  1:3361249/1‑150 (MQ=255)
                             aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                                                                          >  2:1303627/1‑150 (MQ=255)
                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCGGCGATCTCTGGGGGAGCGGGCGGGCCCGCGAAGGGCcgcg                                                                                                                        >  2:2095478/1‑150 (MQ=255)
                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCCGCCCTTCGCGGGCACGCCCGCTCCCCCCGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                                                                    <  2:1774892/150‑1 (MQ=255)
                                    gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACACGGTTGCGGAGATCACTGTGGGAGCGGGCGTGCCCCCCAAGGGCCCCGAAGCg                                                                                                                   >  1:3655548/1‑150 (MQ=38)
                                     aGGCCAAGGCCTACCTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                                                                  <  2:5896658/150‑1 (MQ=255)
                                       gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcc                                                                                                                >  1:456288/1‑150 (MQ=255)
                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCCGCCCCCGCCCACCGGGTTTGGGAGATGTCTGTGGGGGCGGGGGGGGCCGCGCAGGGCCGCGGAGCGGcccc                                                                                                              >  1:53493/1‑150 (MQ=255)
                                                tACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                                                                       <  1:6757438/150‑1 (MQ=255)
                                                tACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                                                                       >  1:1663445/1‑150 (MQ=255)
                                                  cATGGGTTCGCGCGGGATCGCCGTTCGCTGCTCCATTGACCGGGGCTGCTGTGCCGCCCTTCGCGGGCACGCCCGCACCCCCCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                                                                     <  2:5748181/150‑1 (MQ=38)
                                                  cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                                                                     >  2:2899171/1‑150 (MQ=255)
                                                  cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCGTc                                                                                                     >  2:2899235/1‑150 (MQ=255)
                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATAGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCACACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                               >  2:5479588/1‑150 (MQ=255)
                                                         tCGCGCGGGATCGTCGTTCGCTGATCCATTGCCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                                                              <  1:6110134/150‑1 (MQ=255)
                                                           gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGGTGAGATCTCTATGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGt                                                                                            >  1:515263/1‑150 (MQ=255)
                                                           gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTCCCACAGGTTAGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCCAAGGGCCGCGAAGCGGCCCCAAATGCTTCAACCGTGGt                                                                                            >  2:262923/1‑150 (MQ=39)
                                                              cgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCt                                                                                         <  2:456288/150‑1 (MQ=255)
                                                               gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                                                                        >  2:2800011/1‑150 (MQ=255)
                                                                        gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAACGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcag                                                                               >  2:4050931/1‑150 (MQ=255)
                                                                              tGATCCATTGACCGGGGCTGCTGTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                         <  2:2539398/150‑1 (MQ=255)
                                                                                aTCCATTGACCGGGGCTGCTTGGCCGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGt                                                                       <  2:53493/150‑1 (MQ=255)
                                                                                   cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCTTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                    <  2:3116924/150‑1 (MQ=255)
                                                                                   cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                    <  2:5086537/150‑1 (MQ=255)
                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGCGATCTCGGTGGGCGCGGGCGTGCCCCCGAAGGGCCGCGAAGCGGCCCCAAATGCTTAAACCGTGGTTCTTTCCGAGCAGCCCGTGGTAAAg                                                                   >  1:3264992/1‑150 (MQ=38)
                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAg                                                                   >  2:1122641/1‑150 (MQ=255)
                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTCCCCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCAACCCTGGTTCTTGCCGAGCAGCGCGTGGTAAAg                                                                   >  2:4834688/1‑150 (MQ=39)
                                                                                       gACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGGGTGCGGCGATCGCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCATTGCGTCAACCGTGGGTCTTGCCGAGCCGCGCGTGGTAAAGCTc                                                                >  1:1127446/1‑150 (MQ=25)
                                                                                        attggggCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCCCACAAGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                                                         >  1:6011315/4‑76 (MQ=17)
                                                                                        attggggCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCCCACAAGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                                                         <  2:6011315/73‑1 (MQ=17)
                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                               <  2:212317/150‑1 (MQ=255)
                                                                                         ccGGGGCTGCTGTGCAGCCCTTCGCGGGCACGCCCGCTCCCTCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCgg                                                              <  2:2673318/150‑1 (MQ=255)
                                                                                            gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGGTTGCTGGGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCGTGCCGAGCAGCGCGGGGTAAAGCTCGGTAt                                                           >  2:3019457/1‑150 (MQ=38)
                                                                                             ggcggcTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGTTTTGCTGAGATCCCTGTGGCAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATc                                                          <  1:1827511/150‑1 (MQ=255)
                                                                                                  ctTTGCCCCCCTGCGCGGGCACGCCCGCTCCCACGGGTTGGCTCCGATCTCGGGGGGCCCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:318308/150‑1 (MQ=38)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCGGCTCCCAATGGTTTGCTGAGATCTCTGGGGGAGGGGGCTTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:4383261/150‑1 (MQ=39)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTTAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:414711/150‑1 (MQ=255)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:6842061/150‑1 (MQ=255)
                                                                                                    ttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt                                                   <  1:3987326/150‑1 (MQ=255)
                                                                                                    tggccgCCCTTCGCGGGCACGCCCCCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt                                                   <  1:3526616/145‑1 (MQ=255)
                                                                                                        aGCCCTTCGCGGGCACGCCCGCTCCCACCGGGTTGCGGAGATCTCGGGGGGAGCGGGCGGGCCCCCGAAGGGCCGCGAAGCGGCCCCCATTGCGTCAACCGTGGGTCGTGCCGAGCCGCGCGGGGTAAAGCGCGGTATCCCCCAGTATcc                                               >  2:78468/1‑150 (MQ=255)
                                                                                                             ttCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccacc                                          >  2:4523687/1‑150 (MQ=255)
                                                                                                              tCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccacca                                         <  2:474616/150‑1 (MQ=255)
                                                                                                                    gCACGCCCGCTCCCACAGGTTTGCTGAGAACTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGt                                   <  1:5965111/150‑1 (MQ=255)
                                                                                                                      aCGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                 >  2:1818035/1‑150 (MQ=255)
                                                                                                                      aCGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                 >  1:3799057/1‑150 (MQ=255)
                                                                                                                      aCGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCg                                 >  1:6829534/1‑150 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGtt                               >  1:5771795/1‑150 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGtt                               >  2:2833583/1‑150 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTAt                             >  1:4097972/1‑150 (MQ=255)
                                                                                                                           cGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATc                            >  2:2863686/1‑150 (MQ=255)
                                                                                                                            gCTCCCACCGGTTGGCTGAGAACTCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCGTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                           <  1:404426/150‑1 (MQ=39)
                                                                                                                              tCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGTGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTg                         >  1:2512329/1‑150 (MQ=255)
                                                                                                                              tCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGGCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTg                         >  1:6488320/1‑150 (MQ=255)
                                                                                                                              tCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTg                         >  1:1933247/1‑150 (MQ=255)
                                                                                                                               cccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgc                        <  2:6138206/150‑1 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       <  1:2765677/150‑1 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  2:6576967/1‑150 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  2:1131493/1‑150 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       <  2:4814911/150‑1 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccac                                           >  2:3770360/1‑128 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  1:3582678/1‑150 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  2:427110/1‑150 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  1:903786/1‑150 (MQ=255)
                                                                                                                                     ggTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCAcc                  >  2:3215305/1‑150 (MQ=255)
                                                                                                                                      gTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCa                 >  1:3327281/1‑150 (MQ=255)
                                                                                                                                         tGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCt              >  2:5949149/1‑150 (MQ=255)
                                                                                                                                          gCTGAGATCGCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCtt             <  1:155127/150‑1 (MQ=255)
                                                                                                                                           cTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttg            <  2:1609166/150‑1 (MQ=255)
                                                                                                                                           cTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCCAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCCAGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCAGTCGTTAAACTGCAGCACCAGcaag            >  1:6445504/1‑147 (MQ=38)
                                                                                                                                           cTGAGATCTCGGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttg            <  1:2362804/150‑1 (MQ=255)
                                                                                                                                            tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  1:2715439/1‑150 (MQ=255)
                                                                                                                                            tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  2:4047944/1‑150 (MQ=255)
                                                                                                                                             gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgtt          >  2:2235805/1‑150 (MQ=255)
                                                                                                                                             gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgtt          >  1:201550/1‑150 (MQ=255)
                                                                                                                                              agaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgttg         >  1:3396340/1‑150 (MQ=255)
                                                                                                                                                    tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  2:1804837/1‑150 (MQ=255)
                                                                                                                                                    tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  1:6319845/1‑150 (MQ=255)
                                                                                                                                                    tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCGTGTTGCCGGt    >  1:4637017/1‑150 (MQ=255)
                                                                                                                                                     ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   >  2:4966695/1‑150 (MQ=255)
                                                                                                                                                      tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  <  1:2127188/150‑1 (MQ=255)
                                                                                                                                                      tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  >  1:1838902/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352956‑353254

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