Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,804,832 C→G 5.1% A201G (GCC→GGC)  rbsR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,804,8320CG 5.1% 282.1 / 3.2 90A201G (GCC→GGC) rbsRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base G (2/3);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCTTTG  >  NC_002947/2804684‑2804978
                                                                                                                                                    |                                                                                                                                                  
gATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGcc                                                                                                                                                   >  1:743225/1‑150 (MQ=255)
   cccGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGCCGTCGCCTATATCGGCGGCCCCCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATgg                                                                                                                                               <  1:1345614/149‑1 (MQ=255)
      gACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGcc                                                                                                                                             >  1:621391/1‑150 (MQ=255)
         cTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGaa                                                                                                                                          <  1:590647/150‑1 (MQ=255)
         cTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGaa                                                                                                                                          <  2:776446/150‑1 (MQ=255)
            gTGCGCATTGACCACGAGCAGGGGGCGTATCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCGGCCATGGCACGGCGGAAGcc                                                                                                                                       <  2:1821860/150‑1 (MQ=255)
             tGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCg                                                                                                                                      <  1:396505/150‑1 (MQ=255)
                gcATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTg                                                                                                                                   >  1:1610246/1‑150 (MQ=255)
                 cATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGCACCCGGGGTTACCCCGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGtcg                                                                                                                                  >  2:898006/1‑148 (MQ=255)
                   ttGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCCGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTgcgc                                                                                                                                >  1:409365/1‑150 (MQ=255)
                    tGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGcc                                                                                                                               <  2:686617/150‑1 (MQ=255)
                       ccACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCGGCCATGGCACGGCGGAAGCCgctca                                                                                                                                  <  2:1855860/144‑5 (MQ=255)
                       ccACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCGGCCATGGCACGGCGGAAGCCgctca                                                                                                                                  >  1:1855860/1‑140 (MQ=255)
                       ccACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGt                                                                                                                            <  2:291806/150‑1 (MQ=255)
                       ccACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGt                                                                                                                            <  2:138234/150‑1 (MQ=255)
                        cACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTACAGGCATCCGTCTGTTTCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTg                                                                                                                           >  1:331073/1‑150 (MQ=255)
                         aCGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGc                                                                                                                          >  1:878311/1‑150 (MQ=255)
                            aGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCgg                                                                                                                       <  1:1403581/150‑1 (MQ=255)
                            aGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCgg                                                                                                                       <  1:510058/150‑1 (MQ=255)
                               aGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCa                                                                                                                    >  2:559129/1‑150 (MQ=255)
                               aGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCa                                                                                                                    >  2:558833/1‑150 (MQ=255)
                                   ggCGTACCTGGCCACGCGGCATCTGCAGGAGCTTGGCCATCGCGACGTCGCGTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTCCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCg                                                                                                                >  2:500556/1‑150 (MQ=255)
                                    gCGTACCTGGCCATCCGGCATCTCCTTGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCgt                                                                                                               <  1:392143/150‑1 (MQ=255)
                                    gCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCgt                                                                                                               <  1:184314/150‑1 (MQ=255)
                                    gCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCgt                                                                                                               <  2:621391/150‑1 (MQ=255)
                                       tACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTa                                                                                                            <  1:158936/150‑1 (MQ=255)
                                          cTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTactg                                                                                                         >  1:16530/1‑150 (MQ=255)
                                          cTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTactg                                                                                                         >  1:1632833/1‑150 (MQ=255)
                                          cTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACATAACTGCCCCTGAGCGGCTTCCGCCGTGCCATGGCAGCAGCTGGTGCGCCAGTGCCTTGCAGCCGTGGacgg                                                                                                         >  1:1860795/1‑148 (MQ=255)
                                           tGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTactgc                                                                                                        >  1:343031/1‑150 (MQ=255)
                                            ggCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTactgca                                                                                                       <  1:1742965/150‑1 (MQ=255)
                                                  ccGGCATCTGCTGGAGCTTGGCCATGGCGCCGGCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGACTGAACGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGa                                                                                                 >  1:527691/1‑150 (MQ=255)
                                                  ccGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGa                                                                                                 >  2:1145548/1‑150 (MQ=255)
                                                  ccGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGa                                                                                                 >  1:1251208/1‑150 (MQ=255)
                                                  ccGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGCCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGa                                                                                                 <  1:1564865/150‑1 (MQ=255)
                                                            ctgGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGCGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCa                                                                                         <  1:1592611/148‑1 (MQ=255)
                                                            ctgGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGCGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCa                                                                                         >  2:1592611/1‑148 (MQ=255)
                                                              ggAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCCCTCCATTCACCACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGccc                                                                                     <  2:527691/150‑1 (MQ=255)
                                                               gAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCg                                                                                    >  1:969446/1‑150 (MQ=255)
                                                                ggCTTGGCCATCGCGATGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCgg                                                                                   <  2:272075/149‑1 (MQ=255)
                                                                   ttGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGACGCAACCGATGTTACCTCGCTTCGCCGTTCGGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTg                                                                                >  2:65565/1‑150 (MQ=255)
                                                                            gcgACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCa                                                                                         >  2:658102/1‑132 (MQ=255)
                                                                            gcgACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCa                                                                                         <  1:658102/132‑1 (MQ=255)
                                                                             cgACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgc                                                                      >  2:1244184/1‑150 (MQ=255)
                                                                                  tCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcgg                                                                 >  2:795620/1‑150 (MQ=255)
                                                                                  tCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcgg                                                                 >  1:1023895/1‑150 (MQ=255)
                                                                                    gCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcggcg                                                               >  2:541767/1‑150 (MQ=255)
                                                                                    cccggtttCGGCGGGGCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcggcg                                                               >  2:453545/8‑150 (MQ=255)
                                                                                      cTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCa                                                             >  1:1402894/1‑150 (MQ=255)
                                                                                      acggagcggcggCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCCGCCCGGGTGGCCATGCAGCGGCGGCGCa                                                             >  1:65948/7‑150 (MQ=255)
                                                                                       tataTCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAg                                                            <  2:1610246/150‑1 (MQ=255)
                                                                                       tataTCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAg                                                            <  1:1118170/150‑1 (MQ=255)
                                                                                           tCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGc                                                        >  2:1773907/1‑150 (MQ=255)
                                                                                           tCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGc                                                        >  2:1731307/1‑150 (MQ=255)
                                                                                                ggcCGCGCCGACACAGGGGTTCCCCACGGGCGCTTGAGCGGCGGCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTGCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGaa                                                   <  2:409365/150‑1 (MQ=255)
                                                                                                 gcgccggcgAAATCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  <  1:1462437/147‑1 (MQ=255)
                                                                                                 gcCCCGCGGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGGCATGGCCTAAGCCGGTGCGCGCGTGCCGGGCAGCCGTGTACTGGACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  <  2:913162/150‑1 (MQ=255)
                                                                                                 gcCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  <  2:969446/150‑1 (MQ=255)
                                                                                                 gcCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  >  1:1561735/1‑150 (MQ=255)
                                                                                                 gcCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  <  1:1585919/150‑1 (MQ=255)
                                                                                                 gcCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  <  2:282828/150‑1 (MQ=255)
                                                                                                   cccGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGTGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGc                                                <  1:972061/150‑1 (MQ=255)
                                                                                                    ccgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCa                                               >  1:58669/1‑150 (MQ=255)
                                                                                                    ccgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGGGCCGGGCCGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGcc                                               >  2:785422/1‑149 (MQ=255)
                                                                                                     cgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCTGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCaa                                              >  1:805989/1‑150 (MQ=255)
                                                                                                      gccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGcc                                                                                                                 <  1:1793339/82‑1 (MQ=255)
                                                                                                      gccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGcc                                                                                                                 >  2:1793339/1‑82 (MQ=255)
                                                                                                      gccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAg                                             <  1:559129/150‑1 (MQ=255)
                                                                                                      gccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAg                                             <  1:558833/150‑1 (MQ=255)
                                                                                                         gAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGGGTGCTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCgg                                          >  1:118368/1‑150 (MQ=255)
                                                                                                               ggggTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCg                                    <  2:16530/150‑1 (MQ=255)
                                                                                                                 ggTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGc                                  <  2:878311/150‑1 (MQ=255)
                                                                                                                   ttACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGa                                <  1:265637/150‑1 (MQ=255)
                                                                                                                     aCCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCTGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGTCCCTCGGCGAtt                              >  2:1056143/1‑150 (MQ=255)
                                                                                                                       ccAGCTGCTCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGcc                                        >  2:1918595/1‑138 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCAc                                                                                           >  2:1901133/1‑87 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCAc                                                                                           <  1:1901133/87‑1 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGt                                                                                 <  1:151064/97‑1 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGt                                                                                 >  2:151064/1‑97 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAg                                             >  2:139755/1‑133 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAg                                             <  1:139755/133‑1 (MQ=255)
                                                                                                                       ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGcc                                        <  1:1918595/138‑1 (MQ=255)
                                                                                                                           cTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGcc                        <  2:1561735/150‑1 (MQ=255)
                                                                                                                            tGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCg                       <  1:183202/150‑1 (MQ=255)
                                                                                                                                 cTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAAc                  >  1:622192/1‑150 (MQ=255)
                                                                                                                                      cGGCTTCGGCCATGGCACGGCGGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACAt             >  2:136283/1‑150 (MQ=255)
                                                                                                                                       ggCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATg            >  2:585654/1‑150 (MQ=255)
                                                                                                                                         cTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGAt          <  2:1402894/150‑1 (MQ=255)
                                                                                                                                            ccgccgTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCgg       <  2:58669/150‑1 (MQ=255)
                                                                                                                                             cgccgTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGc      <  1:1145548/150‑1 (MQ=255)
                                                                                                                                                 gTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCTTTg  <  2:622192/150‑1 (MQ=255)
                                                                                                                                                    ccATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  >  2:1787536/1‑99 (MQ=255)
                                                                                                                                                    ccATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                  <  1:1787536/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCTTTG  >  NC_002947/2804684‑2804978

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