Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,070 G→A 8.2% intergenic (‑111/+17) PP_0386 ← / ← rpoD sensory box protein/RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,0700GA 8.2% 267.3 / 10.1 97intergenic (‑111/+17)PP_0386/rpoDsensory box protein/RNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (49/38);  new base A (6/2);  total (56/41)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTG  >  NC_002947/470922‑471216
                                                                                                                                                    |                                                                                                                                                  
aTCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGTCCTgg                                                                                                                                                   >  1:483480/1‑150 (MQ=255)
     gATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgg                                                                                                                                              <  2:483480/150‑1 (MQ=255)
     gATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgg                                                                                                                                              <  2:1190124/150‑1 (MQ=255)
           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAATGAAAAACCCCGGCCTGGGc                                                                                                                                                 >  1:414451/1‑141 (MQ=255)
           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGc                                                                                                                                                 <  2:414451/141‑1 (MQ=255)
             aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGt                                                                                                                                      >  1:1022155/1‑150 (MQ=255)
             aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGt                                                                                                                                      <  2:1281693/150‑1 (MQ=255)
              cAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtc                                                                                                                                     <  1:45753/150‑1 (MQ=255)
               aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtc                                                                                                                                     <  1:1034738/149‑1 (MQ=255)
               aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtc                                                                                                                                     >  2:1034738/1‑149 (MQ=255)
                gAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGATGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtcat                                                                                                                                   >  2:1207643/1‑150 (MQ=255)
                           ttCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                        <  2:320732/150‑1 (MQ=255)
                                  aCATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAg                                                                                                                 <  2:48832/150‑1 (MQ=255)
                                      ggCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAATTGTGGCTTCATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCa                                                                                                             <  1:1349768/150‑1 (MQ=255)
                                         aGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                                                                                                                                        >  2:959983/1‑120 (MQ=255)
                                         aGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                                                                                                                                        <  1:959983/120‑1 (MQ=255)
                                          gCCTCGCACGCAGGTGGTCTGAACACGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCGCTAAAGAAAAGCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTg                                                                                                         <  2:1197834/150‑1 (MQ=255)
                                            cTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCGGGGtt                                                                                                                                          >  2:307857/1‑115 (MQ=255)
                                            cTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCGGGGtt                                                                                                                                          <  1:307857/115‑1 (MQ=255)
                                             tCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTc                                                                                                      >  2:8628/1‑150 (MQ=255)
                                             tCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCAGATGGGGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTc                                                                                                      >  1:1412608/1‑150 (MQ=255)
                                                 aCTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                                 >  1:1359158/1‑150 (MQ=255)
                                                   tCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCg                                                                                                >  2:1486745/1‑150 (MQ=255)
                                                       gtggtCTGAAGTTGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTc                                                                                            <  1:1486745/150‑1 (MQ=255)
                                                       gtggtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTCAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACtc                                                                                                                              <  1:1572445/116‑1 (MQ=255)
                                                       gtggtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                                                                                                                                             <  2:618935/101‑1 (MQ=255)
                                                       gtggtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                                                                                                                                             >  1:618935/1‑101 (MQ=255)
                                                       gtggtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACtc                                                                                                                              >  2:1572445/1‑116 (MQ=255)
                                                             tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgggg                                                                                                                                            <  2:287175/96‑1 (MQ=255)
                                                             tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgggg                                                                                                                                            >  1:287175/1‑96 (MQ=255)
                                                             tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgggggg                                                                                      >  1:3598/1‑150 (MQ=255)
                                                              gAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGgc                                                                                     >  2:615903/1‑150 (MQ=255)
                                                               aaGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCa                                                                                                             >  2:918036/1‑125 (MQ=255)
                                                               aaGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCa                                                                                                             <  1:918036/125‑1 (MQ=255)
                                                                aGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTAGCTGCGAGTCGGGTGgcgc                                                                                   <  1:1207643/150‑1 (MQ=255)
                                                                   cGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGc                                                                                <  1:933037/150‑1 (MQ=255)
                                                                    gAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGTGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCt                                                                               >  2:660051/1‑150 (MQ=255)
                                                                    gAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCt                                                                               >  1:1074448/1‑150 (MQ=255)
                                                                    gAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCt                                                                               >  1:416011/1‑150 (MQ=255)
                                                                       tgtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgc                                                                            >  1:748944/1‑150 (MQ=255)
                                                                          aGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGc                                                                                                        >  1:716735/1‑119 (MQ=255)
                                                                          aGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGc                                                                                                        <  2:716735/119‑1 (MQ=255)
                                                                                gCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaag                                                                                                                    <  2:1495808/101‑1 (MQ=255)
                                                                                gCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaag                                                                                                                    >  1:1495808/1‑101 (MQ=255)
                                                                                 cATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCctt                                                                  <  1:8628/150‑1 (MQ=255)
                                                                                  aTTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttc                                                                 <  1:1244354/150‑1 (MQ=255)
                                                                                   ttGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcg                                                                >  1:881040/1‑150 (MQ=255)
                                                                                   ttGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcg                                                                >  2:1534960/1‑150 (MQ=255)
                                                                                      ccGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgctt                                                             >  2:885316/1‑150 (MQ=255)
                                                                                        gCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGc                                                                                                        <  1:1345568/105‑1 (MQ=255)
                                                                                        gCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGc                                                                                                        >  2:1345568/1‑105 (MQ=255)
                                                                                               tgCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGAc                                                    <  2:1022155/150‑1 (MQ=255)
                                                                                                  aaCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGa                                                 >  2:561955/1‑150 (MQ=255)
                                                                                                       aaGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCt                                           <  2:1359158/150‑1 (MQ=255)
                                                                                                          gCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCa                                         >  2:1138877/1‑150 (MQ=255)
                                                                                                           cGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAc                                        >  2:1566491/1‑150 (MQ=255)
                                                                                                              gCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgcg                                     >  1:954856/1‑150 (MQ=255)
                                                                                                              gCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCAGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgcg                                     >  1:25613/1‑150 (MQ=255)
                                                                                                               cTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                          <  1:932019/96‑1 (MQ=255)
                                                                                                               cTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                          >  2:932019/1‑96 (MQ=255)
                                                                                                                      tCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGa                                                          >  2:209675/1‑121 (MQ=255)
                                                                                                                      tCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGa                                                          <  1:209675/121‑1 (MQ=255)
                                                                                                                      tCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCa                             >  2:101827/1‑150 (MQ=255)
                                                                                                                      tCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCa                             >  1:11971/1‑150 (MQ=255)
                                                                                                                        aCGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGGGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                           >  1:552643/1‑150 (MQ=255)
                                                                                                                        aCGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                           >  2:1148706/1‑150 (MQ=255)
                                                                                                                        aCGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                           <  1:885316/150‑1 (MQ=255)
                                                                                                                            gCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGACGAGGCAGGAACGCAGGGGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTCCGCTTCGATCTGACGGATCCGCTCGCGCGTTACGTCAACCTga                       >  2:648719/1‑149 (MQ=255)
                                                                                                                             ccACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCt                      >  1:39105/1‑150 (MQ=255)
                                                                                                                              cACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCtt                     <  2:1074448/150‑1 (MQ=255)
                                                                                                                                atatgacaaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                   >  2:1229265/8‑150 (MQ=255)
                                                                                                                                 tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                  <  2:748944/150‑1 (MQ=255)
                                                                                                                                 tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCAAACTGAATGCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                  >  2:457969/1‑150 (MQ=255)
                                                                                                                                   aaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCa                >  2:1055892/1‑150 (MQ=255)
                                                                                                                                    aGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGGGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAc               <  2:552643/150‑1 (MQ=255)
                                                                                                                                    aGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAc               <  2:416011/150‑1 (MQ=255)
                                                                                                                                      caaaaCCCCGGCCCAGGCCGGGGTTTTGTCCTCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACct             >  1:679253/2‑150 (MQ=255)
                                                                                                                                        aaaCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACctct           >  2:1490006/1‑150 (MQ=255)
                                                                                                                                        aaaCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACctct           >  1:1245730/1‑150 (MQ=255)
                                                                                                                                        aaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCTCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACctct           >  2:1208514/1‑150 (MQ=255)
                                                                                                                                          aCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc         <  1:561955/150‑1 (MQ=255)
                                                                                                                                           ccccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgcttc                                                            <  1:1016281/98‑1 (MQ=255)
                                                                                                                                           ccccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgcttc                                                            >  2:1016281/1‑98 (MQ=255)
                                                                                                                                           ccccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCg        >  2:478744/1‑150 (MQ=255)
                                                                                                                                           ccccGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgc                                                                            <  2:1625351/82‑1 (MQ=255)
                                                                                                                                           ccccGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgc                                                                            >  1:1625351/1‑82 (MQ=255)
                                                                                                                                            cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGTAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa       >  2:658033/1‑150 (MQ=255)
                                                                                                                                            cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa       >  2:1053005/1‑150 (MQ=255)
                                                                                                                                            cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa       >  1:718948/1‑150 (MQ=255)
                                                                                                                                            cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa       <  2:881040/150‑1 (MQ=255)
                                                                                                                                             ccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAg      >  2:127944/1‑150 (MQ=255)
                                                                                                                                             ccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAg      >  1:96394/1‑150 (MQ=255)
                                                                                                                                              cGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg     <  1:1566491/150‑1 (MQ=255)
                                                                                                                                              cGGCCTGGGACGTGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg     <  1:660051/150‑1 (MQ=255)
                                                                                                                                                gCCTGGGCCGGGGTTTTGTCATCACTCGGCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgtgt   <  1:325997/150‑1 (MQ=255)
                                                                                                                                                 ccTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGgtg  >  1:1224162/1‑150 (MQ=255)
                                                                                                                                                   tGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAc               >  1:830245/1‑135 (MQ=255)
                                                                                                                                                   tGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAc               <  2:830245/135‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTG  >  NC_002947/470922‑471216

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