Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,897,179 T→C 6.5% intergenic (‑32/+112) PP_mr24 ← / ← rdgC mini‑ykkC regulatory region/recombination‑associated protein RdgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,897,1790TC 6.5% 177.7 / 2.5 61intergenic (‑32/+112)PP_mr24/rdgCmini‑ykkC regulatory region/recombination‑associated protein RdgC
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/2);  total (30/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.21e-01

CCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACC  >  NC_002947/1897034‑1897320
                                                                                                                                                 |                                                                                                                                             
gcTCGATCAGGGCGGCTAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCAGGGTGCCTGTCTGCGGCCTTGCg                                                                                                                                            <  2:380378/148‑1 (MQ=255)
 cTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCgg                                                                                                                                           <  7:166078/149‑1 (MQ=255)
 cTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTCGCgg                                                                                                                                           >  4:288341/1‑149 (MQ=255)
   cGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGGGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGc                                                                                                                                         <  7:86810/149‑1 (MQ=255)
    gAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGcccgc                                                                                                                                         <  2:374036/148‑5 (MQ=255)
    gAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGc                                                                                                                                         <  7:38148/148‑1 (MQ=255)
    gAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGc                                                                                                                                         >  7:307704/1‑148 (MQ=255)
     aaCAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCTCCATGGTGCCTGTCTGCGGCCTTGCGGGCaa                                                                                                                                       <  3:300019/149‑1 (MQ=255)
         ggCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGAca                                                                                                                                   <  2:31199/149‑2 (MQ=255)
              gCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACcg                                                                                                                                  <  8:238425/145‑1 (MQ=255)
              gCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACcg                                                                                                                                  >  7:238425/1‑145 (MQ=255)
                gAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTCGCGGGTAAACCCGCGCCTg                                                                                                                            >  4:190846/1‑149 (MQ=255)
                 aaaaaCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGATCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGc                                                                                                                           <  1:232140/149‑1 (MQ=255)
                         gatgatCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACAATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGt                                                                                                                    >  1:541/1‑148 (MQ=255)
                              tCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCAc                                                                                                              >  5:180874/1‑149 (MQ=255)
                              tCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGACCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCAc                                                                                                              >  5:137475/1‑149 (MQ=255)
                                  ggACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGc                                                                                                          <  7:308819/149‑1 (MQ=255)
                                  ggACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGc                                                                                                          <  3:402782/149‑1 (MQ=255)
                                   gACATGAGGGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCa                                                                                                         <  1:263282/149‑1 (MQ=255)
                                    aCATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTCGCGGGTAAACCCGCGCCTGCAAGTGGTAACCCCTTGCAg                                                                                                        <  3:190846/149‑1 (MQ=38)
                                        gAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGgcgc                                                                                                    >  7:88083/1‑149 (MQ=255)
                                         aGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGgcgc                                                                                                    <  5:406802/148‑1 (MQ=255)
                                         aGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGcgcg                                                                                                   <  5:169839/149‑1 (MQ=255)
                                               tctcCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                      >  5:4409/1‑124 (MQ=255)
                                               tctcCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                      <  6:4409/124‑1 (MQ=255)
                                                  ccATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTAccc                                                                                          >  1:31976/1‑149 (MQ=255)
                                                        tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACcgcg                                                                                                                                <  6:196107/105‑1 (MQ=255)
                                                        tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACcgcg                                                                                                                                >  5:196107/1‑105 (MQ=255)
                                                        tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGcgcg                                                                                                   <  3:283810/134‑1 (MQ=255)
                                                        tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGcgcg                                                                                                   >  4:283810/1‑134 (MQ=255)
                                                        tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTCGCGGGTAAACCCGCGCCTGCAAGTGGTAACCCCTTGCAGGCGCGGTCTTGCCcgcaagg                                                                                    <  3:288341/149‑5 (MQ=21)
                                                              tgatgGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCg                                                                              >  2:247914/1‑149 (MQ=255)
                                                                 tgGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTGTCCCGCGAAGAGGCCGGTg                                                                           >  4:268488/1‑149 (MQ=255)
                                                                 tgGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCAACATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTg                                                                           >  8:385664/1‑149 (MQ=255)
                                                                       cGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCa                                                                     <  8:307704/149‑1 (MQ=255)
                                                                         tCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATCGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGAAGGCGCGGGTTTACCAGCGAAGAGGCCGGTGCTCGCAg                                                                    >  3:218942/1‑148 (MQ=255)
                                                                                    gTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAg                                                        >  5:331317/1‑149 (MQ=255)
                                                                                           ggTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCa                                                 >  7:175542/1‑149 (MQ=255)
                                                                                                 cggtGAGGACGACCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTgg                                                                       <  4:349862/119‑1 (MQ=255)
                                                                                                      aGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCa                                      <  6:331317/149‑1 (MQ=255)
                                                                                                       ggTCGTCCTCACCGTGGGCACCATGGTGCCTATCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCa                                      >  6:102553/1‑148 (MQ=255)
                                                                                                          cgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCa                                                         >  6:314870/1‑126 (MQ=255)
                                                                                                          cgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCa                                                         <  5:314870/126‑1 (MQ=255)
                                                                                                              cTCACCGTGGGCACCATGGTGCCTATCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAAt                               <  5:102553/148‑1 (MQ=255)
                                                                                                               tCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTAc                                                                                            <  1:46285/86‑1 (MQ=255)
                                                                                                               tCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTAc                                                                                            >  2:46285/1‑86 (MQ=255)
                                                                                                               tCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAAt                               >  2:177734/1‑147 (MQ=255)
                                                                                                                cACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAAt                               >  5:369585/1‑146 (MQ=255)
                                                                                                                  ccGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGAc                          >  3:92997/1‑149 (MQ=255)
                                                                                                                     tGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGcc                                                     >  4:290044/1‑119 (MQ=255)
                                                                                                                     tGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGcc                                                     <  3:290044/119‑1 (MQ=255)
                                                                                                                        gCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTg                    >  4:16680/1‑149 (MQ=255)
                                                                                                                         cACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCt                                                                         >  7:106163/1‑95 (MQ=255)
                                                                                                                         cACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCt                                                                         <  8:106163/95‑1 (MQ=255)
                                                                                                                            cATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCa                                                    >  2:349608/1‑113 (MQ=255)
                                                                                                                            cATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCa                                                    <  1:349608/113‑1 (MQ=255)
                                                                                                                                gTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAAtc             <  2:31976/148‑1 (MQ=255)
                                                                                                                                    cTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTc        <  1:177734/149‑1 (MQ=255)
                                                                                                                                       tCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTca       <  3:16680/147‑1 (MQ=255)
                                                                                                                                          gCGGTCTTGCCCGCAAGGCCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCaccacc  >  5:232776/1‑149 (MQ=255)
                                                                                                                                           cGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCaccacc  >  2:16524/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACC  >  NC_002947/1897034‑1897320

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: