Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,621 T→C 6.8% intergenic (+44/+47) ssb → / ← PP_0486 single‑stranded DNA‑binding protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,6210TC 6.8% 552.9 / 20.1 177intergenic (+44/+47)ssb/PP_0486single‑stranded DNA‑binding protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (78/86);  new base C (9/3);  total (88/89)
Fisher's exact test for biased strand distribution p-value = 7.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 179 total aligned reads displayed.

GCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGCACA  >  NC_002947/571476‑571770
                                                                                                                                                 |                                                                                                                                                     
gCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAAGAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCTGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGcc                                                                                                                                                   <  2:2061421/150‑1 (MQ=255)
  ccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCt                                                                                                                                                 >  1:4387093/1‑150 (MQ=255)
  ccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCt                                                                                                                                                 >  1:1399679/1‑150 (MQ=255)
   cAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGATGCCCTg                                                                                                                                                >  2:3564397/1‑150 (MQ=255)
    aGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTgg                                                                                                                                               <  1:2227426/150‑1 (MQ=255)
    aGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTgg                                                                                                                                               <  2:421890/150‑1 (MQ=255)
     gcgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTggg                                                                                                                                              <  1:1065458/150‑1 (MQ=255)
     gcgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTggg                                                                                                                                              <  1:3447366/150‑1 (MQ=255)
     gcgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTggg                                                                                                                                              <  1:3447357/150‑1 (MQ=255)
      cgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCCGGGCAGTTGCCCTGGGt                                                                                                                                             <  2:2814338/150‑1 (MQ=255)
           gCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGtttttt                                                                                                                                        >  1:4418719/1‑150 (MQ=255)
             aaCAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAAAAACCCAGGGCAACTGCCCTGGGtttctggt                                                                                                                                      <  2:3457901/150‑6 (MQ=255)
               cAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGGAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATg                                                                                                                                    >  2:1616634/1‑150 (MQ=255)
                 ggCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATggg                                                                                                                                  >  1:2284701/1‑150 (MQ=255)
                  gCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgc                                                                                                                                 <  1:3488168/150‑1 (MQ=255)
                  gCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgc                                                                                                                                 >  1:216206/1‑150 (MQ=255)
                   ccccTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgca                                                                                                                                >  2:2733453/1‑150 (MQ=255)
                    cccTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgcag                                                                                                                               <  1:4243546/150‑1 (MQ=255)
                       tCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCAc                                                                                                                            >  2:1601131/1‑150 (MQ=255)
                          gcgcccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACAAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACcgc                                                                                                                         >  2:413112/1‑150 (MQ=255)
                              ccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTAt                                                                                                                     <  1:3588552/150‑1 (MQ=255)
                                agcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGTCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCa                                                                                                                   <  2:1512453/150‑1 (MQ=255)
                                agcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCa                                                                                                                   <  1:4138620/150‑1 (MQ=255)
                                 gcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAg                                                                                                                  <  1:1502222/150‑1 (MQ=255)
                                 gcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAg                                                                                                                  <  2:1749053/150‑1 (MQ=255)
                                  cgccgcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAgg                                                                                                                 <  2:3135117/150‑1 (MQ=255)
                                  cgccgcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAgg                                                                                                                 <  2:553917/150‑1 (MQ=255)
                                   gccgcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGc                                                                                                                >  1:1737518/1‑150 (MQ=255)
                                      gcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATc                                                                                                             <  1:2047072/150‑1 (MQ=255)
                                        aaCAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGc                                                                                                           >  2:774278/1‑150 (MQ=255)
                                            gCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTc                                                                                                       <  2:1378881/150‑1 (MQ=255)
                                             ccTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCa                                                                                                      <  2:3734573/150‑1 (MQ=255)
                                             ccTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCa                                                                                                      >  2:2482809/1‑150 (MQ=255)
                                                gcgcCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTa                                                                                                   <  2:1702599/150‑1 (MQ=255)
                                                 cgccgcAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAg                                                                                                  <  1:78701/150‑1 (MQ=255)
                                                       aGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGaacaac                                                                                            <  1:3381547/150‑1 (MQ=255)
                                                        gCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAAcc                                                                                           >  2:4367561/1‑150 (MQ=255)
                                                          cggcggCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCGGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCgg                                                                                         <  1:2649602/150‑1 (MQ=255)
                                                          cggcggCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCgg                                                                                         >  1:1551433/1‑150 (MQ=255)
                                                             cggcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGTCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACACCCGGGGc                                                                                      <  2:802392/150‑1 (MQ=255)
                                                               gcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTCGAACAACCGGGGCTc                                                                                    >  1:2610428/1‑150 (MQ=255)
                                                                  gACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCgg                                                                                 <  1:3069341/150‑1 (MQ=255)
                                                                     ttCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGa                                                                              >  2:969425/1‑150 (MQ=255)
                                                                     ttCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGa                                                                              <  2:298885/150‑1 (MQ=255)
                                                                     ttCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGa                                                                              >  2:1441635/1‑150 (MQ=255)
                                                                      tCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGaa                                                                             <  2:1839791/150‑1 (MQ=255)
                                                                        gACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATc                                                                           <  1:2290677/150‑1 (MQ=255)
                                                                         aCAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATcc                                                                          >  1:1554205/1‑150 (MQ=255)
                                                                          cAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCca                                                                         <  1:1908187/150‑1 (MQ=255)
                                                                              ttCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcagcag                                                                     >  1:3991044/1‑150 (MQ=255)
                                                                               tCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAgg                                                                    <  1:4029891/150‑1 (MQ=255)
                                                                               tCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAgg                                                                    <  2:4308806/150‑1 (MQ=255)
                                                                                  aTGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCtt                                                                 <  1:3265089/150‑1 (MQ=255)
                                                                                     aCGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGGTCCGGCGAATCCAGCAGGCTTGcc                                                              <  1:2890799/150‑1 (MQ=255)
                                                                                       gATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAg                                                            >  1:4050969/1‑150 (MQ=255)
                                                                                        atatTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGc                                                           <  1:4296760/150‑1 (MQ=255)
                                                                                          atTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCtt                                                         >  1:2423660/1‑150 (MQ=255)
                                                                                            tCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTgg                                                       >  1:2582451/1‑150 (MQ=255)
                                                                                              cGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTc                                                     >  1:1797015/1‑150 (MQ=255)
                                                                                               gTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCa                                                    >  2:571295/1‑150 (MQ=255)
                                                                                                ttCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAg                                                   >  1:1295964/1‑150 (MQ=255)
                                                                                                 tCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAgg                                                  <  2:3892173/150‑1 (MQ=255)
                                                                                                      tCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCaaa                                             >  1:736024/1‑150 (MQ=255)
                                                                                                        agagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCACAGc                                           >  1:3861593/1‑150 (MQ=255)
                                                                                                        agagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGc                                           >  1:1591543/1‑150 (MQ=255)
                                                                                                         gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                          >  2:4338303/1‑150 (MQ=255)
                                                                                                         gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                          <  1:1094196/150‑1 (MQ=255)
                                                                                                           gCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCa                                        <  1:4350922/150‑1 (MQ=255)
                                                                                                            cGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAg                                       <  1:1967161/150‑1 (MQ=255)
                                                                                                                 gggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                                  <  2:22335/150‑1 (MQ=255)
                                                                                                                 gggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                                  <  2:3869960/150‑1 (MQ=255)
                                                                                                                     gTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTggcg                              >  1:2680681/1‑150 (MQ=255)
                                                                                                                     gTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTggcg                              >  1:2792393/1‑150 (MQ=255)
                                                                                                                       tCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGAACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTggcggc                            >  2:3360416/1‑150 (MQ=255)
                                                                                                                         aGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgca                                                                                                                                >  2:1116485/1‑48 (MQ=255)
                                                                                                                           cAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTggg                        <  2:3828219/150‑1 (MQ=255)
                                                                                                                            aaGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGt                       <  1:3686295/150‑1 (MQ=255)
                                                                                                                            aaGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGt                       <  2:3375919/150‑1 (MQ=255)
                                                                                                                             aGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTg                      >  1:798426/1‑150 (MQ=255)
                                                                                                                              ataaaaaaaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGc                     >  1:779870/7‑150 (MQ=255)
                                                                                                                              ataaaaaaaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCACTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGc                     >  1:1670508/7‑150 (MQ=255)
                                                                                                                               taaaaaaaCCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGcc                    >  2:1228627/6‑150 (MQ=255)
                                                                                                                                ccAGAAACCCAGGGCAGTTGCCCTGGGATTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCa                   <  2:2559892/150‑1 (MQ=255)
                                                                                                                                 cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCAt                  >  1:2362945/1‑150 (MQ=255)
                                                                                                                                 cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCAt                  >  2:2372566/1‑150 (MQ=255)
                                                                                                                                   gacaCCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCa                <  2:1976285/147‑1 (MQ=255)
                                                                                                                                     aaCCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAgg              <  1:143898/150‑1 (MQ=255)
                                                                                                                                      aCCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg             <  1:1832908/150‑1 (MQ=255)
                                                                                                                                      aCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg             >  2:4280357/1‑150 (MQ=255)
                                                                                                                                       cccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg           >  1:1252277/1‑150 (MQ=255)
                                                                                                                                       cccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc            >  2:3413113/1‑150 (MQ=255)
                                                                                                                                        ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg           >  1:3532469/1‑150 (MQ=255)
                                                                                                                                        ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGt                                                                                                    <  1:2820163/61‑1 (MQ=255)
                                                                                                                                            ggCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATgc       >  1:3594038/1‑150 (MQ=255)
                                                                                                                                            ggCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATgc       <  1:3564397/150‑1 (MQ=255)
                                                                                                                                            ggCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATgc       <  2:2493078/150‑1 (MQ=255)
                                                                                                                                               aGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGca    <  2:4129027/150‑1 (MQ=255)
                                                                                                                                               aGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGca    >  2:288168/1‑150 (MQ=255)
                                                                                                                                                gTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGcac   >  1:885945/1‑150 (MQ=255)
                                                                                                                                                 ttGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGcaca  <  2:4418719/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGCACA  >  NC_002947/571476‑571770

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