Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,105 T→G 9.1% 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,1050TG 9.1% 612.1 / 3.7 256intergenic (+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 (117/104);  new base G (14/18);  total (133/123)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTT  >  NC_002947/352956‑353254
                                                                                                                                                     |                                                                                                                                                     
gcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATctct                                                                                                                                                       <  1:4683490/150‑1 (MQ=255)
  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                                                                     <  2:6842872/150‑1 (MQ=255)
  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCCCCCGGGTTGCGGAGATCTCtgt                                                                                                                                                     >  1:7463741/1‑150 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  2:4731494/150‑1 (MQ=255)
    gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATGGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                                                                   <  2:5371627/150‑1 (MQ=255)
      aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGa                                                                                                                                                 <  1:3449320/150‑1 (MQ=255)
      aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCAGTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATATCTGTGGGa                                                                                                                                                 >  2:1333032/1‑150 (MQ=255)
         ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCg                                                                                                                                              >  2:779515/1‑150 (MQ=255)
          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                                             <  1:2385754/150‑1 (MQ=255)
          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCACTGTGGGAGCgg                                                                                                                                             >  2:4557273/1‑150 (MQ=255)
          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCGGCCCCAGGTTTGCGGAGATGACTGGGGGAGCgg                                                                                                                                             >  1:2976111/1‑150 (MQ=255)
            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                           >  2:6733394/1‑150 (MQ=255)
            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                                                           <  2:1557215/150‑1 (MQ=255)
              cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                                                         <  2:1595936/150‑1 (MQ=255)
               tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCGCTGTGGGAGCGGGCGTg                                                                                                                                        >  2:7314240/1‑150 (MQ=255)
                ttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                                                                       <  1:3885940/150‑1 (MQ=255)
                  ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGGGTGCTGAGATCTCTGTGGGGGCGGGCGTGccc                                                                                                                                     >  1:6031231/1‑150 (MQ=255)
                    cGAGTACACAGTGCCGGAGGCCAAGGCCTACAGGGCTTCGCGCGGGATCGTCGTTCGCTGATCCAGTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                                                                   <  2:6103583/150‑1 (MQ=255)
                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCGGAGATCTCTGGGGGAGCGGGCGGGCCcgcgggccgg                                                                                                                            >  1:1214951/1‑144 (MQ=38)
                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGGTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                                                                    >  2:2679250/1‑150 (MQ=255)
                                    gcggCCACGGCCTCCCTGGCTTCGCGCGGGATCGTCGTTCGCGGATCCATTGCCCGGGGCGGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCACCGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCg                                                                                                                   <  1:729344/148‑1 (MQ=38)
                                        ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                                                               >  2:2566648/1‑150 (MQ=255)
                                        ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCACCGGTTTTCTGAGATCTATGTTGGAGCGGGCGTGCCCGCGCAGGGCCGCCAAGTGGccc                                                                                                               >  1:7427053/1‑150 (MQ=255)
                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcccc                                                                                                              <  2:7252728/150‑1 (MQ=255)
                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                                                                 >  2:7171552/1‑147 (MQ=255)
                                           aGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATGGACCGGGGCGGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTTCGCGGcccttc                                                                                                            <  2:302406/150‑4 (MQ=38)
                                                tACTTGGCTTCGCGCGGTATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                                                                       <  1:7854073/150‑1 (MQ=255)
                                                tACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCaagttct                                                                                                       >  1:2697598/1‑145 (MQ=255)
                                                tACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCGGCTTTGCAGCCCTTCGGGGGCACGCCCGCTCCCCCAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCGCCAATTGCt                                                                                                       <  2:8288735/150‑1 (MQ=38)
                                                 acctGGGTTCGCGGGGGAGCGTCGTTCGCTGAGCCATTGCCCGGGGGTGCTTGGCCGCCCTGCGGGGGCCCGCCCGCTCCCCCTGGGTTGCTGCGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtt                                                                                                      <  1:1922354/147‑1 (MQ=25)
                                                   aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGGCCGCTCCCTCGGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGTAGCGGCCCCAATTGCTTCa                                                                                                    >  2:7501443/1‑150 (MQ=39)
                                                   aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCa                                                                                                    >  2:479284/1‑150 (MQ=255)
                                                     ggCTTCGCGCGGGATCGTCGGTCGCTGATCCATTGACCGGGCCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCCCCGGGTGGCGGAGCTCGCGGGGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCCAGTGCTTCAAc                                                                                                  >  1:8773848/1‑150 (MQ=255)
                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                               <  2:317189/150‑1 (MQ=255)
                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCGCCCACACGGTTGGGGAGATCTCTGTGGGAGCGGGCGTGCCCCCCGAGGGCCCCGCACCGGCCCCACATGCGTCCAccct                                                                                               >  1:149882/1‑148 (MQ=25)
                                                         tCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                                                              <  2:1079380/150‑1 (MQ=255)
                                                             gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGTTTGCTGAGATATCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCCCGAAGCGGCCCCAATTGCTTCAAACGGGGTTc                                                                                          >  2:1704702/1‑150 (MQ=38)
                                                              cgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCCACAGGGTTGGGGGGCTCTCGGGGGGAGCGGGCGTGCCCGCCGAGGGCCGCGCAGCGGCCCCCATTTCTTCAACCGTGGttgt                                                                                         >  1:4854742/1‑148 (MQ=25)
                                                                ggATCGTCGTTCGCTGATCCATTGACCGGGGCGGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                                                                       <  2:5130881/150‑1 (MQ=255)
                                                                 gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                                                                      <  2:1883236/150‑1 (MQ=255)
                                                                    cgtcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCACTTGCTTCAACCGTGGTTCTTGCCGa                                                                                   <  2:4996871/150‑1 (MQ=255)
                                                                      tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagc                                                                                 >  2:7011453/1‑150 (MQ=255)
                                                                            gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                                                                           >  2:1038613/1‑150 (MQ=255)
                                                                            gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                                                                           >  2:2279194/1‑150 (MQ=255)
                                                                             cTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGATCTTTCCGAGCAgcgcct                                                                          >  2:7049598/1‑148 (MQ=255)
                                                                              tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                         >  2:5718879/1‑150 (MQ=255)
                                                                              tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                         <  2:971698/150‑1 (MQ=255)
                                                                              tGATCAATTGACGGGGGGTGCTTGGCCGCCAGTCGCGGGCACGCCCGCCCCCACTGGTTTGCTGAGATCTGTGGGGGCGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                                                         <  1:1120848/150‑1 (MQ=255)
                                                                                  ccATTGACCGGGGCTGCTTTGCCGCCCTTCGCGGGCACGCCCGCACCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaa                                                                     <  1:4156237/150‑1 (MQ=255)
                                                                                  ccATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGGTTGCGGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaa                                                                     >  1:2033939/1‑150 (MQ=38)
                                                                                   cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                                                                    <  1:6973898/150‑1 (MQ=255)
                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCAGAGGTTTGCTGGGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGGAGCGGCCCCCCTTGCTTCAACCGTGGTTGTTGGCGCGCCGCGCGTGGTAAAGCTCg                                                               >  1:5665209/1‑150 (MQ=37)
                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                                                               <  2:4179130/150‑1 (MQ=255)
                                                                                         cgggggCTGCTGTGCCGCCCGGCGCGGGAACGCCCGCTCCCACCGGTTGGCTGATACCTGTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCgg                                                              <  1:1701287/148‑1 (MQ=37)
                                                                                         ccGGGGCTGCTTTGCCGCCCTGGGCGGCCCCCCCCGCGCCCACCGGTTGGCTGAGACCTCGGTGGGTGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCgg                                                              <  1:222213/150‑1 (MQ=37)
                                                                                          cGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTCCCACAGGGTTGCGGAGATCGCTGTGGGAGCGGGGGGGCCCGCGAAGGGCCGCGGATCGGCCCCCATTGCTGCAACCCTGGGTGGTGCCGAGCAGCGCCTGGTAAAGCTCggg                                                             >  2:2410346/1‑149 (MQ=25)
                                                                                           agggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                                                              >  1:7157611/2‑116 (MQ=255)
                                                                                            gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTAt                                                           >  2:7103290/1‑150 (MQ=255)
                                                                                               ctgctTTGCCGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCGGCGATCTCGGGGGGTGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc                                                        <  1:375956/150‑1 (MQ=38)
                                                                                               ctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc                                                        <  2:3768431/150‑1 (MQ=255)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  1:5824117/150‑1 (MQ=255)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTGCGGCGATCTCGGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCCGCCGCCCCAATTGCGTCAACCGTGGTTCATGCCGAGCAGCGCGTGGTAAAGCGCGGTATCCCCgc                                                     >  2:3414302/1‑149 (MQ=37)
                                                                                                  ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                     <  2:6442428/150‑1 (MQ=255)
                                                                                                       cAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATc                                                <  2:7406301/150‑1 (MQ=255)
                                                                                                           ccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca                                            <  2:5393416/150‑1 (MQ=255)
                                                                                                               cgcgGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCcaccac                                        <  2:6544544/150‑1 (MQ=255)
                                                                                                                gcgGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCaccacc                                       >  2:4283862/1‑150 (MQ=255)
                                                                                                                     cACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTc                                  <  2:851447/150‑1 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGtt                               >  1:6319393/1‑150 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGtt                               <  2:31117/150‑1 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGtt                               >  1:4942625/1‑150 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTAt                             >  2:4656679/1‑150 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTAt                             >  2:4426978/1‑150 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTAt                             >  2:6320113/1‑150 (MQ=255)
                                                                                                                           cGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATc                            <  2:3197622/150‑1 (MQ=255)
                                                                                                                            gCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATcc                           >  1:6367608/1‑150 (MQ=255)
                                                                                                                             cTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCt                          <  1:4283862/150‑1 (MQ=255)
                                                                                                                               cccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgc                        <  2:4053727/150‑1 (MQ=255)
                                                                                                                                ccACCGGTTTGCTTATATCTCTGTGGGTGCGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTTGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCATgca                       <  1:1333032/150‑1 (MQ=39)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       <  1:6533673/150‑1 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  1:5281207/1‑150 (MQ=255)
                                                                                                                                ccACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgca                       >  1:928872/1‑150 (MQ=255)
                                                                                                                                 cacaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcag                      >  2:8339316/1‑150 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  1:8199926/1‑150 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  2:630560/1‑150 (MQ=255)
                                                                                                                                  acaGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagc                     >  1:5587901/1‑150 (MQ=255)
                                                                                                                                   caGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTgcagca                    >  2:5880352/1‑150 (MQ=255)
                                                                                                                                     ggTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCAcc                  <  2:1061114/150‑1 (MQ=255)
                                                                                                                                     ggTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCAcc                  <  2:5032277/150‑1 (MQ=255)
                                                                                                                                            tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  1:4998703/1‑150 (MQ=255)
                                                                                                                                            tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  1:2797965/1‑150 (MQ=255)
                                                                                                                                            tGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgt           >  1:8314332/1‑150 (MQ=255)
                                                                                                                                             gagaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgtt          >  2:8611460/1‑150 (MQ=255)
                                                                                                                                              agaTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCttgttg         >  2:2009563/1‑150 (MQ=255)
                                                                                                                                                aTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGcc       >  2:2897152/1‑150 (MQ=255)
                                                                                                                                                aTCTCTGTGGGAGCAGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGcc       >  1:1078124/1‑150 (MQ=255)
                                                                                                                                                 tctcTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCg      >  2:3700999/1‑150 (MQ=255)
                                                                                                                                                   tctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGt    >  2:3892186/1‑150 (MQ=255)
                                                                                                                                                    ctGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGtt   <  2:6031408/150‑1 (MQ=255)
                                                                                                                                                     tgtgGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGttt  <  2:5618546/150‑1 (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: