Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 798,667 T→G 8.2% intergenic (+19/+8) fklB‑I → / ← PP_0685 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947798,6670TG 8.2% 641.5 / 32.0 226intergenic (+19/+8)fklB‑I/PP_0685FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (113/91);  new base G (3/17);  total (116/110)
Fisher's exact test for biased strand distribution p-value = 6.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

GCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTTCGTCG  >  NC_002947/798519‑798816
                                                                                                                                                    |                                                                                                                                                     
 cTGATGAACGCTGGCAGCAAAGGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGc                                                                                                                                                     <  1:4080536/150‑1 (MQ=255)
  tgatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGcc                                                                                                                                                    <  1:1305612/150‑1 (MQ=255)
   gatgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCg                                                                                                                                                   >  2:2937764/1‑150 (MQ=255)
    atgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGACGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCgg                                                                                                                                                  >  2:951849/1‑150 (MQ=255)
     tgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                                 <  2:2656592/150‑1 (MQ=255)
     tgaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                                 <  1:2315341/150‑1 (MQ=255)
       aaCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCt                                                                                                                                               <  2:4207488/150‑1 (MQ=255)
        aCGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTc                                                                                                                                              >  1:1009591/1‑150 (MQ=255)
          gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCaa                                                                                                                                            <  2:1729564/150‑1 (MQ=255)
          gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCaa                                                                                                                                            <  1:379510/150‑1 (MQ=255)
          gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                            <  1:1124411/150‑3 (MQ=255)
          gCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                            <  2:3764565/150‑3 (MQ=255)
           cTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTccag                                                                                                                                           <  1:1522566/150‑4 (MQ=255)
            tGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAAtt                                                                                                                                          >  1:1482062/1‑150 (MQ=255)
              gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTcc                                                                                                                                        >  1:4487993/1‑150 (MQ=255)
               cagcaATTGGCGCGTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCATGGCTTTGGCAGCATCCCGCCGCACTGCTTTCTGTTGTTCTACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                       <  1:2582075/150‑1 (MQ=255)
               cagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACTTCGAGCTGCTCGACTTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                       <  1:2993688/150‑1 (MQ=255)
               cagaaATTGGCGCCTGTACGTGCCAAGCGTGCTGGCTTACGGCGCCCAAGGCTTTGGCAGCATCCCGCCGCACCGCGTTCTGGTGTTCGACGTCGCGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                                       <  1:2993557/150‑1 (MQ=255)
                agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAg                                                                                                                                      >  2:3295147/1‑150 (MQ=255)
                 gcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCCATTCCAGt                                                                                                                                     >  2:3657838/1‑150 (MQ=255)
                  caAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTc                                                                                                                                    >  1:65259/1‑150 (MQ=255)
                     aTGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGt                                                                                                                                 >  2:393110/1‑150 (MQ=255)
                        gcgcCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccc                                                                                                                              >  2:1269863/1‑150 (MQ=255)
                          gcCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccccc                                                                                                                            >  2:2472162/1‑150 (MQ=255)
                          gcCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccccc                                                                                                                            >  1:3129587/1‑150 (MQ=255)
                              gTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGGGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCTACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAgg                                                                                                                        <  2:849071/150‑1 (MQ=255)
                                aCGTGCCACGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGgc                                                                                                                      <  2:176845/150‑1 (MQ=255)
                                  gTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGgcgc                                                                                                                    >  2:1393219/1‑150 (MQ=255)
                                    gCCAAGCGTGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTTTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCaa                                                                                                                  <  1:2925627/150‑1 (MQ=255)
                                    gCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACCGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCAGGGCGCAAc                                                                                                                 <  2:3203251/150‑1 (MQ=255)
                                       aaGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCACCGc                                                                                                               >  1:3105308/1‑150 (MQ=255)
                                       aaGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGc                                                                                                               >  1:1406159/1‑150 (MQ=255)
                                       aaGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGc                                                                                                               >  1:603392/1‑150 (MQ=255)
                                        aGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGcc                                                                                                              >  1:1165848/1‑150 (MQ=255)
                                               cggcctccGGCGCCCACGGCGGTGGCAGCATCCCGCCGCACAGCCGACTGGGGTACGACGTCCAGCCGCTCGCCGTTCTGCAAGGCCCTGGCGGCCGCCCCTGCCGGCCCCAATTCCCGTCTGTCCCCCCAGGGCGCAACGCCCGCGCAt                                                                                                       <  2:1226791/143‑1 (MQ=255)
                                                ggCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGTCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATa                                                                                                      <  2:3612363/150‑1 (MQ=255)
                                                ggCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATa                                                                                                      <  2:3535423/150‑1 (MQ=255)
                                                   ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                                   >  1:3051985/1‑150 (MQ=255)
                                                   ttACGGCGCCAACGGCGTGGGCAGAATCCCGCCGCACAGCGTTCTGGGGTTCGCGGTCGAGCTGCACGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGca                                                                                                   <  1:1024667/150‑1 (MQ=255)
                                                       ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCCACGCCCGCGCATAGcagaac                                                                                               >  1:1933976/1‑150 (MQ=255)
                                                          ccccATGGGGTTGGCAGCAGCCCGCCGCACAGGTTTCTGGTTTTCGACGTCGAGCTGCTCGACTTTCTGTAATGCCCTGGAGGCCGCCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaacaga                                                                                            <  1:2643594/149‑1 (MQ=255)
                                                            ccAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaa                                                                                          <  1:1466863/150‑1 (MQ=255)
                                                             cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaac                                                                                         <  2:678232/150‑1 (MQ=255)
                                                             cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaac                                                                                         <  1:4518559/150‑1 (MQ=255)
                                                              aaGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaaca                                                                                        >  1:3836388/1‑150 (MQ=255)
                                                               aGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaacag                                                                                       >  1:1666255/1‑150 (MQ=255)
                                                                ggCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaacaga                                                                                      <  2:1673991/150‑1 (MQ=255)
                                                                   gTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTc                                                                                   <  2:1365680/150‑1 (MQ=255)
                                                                   gTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTc                                                                                   <  2:227963/150‑1 (MQ=255)
                                                                       gcagcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTa                                                                               >  2:1564134/1‑150 (MQ=255)
                                                                       gcagcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTa                                                                               <  2:3469872/150‑1 (MQ=255)
                                                                       gcagcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTa                                                                               <  2:640937/150‑1 (MQ=255)
                                                                       gcagcaCCCCGCCGCACAGGGTTCTGGTGTTCGCCGTCGCGGTGCTCGACGTTCTGAAAGCCCCTGGAGGCCGGCCCTGCCGGCCACAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTa                                                                               <  1:3795166/150‑1 (MQ=255)
                                                                        cagctccccGCCGCACCGCGTTGTGGTGTTCGACGTCGCGGGGCTCGACGTTGTGTAATGCCGTGGCGGCCGGCCCGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAc                                                                              <  1:2712024/144‑1 (MQ=255)
                                                                        cagcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATTCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAc                                                                              >  2:3276043/1‑150 (MQ=255)
                                                                            aTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGt                                                                          >  2:611695/1‑150 (MQ=255)
                                                                             tCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACACATTCCGTACCAGtt                                                                         >  1:3178658/1‑150 (MQ=255)
                                                                               ccgccgCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTcc                                                                       >  2:1839241/1‑150 (MQ=255)
                                                                                cgccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCt                                                                      >  2:4286548/1‑150 (MQ=255)
                                                                                cgccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCt                                                                      >  2:3850281/1‑150 (MQ=255)
                                                                                cgccgcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCt                                                                      >  2:3424676/1‑150 (MQ=255)
                                                                                 gccgcACCGCTTTGTGGTGTGCGCCGTCGCGCGGCGCGCCGTTCTGTACGGCCCGGGCGGCCGGCCCGGCCGGCCGCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCtt                                                                     <  1:898281/150‑1 (MQ=255)
                                                                                    gcACAGCGTTCTGGTGTTCGACGTCGCGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAg                                                                  <  2:3304991/150‑1 (MQ=255)
                                                                                      acaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCATTGAGCa                                                                >  1:3226588/1‑150 (MQ=255)
                                                                                      acaGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTACGTACCAGTTCCTTGAGCa                                                                >  1:3114293/1‑150 (MQ=255)
                                                                                         gCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGAGTCCGTCCCAGTTCCTTGAGCAccc                                                             >  2:2576765/1‑150 (MQ=255)
                                                                                            ttCTGGTGTTCGTCGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACTCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCg                                                          <  2:3366082/150‑1 (MQ=255)
                                                                                            ttCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCg                                                          <  2:2058562/150‑1 (MQ=255)
                                                                                                   gtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACt                                                   <  2:4121728/150‑1 (MQ=255)
                                                                                                    ttCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTc                                                  <  1:2946146/150‑1 (MQ=255)
                                                                                                      gaaCGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGa                                                >  1:1240540/3‑150 (MQ=255)
                                                                                                      gaaCGTCGAGCAGCTCGACGTCCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGa                                                >  1:4269932/3‑150 (MQ=255)
                                                                                                      ccaCGCCGTGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGa                                                <  1:1620244/148‑1 (MQ=255)
                                                                                                       gACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                               <  2:2543091/150‑1 (MQ=255)
                                                                                                       gACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                               >  1:3424580/1‑150 (MQ=255)
                                                                                                       aaCGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGaa                                               >  2:1984842/2‑150 (MQ=255)
                                                                                                        aCGTCGAGCAGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAAt                                              >  1:3574338/1‑150 (MQ=255)
                                                                                                         cGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAAtt                                             <  2:1460059/150‑1 (MQ=255)
                                                                                                           tCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCa                                           >  1:961025/1‑150 (MQ=255)
                                                                                                              agcagcTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGcc                                        >  1:4472271/1‑150 (MQ=255)
                                                                                                               gctgctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGccc                                       >  1:2547444/1‑150 (MQ=255)
                                                                                                                 tgctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGc                                     <  1:3711518/150‑1 (MQ=255)
                                                                                                                   ctCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTc                                   >  2:2784159/1‑150 (MQ=255)
                                                                                                                     ctaCGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCACAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCAt                                 <  2:130201/148‑1 (MQ=255)
                                                                                                                     cGACGTTATGTAATCCCCTGGAGGCCGGCCCTGCCGGCATCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCAt                                 <  1:4428838/150‑1 (MQ=255)
                                                                                                                           tCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                           <  2:3327657/150‑1 (MQ=255)
                                                                                                                             tGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGt                         <  2:1962976/150‑1 (MQ=255)
                                                                                                                                   gCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGGTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCt                   >  2:4427774/1‑150 (MQ=255)
                                                                                                                                     cccggACGCCGGCCCCCCCGGCCACAATACCAGTCTGTCCCCCCACGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCCgg                 <  1:2048006/147‑1 (MQ=255)
                                                                                                                                     aattgAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTgg                 >  1:2483112/5‑150 (MQ=255)
                                                                                                                                      cTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGt                >  1:1973515/1‑150 (MQ=255)
                                                                                                                                      attgAGGCCGGCAGGGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGc                                                                                                    <  2:3128620/63‑1 (MQ=37)
                                                                                                                                       ttgAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTc               >  1:3313237/3‑150 (MQ=255)
                                                                                                                                        ggAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCa              >  2:3446573/1‑150 (MQ=255)
                                                                                                                                            gccggccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCa          >  2:964882/1‑150 (MQ=255)
                                                                                                                                            gccggccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCa          >  2:1049618/1‑150 (MQ=255)
                                                                                                                                             ccggccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAg         <  2:1009591/150‑1 (MQ=255)
                                                                                                                                              cggccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGt        >  1:4180087/1‑150 (MQ=255)
                                                                                                                                                  ccTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTtcgt    >  2:1621968/1‑150 (MQ=255)
                                                                                                                                                    tGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTtcgtcg  >  2:732077/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTGATGAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACGCAGTTCGTCG  >  NC_002947/798519‑798816

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: