Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,718 C→G 6.8% A141P (GCC→CCC)  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,7180CG 6.8% 187.0 / 3.1 59A141P (GCC→CCC) PP_0162ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base G (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCCTGCCTGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGG  >  NC_002947/185572‑185866
                                                                                                                                                  |                                                                                                                                                    
gcctgcctGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCt                                                                                                                                                    >  6:269803/1‑149 (MQ=255)
  ctgcctgACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGa                                                                                                                                                  <  7:195891/149‑1 (MQ=255)
       tgACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTAAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGAACGGAGATCCCCATCTCTTCGGCGATGGCCTGATGCt                                                                                                                                              >  6:48093/1‑148 (MQ=255)
        gACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCt                                                                                                                                                    <  4:257249/141‑1 (MQ=255)
        gACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCt                                                                                                                                                    >  3:257249/1‑141 (MQ=255)
         aCCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCa                                                                                                                                           >  8:78728/1‑149 (MQ=255)
          ccGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCATCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTAAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCa                                                                                                                                           <  5:48093/148‑1 (MQ=255)
          ccGTGTTTGTATCAGGGCGGGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAg                                                                                                                                          <  4:296186/149‑1 (MQ=255)
          acttttttGTATCAGGGCGTGCGACGGTAATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCa                                                                                                                                           <  8:91005/145‑1 (MQ=255)
                ttGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCt                                                                                                                                     <  7:274354/148‑1 (MQ=255)
                       aGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGAGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGCTCCCCATCTGTTCGGCTATGGCCTGATGCTTCAGCCCCTCGCATTGa                                                                                                                             <  5:78449/149‑1 (MQ=255)
                          gCGTGCGACGGTCATGGATTCAGCAGTGTTCTTCATAGCGCAGCACGTAGCAGTTATAGAGCGCATTGGCGATGTATCGCTCAAGGGTGCGCAGGGATATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGc                                                                                                                           <  7:171213/148‑1 (MQ=255)
                           cGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCca                                                                                                                         <  2:183490/149‑1 (MQ=255)
                                        tGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGcgcg                                                                                                            <  5:269803/149‑1 (MQ=255)
                                          gATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGAc                                                                                                          >  6:78802/1‑149 (MQ=255)
                                               agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                     >  6:247877/1‑149 (MQ=255)
                                                 cagTCGTCTTCATAGCGCAGCACCTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCtt                                                                                                   >  2:63344/1‑149 (MQ=255)
                                                  agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                  >  6:46455/1‑149 (MQ=255)
                                                  agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                  <  1:133492/149‑1 (MQ=255)
                                                    tcgtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGcc                                                                                                                                                     >  7:268759/1‑96 (MQ=255)
                                                      gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                              >  1:322943/1‑149 (MQ=255)
                                                      gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                              <  4:60362/149‑1 (MQ=255)
                                                      gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                              >  4:332920/1‑149 (MQ=255)
                                                       tcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGttt                                                                                             <  2:322943/149‑1 (MQ=255)
                                                         ttCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGAACCTTGGGttt                                                                                             >  7:54499/1‑147 (MQ=255)
                                                            aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAgg                                                                                         <  3:332920/148‑1 (MQ=255)
                                                            aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAgg                                                                                         <  5:46455/148‑1 (MQ=255)
                                                                 gcagcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCACGAATGCGCGACGTACCTAGGGATTGAGGCCGTc                                                                                    >  3:297832/1‑148 (MQ=255)
                                                                    gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                <  1:219632/149‑1 (MQ=255)
                                                                    gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                <  5:78802/149‑1 (MQ=255)
                                                                           gCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCGCATCTGTTCGGCGCTGGCCTTATGCTTCAGCCCCTCGCATTGAGCCAGCATGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGt                                                                          <  1:158109/148‑1 (MQ=255)
                                                                                 aTAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATAc                                                                   >  7:134791/1‑149 (MQ=255)
                                                                                  tAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACg                                                                  >  2:257864/1‑149 (MQ=255)
                                                                                           ttGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATCCTTCAGCCCCACGCATTTAGCCAGCAGGAATCCCCGACGTACCTTGGGTTTGAGCCCGACGAGCATGCCTCCAATACGCTCAAGCaa                                                         >  3:277777/1‑149 (MQ=255)
                                                                                              gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTgg                                                                                                 <  2:257045/106‑1 (MQ=255)
                                                                                              gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTc                                                      >  3:178895/1‑149 (MQ=255)
                                                                                              gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGCCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTgg                                                                                                 >  1:257045/1‑106 (MQ=255)
                                                                                                gATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                    <  4:178895/149‑1 (MQ=255)
                                                                                                   gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                          <  7:176606/108‑1 (MQ=255)
                                                                                                   gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                          >  8:176606/1‑108 (MQ=255)
                                                                                                      gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                    >  7:23314/1‑143 (MQ=255)
                                                                                                      gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                    <  8:23314/143‑1 (MQ=255)
                                                                                                                    gcACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCTCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGtctc                                 >  8:139020/1‑148 (MQ=255)
                                                                                                                    gcACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGtctc                                 >  7:209207/1‑148 (MQ=255)
                                                                                                                     cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                               >  3:152935/1‑149 (MQ=255)
                                                                                                                     cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                               <  4:297832/149‑1 (MQ=255)
                                                                                                                              ccccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCCCGACGAACCAGGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                      >  8:184928/1‑149 (MQ=255)
                                                                                                                               cccATCTGTTCGGCGATGGCGTGATGCTTCTGCCCCTCGCATTGAGCCAGCAGGAATTCGCGACGTACCTTGTGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTGCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                     <  4:277777/149‑1 (MQ=255)
                                                                                                                                ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAcc                    <  1:63344/149‑1 (MQ=255)
                                                                                                                                ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                     >  1:284068/1‑148 (MQ=255)
                                                                                                                                 cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                          >  2:219328/1‑142 (MQ=255)
                                                                                                                                 cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                          <  1:219328/142‑1 (MQ=255)
                                                                                                                                  aTCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAg                  <  8:209207/149‑1 (MQ=255)
                                                                                                                                     ctcaatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTg                <  7:178506/142‑1 (MQ=255)
                                                                                                                                     ctcaatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTg                >  8:178506/7‑148 (MQ=255)
                                                                                                                                       ttCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCATGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGt                 <  1:257864/145‑1 (MQ=255)
                                                                                                                                         atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc           >  6:249693/3‑149 (MQ=255)
                                                                                                                                         atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc           <  5:249693/147‑1 (MQ=255)
                                                                                                                                             gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                      <  6:209370/118‑1 (MQ=255)
                                                                                                                                             gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                      >  5:209370/1‑118 (MQ=255)
                                                                                                                                                  ccTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCTAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATgg  <  7:78728/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCTGCCTGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGG  >  NC_002947/185572‑185866

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: ‑