Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,066 T→C 10.7% intergenic (‑110/+18) 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,0660TC10.7% 270.9 / 19.8 104intergenic (‑110/+18)PP_0386/rpoDsensory box protein/RNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base T (50/42);  new base C (8/3);  total (59/45)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

Only 100 of 104 total aligned reads displayed.

TCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGT  >  NC_002947/470920‑471212
                                                                                                                                                  |                                                                                                                                                    
tCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑gg                                                                                                                                                  <  3:360476/149‑1 (MQ=255)
   ggATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGcc                                                                                                                                               <  5:357667/149‑1 (MQ=255)
    gATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCg                                                                                                                                              >  5:287708/1‑149 (MQ=255)
        aGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGTCCT‑‑GGGCCggg                                                                                                                                            <  5:325880/147‑1 (MQ=255)
          cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGttt                                                                                                                                        >  3:334688/1‑149 (MQ=255)
           aaCAGAATGTCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGTTGTGCAACTGAAGGCGGGCTTTATGTCACTGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGtttt                                                                                                                                       <  1:38995/149‑1 (MQ=255)
           aaCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGtttt                                                                                                                                       <  2:187217/149‑1 (MQ=255)
            actgAATGGCCGGCTTCAAGGAGATGGCAGCGTCGCACTCAGGTGGTCTGAAGTCGTGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGtttt                                                                                                                                       <  2:178630/145‑1 (MQ=255)
            aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTTGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGCGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTg                                                                                                                                      <  3:158800/149‑1 (MQ=255)
            aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTg                                                                                                                                      >  8:360000/1‑149 (MQ=255)
            aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTCCCGCAGGTGTGCAACGGAAGGCGGGCTTTATCTCACGTGCCACTAAAGAAAAACCCCGGCCT‑‑GGCCCGGGGtttt                                                                                                                                       >  1:178630/1‑148 (MQ=255)
              aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGt                                                                                                                                          <  8:19160/143‑1 (MQ=255)
              aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGt                                                                                                                                          >  7:19160/1‑143 (MQ=255)
              aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGtc                                                                                                                                    <  3:155867/149‑1 (MQ=255)
               gAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCC‑‑AGGCCGGGGtttttctt                                                                                                                                   <  8:34437/149‑5 (MQ=255)
                aaTGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGtcat                                                                                                                                  >  2:162090/1‑149 (MQ=255)
                aaTGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGtcat                                                                                                                                  >  6:209039/1‑149 (MQ=255)
                aaTGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCC‑‑AGGCCGGGGtttttctt                                                                                                                                   >  4:1218/1‑144 (MQ=255)
                    gCCGGCTTCAAGGACATGGCAGCCTAGCACTAAGGTGGTCTTAAGTCGAGTGTAGACGGGCATTGCCGAAGGTGTGCAACTGAAGGCGGGCTTTATGTCCCGGGCCACTAAAGAAAAACACCGGCCT‑‑GGGCCGGGGTTTTGTCATCAc                                                                                                                               >  4:188222/1‑148 (MQ=255)
                     ccGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACtc                                                                                                                             >  4:51114/1‑149 (MQ=255)
                           tCAAGGACATGGCAGCCTCGCAGTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                       >  7:62495/1‑149 (MQ=255)
                           tCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                       >  2:72717/1‑149 (MQ=255)
                            cAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGag                                                                                                                       <  7:360000/148‑1 (MQ=255)
                              aGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGaggaa                                                                                                                    <  5:209039/149‑1 (MQ=255)
                              aGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGagg                                                                                                                      >  1:152380/1‑147 (MQ=255)
                               ggACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACTGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGaggaag                                                                                                                   >  2:101790/1‑149 (MQ=255)
                                  cATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAg                                                                                                                >  5:326535/1‑149 (MQ=255)
                                   aTGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAgc                                                                                                               >  2:270527/1‑149 (MQ=255)
                                     ggCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAgcg                                                                                                              >  4:43693/1‑148 (MQ=255)
                                               cACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGc                                                                                                   >  2:31553/1‑149 (MQ=255)
                                                    aGGTGGTCTGAAGTCGAGTGTAGACTGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                                <  3:43693/147‑1 (MQ=255)
                                                    aGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATTTCACGGGCCACTAAAGAAAAACCCCGGCCTCCGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTg                                                                                                 <  1:101790/148‑1 (MQ=255)
                                                    aGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGACGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                                >  5:210233/1‑147 (MQ=255)
                                                       tggtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTc                                                                                           <  3:51114/149‑1 (MQ=255)
                                                         gtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                         >  1:52618/1‑149 (MQ=255)
                                                         gtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                         <  1:31553/149‑1 (MQ=255)
                                                            tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTc                                                                                                     >  3:303490/1‑134 (MQ=255)
                                                            tGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGCCGGGCTTTATGTCACGGGCCAATAAAGAAAAACCACGGCCA‑‑CGGCCGGGGTTTTCTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTccc                                                                                                   >  3:14807/1‑134 (MQ=255)
                                                                gTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGGCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGgcgc                                                                                  >  6:343925/1‑149 (MQ=255)
                                                                 tCGAGTGTAGACGGGCATTTCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCa                                                                                 <  1:102990/149‑1 (MQ=255)
                                                                  cGAGTGTAGAGGGGCATTGCCGCAGGTGTGCATGGGTAGGCGGGCTTTATGTCTCGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCGGGGGTTTTGTCATCAGTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCa                                                                                 <  2:355902/148‑1 (MQ=255)
                                                                   gAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGc                                                                               >  5:123452/1‑149 (MQ=255)
                                                                     gtgtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGa                                                                                                                        >  5:218821/1‑106 (MQ=255)
                                                                     gtgtAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGa                                                                                                                        <  6:218821/106‑1 (MQ=255)
                                                                            cGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAAc                                                                      <  4:334688/149‑1 (MQ=255)
                                                                            cGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAAc                                                                      <  1:162090/149‑1 (MQ=255)
                                                                                 aTTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCctt                                                                 <  1:72717/149‑1 (MQ=255)
                                                                                     ccGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGt                                                                                                                                          >  2:337972/1‑72 (MQ=255)
                                                                                     ccGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgct                                                             >  4:283183/1‑149 (MQ=255)
                                                                                      cGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgctt                                                            <  2:152380/149‑1 (MQ=255)
                                                                                       gCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgcttc                                                           >  4:201205/1‑149 (MQ=255)
                                                                                           gtgtGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATc                                                       >  1:130832/1‑149 (MQ=255)
                                                                                            tgtgCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCt                                                      <  6:123452/149‑1 (MQ=255)
                                                                                                 aaCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCt                                                      >  7:263614/1‑144 (MQ=255)
                                                                                                 aaCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCt                                                      <  8:263614/144‑1 (MQ=255)
                                                                                                 aaCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACg                                                  >  4:226639/1‑148 (MQ=255)
                                                                                                  aCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGa                                                                                                                        >  4:313091/1‑77 (MQ=255)
                                                                                                  aCGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGa                                                                                                                        <  3:313091/77‑1 (MQ=255)
                                                                                                   cGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGAt                                               >  4:147057/1‑149 (MQ=255)
                                                                                                    ggAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATc                                              >  6:5701/1‑149 (MQ=255)
                                                                                                         gCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCt                                          >  6:232724/1‑148 (MQ=255)
                                                                                                         gCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCt                                          >  5:374003/1‑148 (MQ=255)
                                                                                                            ggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGCCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGAAGGTGCTCGCTGCGAGTAGGCTGGCGCAGATTGCGCAACGCAATCGATTCGATCTGACGGATCCGCTCAcg                                      >  4:359766/1‑149 (MQ=255)
                                                                                                             gCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgc                                     <  3:201205/149‑1 (MQ=255)
                                                                                                               tttATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGc                                           <  6:50144/141‑1 (MQ=255)
                                                                                                               tttATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGc                                           >  5:50144/1‑141 (MQ=255)
                                                                                                               tttATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGt                                   >  6:283484/1‑149 (MQ=255)
                                                                                                                 tATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTa                                 >  2:260071/1‑149 (MQ=255)
                                                                                                                     tCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTc                             <  6:287708/149‑1 (MQ=255)
                                                                                                                      cACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGCGCTCGCCGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTc                             >  5:37628/1‑148 (MQ=255)
                                                                                                                         gggCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgc                                     >  4:405962/1‑137 (MQ=255)
                                                                                                                         gggCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAcgc                                     <  3:405962/137‑1 (MQ=255)
                                                                                                                              aCTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCtt                    <  3:226639/149‑1 (MQ=255)
                                                                                                                                tAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGAc                                                   >  2:146434/1‑116 (MQ=255)
                                                                                                                                tAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGAc                                                   <  1:146434/116‑1 (MQ=255)
                                                                                                                                tAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                  >  2:69039/1‑149 (MQ=255)
                                                                                                                                tAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                  >  8:249090/1‑149 (MQ=255)
                                                                                                                                  aaGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGccc                <  5:283484/149‑1 (MQ=255)
                                                                                                                                  aaGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc                 >  3:323730/1‑148 (MQ=255)
                                                                                                                                  aaGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                  <  2:130832/147‑1 (MQ=255)
                                                                                                                                  aaGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTCAAGGATCCGCTCACGCGTTAAGTCAAACTGCTTGcc                 >  7:366033/1‑148 (MQ=255)
                                                                                                                                     caaaaCCCCGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCGTCGATCTGACGGATCCGCTCACGCGGTACGTCAAACTGCTTGCCCAcc             >  3:124442/2‑149 (MQ=255)
                                                                                                                                     aaaaaCCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTGGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAcc             >  4:91031/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCtt                    >  7:5673/1‑141 (MQ=255)
                                                                                                                                      aaaaCCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCtt                    <  8:5673/141‑1 (MQ=255)
                                                                                                                                      aaaaCCCCGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCAcc             >  6:374895/1‑148 (MQ=255)
                                                                                                                                      aaaaCCCCGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACct            >  7:104244/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCCGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACct            >  7:111944/1‑149 (MQ=255)
                                                                                                                                      aaaaCCCCGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACct            >  3:349066/1‑149 (MQ=255)
                                                                                                                                         aCCCCGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCtt         >  6:352651/1‑149 (MQ=255)
                                                                                                                                          ccccGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGa                                                    >  6:134169/1‑105 (MQ=255)
                                                                                                                                          ccccGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGa                                                    <  5:134169/105‑1 (MQ=255)
                                                                                                                                          ccccGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc        >  4:187254/1‑149 (MQ=255)
                                                                                                                                          ccccGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc        <  5:352651/149‑1 (MQ=255)
                                                                                                                                          ccccGGCCC‑‑AGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc        >  4:239604/1‑149 (MQ=255)
                                                                                                                                             cGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAACCTGCTTGCCCACCTCTTCGAg     >  7:73951/1‑149 (MQ=255)
                                                                                                                                              ggacT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGTTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg    <  3:188222/146‑1 (MQ=255)
                                                                                                                                              ggCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg    <  6:326535/149‑1 (MQ=255)
                                                                                                                                              ggCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg    <  3:283183/149‑1 (MQ=255)
                                                                                                                                              ggCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg    <  5:5701/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCT‑‑GGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGT  >  NC_002947/470920‑471212

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