Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,967 C→G 7.2% intergenic (+315/+164) PP_4318 → / ← PP_4319 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,909,9670CG 7.2% 230.4 / 6.4 82intergenic (+315/+164)PP_4318/PP_4319transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/43);  new base G (3/3);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCA  >  NC_002947/4909819‑4910115
                                                                                                                                                      |                                                                                                                                                     
cTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatcca                                                                                                                                                      >  2:2359596/1‑150 (MQ=37)
 tttCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatccat                                                                                                                                                     <  2:1731779/150‑1 (MQ=37)
    ccGCCATTACGGCGGGTAACTTTTTGTCAAACGCGACAAAAAAAAAAAAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATggg                                                                                                                                                  <  2:1131525/150‑1 (MQ=11)
     cGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAAGGCCCCATCCATggag                                                                                                                                                 <  2:1656144/150‑3 (MQ=18)
     cGCCCTTACGGCGGGTCACTTTTTGTCAAACCAGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCACAGATGGTGATGAT‑‑GCGGCTGGCGCGCTACGGTCTTGCTCCCGGGAGTAAGCGCTGCAGGCCCCATCCATgggg                                                                                                                                                 >  1:2227505/1‑150 (MQ=2)
      gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAGGCCCCATGGATggg                                                                                                                                                  <  1:1088056/147‑1 (MQ=12)
      gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAGGCCCCATGGATggg                                                                                                                                                  >  2:1088056/1‑147 (MQ=255)
               gcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCt                                                                                                                                      >  1:3064975/1‑150 (MQ=32)
                 gggTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTg                                                                                                                                    <  2:149396/150‑1 (MQ=32)
                   gTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGAc                                                                                                                                  <  1:2359596/150‑1 (MQ=32)
                   gTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACCCTCCGGGGTTCCATCGCAACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGTTCCTCAGCGCTTGAc                                                                                                                                  <  2:1513605/150‑1 (MQ=11)
                   gTCACTTTTTGTCAAACGCAACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGAc                                                                                                                                  <  1:2033806/150‑1 (MQ=17)
                    tCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATGGATGGGGCCt                                                                                                                                              >  2:2442998/1‑138 (MQ=16)
                    tCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATGGATGGGGCCt                                                                                                                                              <  1:2442998/138‑1 (MQ=16)
                    tCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACg                                                                                                                                 <  1:143573/150‑1 (MQ=32)
                    tCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACg                                                                                                                                 >  1:515615/1‑150 (MQ=32)
                     cACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCATATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                <  2:2299281/150‑1 (MQ=17)
                      aCTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACggg                                                                                                                               <  1:182886/150‑1 (MQ=32)
                       cTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTACAGGCCCCATCCATGGGGc                                                                                                                                                >  2:2267711/1‑133 (MQ=18)
                       cTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTACAGGCCCCATCCATGGGGc                                                                                                                                                <  1:2267711/133‑1 (MQ=18)
                        tttttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCa                                                                                                                             >  2:3055640/1‑150 (MQ=32)
                          tttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATc                                                                                                                           <  2:692494/150‑1 (MQ=32)
                              tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATg                                                                                                                       <  2:1749661/150‑1 (MQ=32)
                              tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATCGGGCCTCAGCGCTTGACGGGCATCCATg                                                                                                                       >  1:2032397/1‑150 (MQ=17)
                                 aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                   >  2:2935340/1‑120 (MQ=33)
                                 aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                   <  1:2935340/120‑1 (MQ=33)
                                 aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                    >  2:2696452/1‑150 (MQ=32)
                                   cgcgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                  <  1:2696452/150‑1 (MQ=32)
                                   cgcgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                  <  1:1235660/150‑1 (MQ=32)
                                    gcgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                 >  1:479517/1‑150 (MQ=32)
                                      gACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGGTCAGGGGCGAACTCGCTGCATGCTCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                                               >  1:1136051/1‑150 (MQ=11)
                                         aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                   <  2:2888423/112‑1 (MQ=32)
                                         aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgg                                                                                                                                                   >  1:2888423/1‑112 (MQ=33)
                                              gTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                       >  2:2836286/1‑150 (MQ=34)
                                               tAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                   <  2:972100/121‑1 (MQ=33)
                                               tAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                   >  1:972100/1‑121 (MQ=33)
                                                      aaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCa                                                                                                                             <  1:3070387/120‑1 (MQ=33)
                                                      aaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCa                                                                                                                             >  2:3070387/1‑120 (MQ=33)
                                                             tGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                        >  2:1673514/1‑150 (MQ=38)
                                                              ggCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                       <  2:1371271/150‑1 (MQ=38)
                                                              ggCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                       <  2:2533032/150‑1 (MQ=39)
                                                                  ccTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCt                                                                                   <  2:479517/150‑1 (MQ=38)
                                                                   cTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCtt                                                                                  <  2:3064975/150‑1 (MQ=38)
                                                                    tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCCTGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTc                                                                                 <  1:2917538/150‑1 (MQ=37)
                                                                     atcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCt                                                                                <  1:2192954/150‑1 (MQ=38)
                                                                      tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTg                                                                               >  1:1184372/1‑150 (MQ=39)
                                                                       catcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGa                                                                              <  2:2982683/150‑1 (MQ=39)
                                                                       catcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGa                                                                              <  1:1830655/150‑1 (MQ=39)
                                                                        atcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGaa                                                                             >  1:2052739/1‑150 (MQ=38)
                                                                          ccGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGt                                                                           <  2:2032397/150‑1 (MQ=38)
                                                                          ccGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGt                                                                           <  2:1184372/150‑1 (MQ=39)
                                                                           cGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTc                                                                          <  2:211866/150‑1 (MQ=39)
                                                                            ggCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCg                                                                         >  2:2500879/1‑150 (MQ=37)
                                                                               cccTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCCGCTCCCTTCGCAAGCCCTGCGCCCGGCCTAGTGAAGTCGaaa                                                                      >  2:2096673/1‑150 (MQ=11)
                                                                                 cTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATc                                                                    >  1:1030240/1‑150 (MQ=39)
                                                                                  tGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCt                                                                   <  1:3055640/150‑1 (MQ=39)
                                                                                             cgcaggggTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCt                                                          <  1:2445551/146‑1 (MQ=255)
                                                                                                cGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                       <  2:515615/150‑1 (MQ=38)
                                                                                                        ccTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCATTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcg                                               >  1:530338/1‑150 (MQ=255)
                                                                                                          tCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcgcg                                             >  2:1921884/1‑150 (MQ=40)
                                                                                                           cGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCgcgcgc                                            <  2:1030240/150‑1 (MQ=255)
                                                                                                               aCGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAc                                        <  2:2052739/150‑1 (MQ=255)
                                                                                                                 ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      <  1:1673514/150‑1 (MQ=255)
                                                                                                                 ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGTGTCGCCTGAACTACCGCGCGCTTACaa                                      <  2:563335/150‑1 (MQ=255)
                                                                                                                         tCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagc                              >  2:1931351/1‑150 (MQ=255)
                                                                                                                         tCCCGGGAGTACGCCCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCTGTCCCCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagc                              >  1:1235221/1‑150 (MQ=255)
                                                                                                                           ctgggAGTACACGCTGCAGGCCGCATCCATGGGGCCTTAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcag                            >  2:198611/3‑150 (MQ=255)
                                                                                                                             gggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaa                          <  1:2836286/150‑1 (MQ=255)
                                                                                                                             gggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaa                          <  1:1931351/150‑1 (MQ=255)
                                                                                                                                     gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCGCCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCTCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc                  <  2:2227505/150‑1 (MQ=255)
                                                                                                                                     gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCACGCCCTTCTCAAGTCGCAAGCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc                  >  1:2665242/1‑150 (MQ=255)
                                                                                                                                     gcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc                  <  1:2185665/150‑1 (MQ=17)
                                                                                                                                         tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAgcagaagcagaagcagaagcagaa              >  1:1942761/1‑150 (MQ=255)
                                                                                                                                          tgaGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      >  2:2351455/3‑125 (MQ=255)
                                                                                                                                          tgaGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      <  1:2351455/123‑1 (MQ=255)
                                                                                                                                          gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaag             >  2:733668/1‑150 (MQ=255)
                                                                                                                                            aGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagca           >  2:511322/1‑150 (MQ=255)
                                                                                                                                             ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag          >  2:616311/1‑150 (MQ=255)
                                                                                                                                               ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAGAAGCAGAAGCAGAagcagaagcagca        <  2:1942761/150‑1 (MQ=255)
                                                                                                                                               ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc                  >  1:1611642/1‑140 (MQ=255)
                                                                                                                                               ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc                  <  2:1611642/140‑1 (MQ=255)
                                                                                                                                                cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGc                                                                                     >  2:958020/1‑72 (MQ=35)
                                                                                                                                                cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGc                                                                                     <  1:958020/72‑1 (MQ=35)
                                                                                                                                                cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcag       >  1:2073663/1‑150 (MQ=255)
                                                                                                                                                   atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcagaag    >  1:2167843/1‑150 (MQ=255)
                                                                                                                                                     ccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAgcagaagcagcagaagca  <  2:2073663/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCA  >  NC_002947/4909819‑4910115

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