Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,731 C→G 6.7% G136G (GGG→GGC PP_0162 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947185,7310CG 6.7% 242.5 / 6.2 74G136G (GGG→GGCPP_0162ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base G (3/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTC  >  NC_002947/185583‑185877
                                                                                                                                                     |                                                                                                                                                  
cGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGc                                                                                                                                                    <  4:160496/149‑1 (MQ=255)
cGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGc                                                                                                                                                    <  4:228370/149‑1 (MQ=255)
 gTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGcc                                                                                                                                                   >  7:338512/1‑149 (MQ=255)
    tttGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGcccc                                                                                                                                                 <  4:322063/148‑1 (MQ=255)
    tttGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCt                                                                                                                                                >  1:84861/1‑149 (MQ=255)
      tGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCAATGGCCTGATGCTTCAGCCCCTc                                                                                                                                               <  7:80166/148‑1 (MQ=255)
         aTCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCTGCCCCTCGCAt                                                                                                                                           >  5:245221/1‑149 (MQ=255)
          tCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCCGGCGATGGCCTGATGCTTCAGCCCCTCGCa                                                                                                                                            <  1:33355/148‑1 (MQ=255)
            atggCGTGCGACGGTCATTGATTCAGCAGTCGTCTTCATAGCGCAGCAAGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGa                                                                                                                                        <  7:105157/147‑1 (MQ=255)
             gggCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGc                                                                                                                                             <  1:326805/143‑1 (MQ=255)
             gggCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGc                                                                                                                                             >  2:326805/1‑143 (MQ=255)
                 gTGCGACGTTCTTGGATTCAGCAGTCGTATTCATAGCGCTGCAGGTAGCAGTGATTGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCGATCTGTTCGGCGATGGCCTGATGCTTCTGCCCCTCGCATTGAGCcag                                                                                                                                   <  4:294979/149‑1 (MQ=255)
                    cGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCcagcag                                                                                                                                <  6:175912/149‑1 (MQ=255)
                     gACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCcagcag                                                                                                                                >  2:7407/1‑148 (MQ=255)
                            aTGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATgcgc                                                                                                                        <  4:289433/149‑1 (MQ=255)
                               gATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGAc                                                                                                                     <  8:338512/149‑1 (MQ=255)
                                    agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCACCGGTGCGCACGGAGATCCCCATCTGTTCGGCCATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                                >  6:127063/1‑149 (MQ=255)
                                    agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                                >  8:94961/1‑149 (MQ=255)
                                    agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                                >  5:332770/1‑149 (MQ=255)
                                       agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                             <  6:245221/149‑1 (MQ=255)
                                         tcgtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTggg                                                                                                           >  3:362509/1‑149 (MQ=255)
                                              ttCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGttt                                                                                                        >  4:312251/1‑147 (MQ=255)
                                                cATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGa                                                                                                      <  4:258277/147‑1 (MQ=255)
                                                cATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGa                                                                                                      <  3:312251/147‑1 (MQ=255)
                                                 aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGc                                                                                                   >  1:159141/1‑149 (MQ=255)
                                                   aGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCg                                                                                                 >  1:147066/1‑149 (MQ=255)
                                                    gcgcAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGt                                                                                                <  6:332770/149‑1 (MQ=255)
                                                      gcagcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTc                                                                                               >  6:256356/1‑148 (MQ=255)
                                                        agcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAg                                                                                            >  6:23189/1‑149 (MQ=255)
                                                             gTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGc                                                                                       >  6:160878/1‑149 (MQ=255)
                                                                 cAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGc                                                                                   <  1:7407/149‑1 (MQ=255)
                                                                   gTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGa                                                                                                      <  3:256739/128‑1 (MQ=255)
                                                                   gTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGa                                                                                                      >  4:256739/1‑128 (MQ=255)
                                                                      aTAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATAc                                                                              >  4:146373/1‑149 (MQ=255)
                                                                       tAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGc                                                                                                   >  2:280154/1‑127 (MQ=255)
                                                                       tAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGc                                                                                                   <  1:280154/127‑1 (MQ=255)
                                                                       tAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACg                                                                             <  5:148481/149‑1 (MQ=255)
                                                                         gagCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCt                                                                           <  3:146373/149‑1 (MQ=255)
                                                                            cgcATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCaa                                                                        <  5:256356/149‑1 (MQ=255)
                                                                            cgcATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCaa                                                                        <  8:121223/149‑1 (MQ=255)
                                                                              cATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGc                                                                      <  4:74405/149‑1 (MQ=255)
                                                                                  ggCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAAt                                                                   <  4:219153/148‑1 (MQ=255)
                                                                                     gATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                               <  7:94961/149‑1 (MQ=255)
                                                                                      aTGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAAt                                                                   >  7:119594/1‑144 (MQ=255)
                                                                                      aTGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAAt                                                                   <  8:119594/144‑1 (MQ=255)
                                                                                         tAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCcagcag                                                           <  2:84861/149‑1 (MQ=255)
                                                                                             gcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCcagca                                                                                                                                 <  6:27473/75‑1 (MQ=255)
                                                                                             gcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCagcagcagc                                                       >  5:362497/1‑149 (MQ=255)
                                                                                              cTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCATGCAATTCCAGCAGCAGCt                                                      >  4:177504/1‑149 (MQ=255)
                                                                                                   cGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                                 <  5:75185/149‑1 (MQ=255)
                                                                                                          cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                                          <  1:30944/149‑1 (MQ=255)
                                                                                                           aCGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGAc                                         >  2:115128/1‑149 (MQ=255)
                                                                                                               agaTCCCCATCTGTTCGGCGATGGCCTTATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                     <  5:370785/149‑1 (MQ=255)
                                                                                                                  tCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAggt                                   <  5:127063/148‑2 (MQ=255)
                                                                                                                   ccccATATGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                                 >  3:199032/1‑149 (MQ=255)
                                                                                                                    cccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                                 <  4:362509/148‑1 (MQ=255)
                                                                                                                    cccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                                >  3:166953/1‑149 (MQ=255)
                                                                                                                    cccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                                <  5:112970/149‑1 (MQ=255)
                                                                                                                     ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCcag                                                              >  2:55471/1‑118 (MQ=255)
                                                                                                                     ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                                >  2:100729/1‑148 (MQ=255)
                                                                                                                      cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCa                              <  5:160878/149‑1 (MQ=255)
                                                                                                                       aTCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTcgcg                                                  >  7:126481/1‑128 (MQ=255)
                                                                                                                       aTCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTcgcg                                                  <  8:126481/128‑1 (MQ=255)
                                                                                                                           gTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAgaggggacgcg                           >  1:312283/1‑140 (MQ=255)
                                                                                                                           gTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                     <  2:312283/137‑1 (MQ=255)
                                                                                                                              aatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTc                       >  5:125738/4‑149 (MQ=255)
                                                                                                                                ggCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCg                     <  4:199032/149‑1 (MQ=255)
                                                                                                                                  cGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCgg                    >  5:267059/1‑148 (MQ=255)
                                                                                                                                  cGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGc                   <  6:125738/149‑1 (MQ=255)
                                                                                                                                   gaggggcTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGca                  >  3:311674/7‑149 (MQ=255)
                                                                                                                                      gggcTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACAt               <  6:267059/146‑1 (MQ=255)
                                                                                                                                      ggCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACAt               >  6:105940/1‑149 (MQ=255)
                                                                                                                                        ccTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATgg             <  1:100729/149‑1 (MQ=255)
                                                                                                                                         cTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATggg            >  7:31036/1‑149 (MQ=255)
                                                                                                                                         cTGATGCTTCAGCCCCTCGCATGGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATggg            >  7:237189/1‑149 (MQ=255)
                                                                                                                                            aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGcc          <  5:105940/148‑1 (MQ=255)
                                                                                                                                                ttCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGc     >  3:365521/1‑149 (MQ=255)
                                                                                                                                                   aGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGGTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTc  <  2:159141/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTC  >  NC_002947/185583‑185877

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