Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,027,829 C→A 9.7% A201S (GCC→TCC)  PP_3552 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,027,8290CA 9.7% 217.9 / 4.9 74A201S (GCC→TCC) PP_3552two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base A (6/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATGTACAGGCATG  >  NC_002947/4027685‑4027971
                                                                                                                                                  |                                                                                                                                              
tGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGAtt                                                                                                                                            >  3:377300/1‑149 (MQ=255)
  cccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGAttt                                                                                                                                           <  6:263067/148‑1 (MQ=255)
  cccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGAttt                                                                                                                                           <  8:308328/148‑1 (MQ=255)
   ccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCg                                                                                                                                         <  8:230652/149‑1 (MQ=255)
     gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGaa                                                                                                                                       <  3:356033/149‑1 (MQ=255)
     gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGaa                                                                                                                                       <  5:40160/149‑1 (MQ=255)
     gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGaa                                                                                                                                       >  2:244388/1‑149 (MQ=255)
         ccAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCaa                                                                                                                                   >  6:104044/1‑149 (MQ=255)
           aGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGATCTAAATTTAGAACCAACg                                                                                                                                 >  2:69795/1‑149 (MQ=255)
              tCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                                                                                                >  5:4471/1‑147 (MQ=255)
                gCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCg                                                                                                                                         >  2:144157/1‑136 (MQ=255)
                gCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCg                                                                                                                                         <  1:144157/136‑1 (MQ=255)
                                 gggCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCg                                                                                                            >  8:365877/1‑148 (MQ=255)
                                  ggCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTa                                                                                                          >  6:342549/1‑149 (MQ=255)
                                    cACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGCCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTAcc                                                                                                        >  4:35836/1‑149 (MQ=255)
                                           cgcCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCCTACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCa                                                                                                   <  1:325475/149‑1 (MQ=255)
                                             ccAGGCTTTCCTGCGCTTCGATGTGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGt                                                                                                >  3:398197/1‑148 (MQ=255)
                                              cAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGTACACGCTGACCACCTCAACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCa                                                                                              >  1:400278/1‑149 (MQ=255)
                                                ggCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc                                                                                            <  5:255031/149‑1 (MQ=255)
                                                 gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCg                                                                                           >  1:378814/1‑149 (MQ=255)
                                                 gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGATCTC‑‑ACTGATCGCCCTAAAATAAAAAACAAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCg                                                                                           >  6:364674/1‑149 (MQ=255)
                                                     tCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCgtg                                                                                         <  1:244388/147‑1 (MQ=255)
                                                     tCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGACAGGTCATCgtgt                                                                                        <  6:217258/148‑1 (MQ=255)
                                                      ccTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg                                                                                       <  7:361895/148‑1 (MQ=255)
                                                       cTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTgg                                                                                      <  7:365877/148‑1 (MQ=255)
                                                             ttCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                                                                                                         <  4:238811/107‑1 (MQ=255)
                                                             ttCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                                                                                                         >  3:238811/1‑107 (MQ=255)
                                                             ttCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGt                                                                               >  5:67602/1‑149 (MQ=255)
                                                               cGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAAGGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTcc                                                                             >  1:331580/1‑149 (MQ=255)
                                                                gATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTcc                                                                             <  8:162005/148‑1 (MQ=255)
                                                                gATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTcc                                                                             <  7:134071/148‑1 (MQ=255)
                                                                 aTGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGt                                                                                                <  1:277399/128‑1 (MQ=255)
                                                                 aTGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGt                                                                                                >  2:277399/1‑128 (MQ=255)
                                                                 aTGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTg                                                                           <  2:378814/149‑1 (MQ=255)
                                                                    cgcAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTc                                                                        <  6:132279/149‑1 (MQ=255)
                                                                      cAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCa                                                                      <  6:67602/149‑1 (MQ=255)
                                                                         ggTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCTCCTGGCGTACCCGGCAGGTCATc                                                                                            >  7:291538/1‑124 (MQ=255)
                                                                         ggTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCTCCTGGCGTACCCGGCAGGTCATc                                                                                            <  8:291538/124‑1 (MQ=255)
                                                                            gATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                                                                                                         >  3:137190/1‑92 (MQ=255)
                                                                            gATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                                                                                                         <  4:137190/92‑1 (MQ=255)
                                                                               gTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAACCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAg                                                             >  7:319665/1‑149 (MQ=255)
                                                                                  gcgGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCgg                                                          >  8:95654/1‑149 (MQ=255)
                                                                                  gcgGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGATCTC‑‑AATGATCGCCCAAAAATAAAAAAAAAAAAAACGGTAGCCTCCTGTCGCCTGGCCTACACCGCAGGTCATCGTGAGGTCTCCGTCCTGGTCCCGCTCGGAAGCgg                                                          >  4:273349/1‑149 (MQ=255)
                                                                                   cgGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGc                                                         <  5:364674/149‑1 (MQ=255)
                                                                                         acGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGgcg                                                   <  3:35836/149‑1 (MQ=255)
                                                                                           gCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCa                                                 <  5:342549/149‑1 (MQ=255)
                                                                                           gCTGACCACCTCCACCACCGCACCGGTGTAGATCTC‑‑ACTGATCGCCCTAAAATAAAAAAAAAAAAAACGGTAGACTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGgcg                                                   >  6:363123/1‑147 (MQ=255)
                                                                                              gACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTa                                              <  3:140908/149‑1 (MQ=255)
                                                                                                 taccTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGAGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCt                                           >  4:112573/2‑149 (MQ=255)
                                                                                                  accTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCt                                           >  3:285951/1‑148 (MQ=255)
                                                                                                    cTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGAc                                        >  2:269866/1‑149 (MQ=255)
                                                                                                      ccaccaCCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCgg                                                          <  4:54583/129‑1 (MQ=255)
                                                                                                      ccaccaCCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCgg                                                          >  3:54583/1‑129 (MQ=255)
                                                                                                      ccaccaCCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGAcgc                                      >  7:201626/1‑149 (MQ=255)
                                                                                                         ccaccGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGAcgccgc                                   <  2:331580/149‑1 (MQ=255)
                                                                                                           accgcaccgGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTc                                 <  4:377300/149‑1 (MQ=255)
                                                                                                            ccgcaccgGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTcg                                >  3:176539/1‑149 (MQ=255)
                                                                                                              gcaccgGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTcgcg                              <  1:238700/149‑1 (MQ=255)
                                                                                                                 gcgGTGTAGGTCTC‑‑ACTGATCGCCCGGCTGGAGATTTCGAACGAACGGTAGCCTCCTGTCGCCTCGCGTTCCCGGCAGTTCATCGTGGTGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCgg                           <  2:400278/148‑1 (MQ=255)
                                                                                                                  cgGTGTAGGTCTC‑‑ACGGATCGCCCTAATATATATAAAAAAAAAACGTTAGCATCCTCTCGCCTGGCGTACCCGGCAGCTCATCGTGAGGAATCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGa                          >  6:19742/1‑149 (MQ=255)
                                                                                                                  cgGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGa                          <  1:269866/149‑1 (MQ=255)
                                                                                                                   ggTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGAt                         <  8:201626/149‑1 (MQ=255)
                                                                                                                    gtgtAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATg                        >  5:204871/1‑149 (MQ=255)
                                                                                                                           tctc‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAgg                 >  7:170142/1‑149 (MQ=255)
                                                                                                                              c‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGgcgc                                                  >  7:164144/1‑113 (MQ=255)
                                                                                                                              c‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGgcgc                                                  <  8:164144/113‑1 (MQ=255)
                                                                                                                                 aCGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAgg                 >  1:190628/1‑145 (MQ=255)
                                                                                                                                 aCGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAgg                 <  2:190628/145‑1 (MQ=255)
                                                                                                                                   ggATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATGt           <  8:170142/149‑1 (MQ=255)
                                                                                                                                   aaaTCGCCAGCCGGGCGATTTCGAACCAACGGTAGCCTCCTGCCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATg            >  3:171466/3‑148 (MQ=255)
                                                                                                                                        gCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCTTCGTGTGGTAGCCGTCCTGGTCCAGTGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATGTACAg       <  7:232188/148‑1 (MQ=255)
                                                                                                                                         cccGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATGTACAgg      >  8:387241/1‑148 (MQ=255)
                                                                                                                                          ccGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATGTACAGGCa    >  7:185989/1‑149 (MQ=255)
                                                                                                                                            ttcTGGAGATTTCTAAGCAACGGTAGCCTCCTGTGGCCTGGGGTACCCGGCAGGTGATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGGCTGGCGGAGTACATGTCGACGCTCGCGCGGATGCAGGAGGGTATGTACAGGCATg  <  4:398197/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTC‑‑ACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCCTGGTCCAGCGCGGCAGCGGCCTGGCGCAGTACCTGACGCCGCTCGCGCGGATGCAGCAGGGTATGTACAGGCATG  >  NC_002947/4027685‑4027971

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: