Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,715 T→G 6.6% I142L (ATC→CTC)  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,7150TG 6.6% 183.5 / 4.0 61I142L (ATC→CTC) PP_0162ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base G (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ATGGCCTGCCTGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCC  >  NC_002947/185569‑185857
                                                                                                                                                  |                                                                                                                                              
aTGGCCTGCCTGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATg                                                                                                                                               <  3:284049/148‑1 (MQ=255)
   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:332920/1‑149 (MQ=255)
                                                         gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                     <  4:60362/149‑1 (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                                                                       <  5:78802/149‑1 (MQ=255)
                                                                       gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                       <  1:219632/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                                                                                 >  8:176606/1‑108 (MQ=255)
                                                                                                      gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                 <  7:176606/108‑1 (MQ=255)
                                                                                                         gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                           <  8:23314/143‑1 (MQ=255)
                                                                                                         gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                           >  7:23314/1‑143 (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                 <  1:219328/142‑1 (MQ=255)
                                                                                                                                    cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                 >  2:219328/1‑142 (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  <  5:249693/147‑1 (MQ=255)
                                                                                                                                            atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc  >  6:249693/3‑149 (MQ=255)
                                                                                                                                                gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                             <  6:209370/118‑1 (MQ=255)
                                                                                                                                                gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                             >  5:209370/1‑118 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ATGGCCTGCCTGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCC  >  NC_002947/185569‑185857

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