Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,356,043 C→T 8.3% intergenic (+310/‑82) PP_1180 → / → PP_1181 hypothetical protein/two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,356,0430CT 8.3% 172.9 / 5.0 60intergenic (+310/‑82)PP_1180/PP_1181hypothetical protein/two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (2/3);  total (28/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 = 7.82e-01

GTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCTCGCCGAACTCACCGCC  >  NC_002947/1355895‑1356190
                                                                                                                                                    |                                                                                                                                                     
gTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGgc                                                                                                                                                       <  8:323862/149‑1 (MQ=255)
    tttGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGTGCACCGCCCGCGGGGCGCTc                                                                                                                                                   <  6:384858/149‑1 (MQ=255)
      tGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGTCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCgg                                                                                                                                                 >  2:299183/1‑149 (MQ=255)
              ttttCTGCACCGTGACGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTGCAGTCGAGGGGGCGGTAGGGCGTGGCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGATCTCCCAg                                                                                                                                         <  5:377197/149‑1 (MQ=25)
              ttttCAGCACCGTGATGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGTGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGATCTCCCa                                                                                                                                          <  8:231527/148‑1 (MQ=25)
                 tCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCGGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAACATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAggg                                                                                                                                       <  4:344125/148‑1 (MQ=255)
                     cACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGcccc                                                                                                                                   <  4:116068/148‑1 (MQ=255)
                      aCCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGccccc                                                                                                                                  >  3:177538/1‑148 (MQ=255)
                        cGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCAt                                                                                                                              <  5:362375/149‑1 (MQ=255)
                         gTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATc                                                                                                                             >  8:242742/1‑149 (MQ=255)
                          tGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATc                                                                                                                             >  1:37475/1‑148 (MQ=255)
                          tGAAGCCCTTGGGCAGTTTATTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCt                                                                                                                            <  1:124080/149‑1 (MQ=255)
                            aaGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTa                                                                                                                           <  7:242742/148‑1 (MQ=255)
                                ccTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGAtc                                                                                                                                               >  4:63316/1‑125 (MQ=37)
                                ccTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGAtc                                                                                                                                               <  3:63316/125‑1 (MQ=37)
                                 cTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCa                                                                                                                     >  7:80361/1‑149 (MQ=255)
                                 cTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCa                                                                                                                     >  3:251926/1‑149 (MQ=255)
                                  ttGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCa                                                                                                                     >  3:239972/1‑148 (MQ=255)
                                       cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                                                    >  2:347445/1‑113 (MQ=255)
                                       cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                                                    <  1:347445/113‑1 (MQ=255)
                                       cAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGGCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATActct                                                                                                              <  2:279773/149‑1 (MQ=255)
                                        aGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACtctc                                                                                                             <  6:136018/149‑1 (MQ=255)
                                         gTTTTTTTTTGGCGCGAAGCCTGTCACTAAAAACGGTTCAGTCTAGGGGGCGGTAGGGCTTGTCTGAAGATCTGGTGTGCCTGGGATATCGAGCACCGCCCTCGGGGCGCTCGGTCTAGCAGGGCCCCCCATCTACATTAATACTCTCa                                                                                                            <  5:212591/149‑1 (MQ=255)
                                          ttttttttGGGCGCGAAGCCTTTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                                                    <  4:24997/110‑1 (MQ=255)
                                          ttttttttGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCt                                                                                                                                                    >  3:24997/1‑110 (MQ=255)
                                                   ggCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGa                                                                                                  <  8:87840/149‑1 (MQ=255)
                                                           gCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATcc                                                                                          <  2:37475/149‑1 (MQ=255)
                                                            ccTGTCAGAAATAAGGGTTCTGTCGAGGGGGCGGTAGGGCGTGTCTGATGATCTGGGGTGCCTGTGATATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCTTACTCTCAGCTGCCCTGACACGATccc                                                                                         <  2:250264/149‑1 (MQ=255)
                                                             cTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATccc                                                                                         >  4:245758/1‑148 (MQ=255)
                                                               gTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCaa                                                                                       >  7:328565/1‑148 (MQ=255)
                                                                  acaAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGTTGCCTGTTATATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGaa                                                                                    <  1:299183/148‑1 (MQ=255)
                                                                      aaaaCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACtctc                                                                                                             >  5:346115/1‑119 (MQ=255)
                                                                      aaaaCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACtctc                                                                                                             <  6:346115/119‑1 (MQ=255)
                                                                        aaCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGGCTCGCAGGGCCCCCCACTCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAAcaacaa                                                                              <  1:274341/149‑1 (MQ=255)
                                                                         aCGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAgg                                                                            <  4:251926/149‑1 (MQ=255)
                                                                                    gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                                                 >  8:231204/1‑149 (MQ=255)
                                                                                    gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                                                 <  8:98378/149‑1 (MQ=255)
                                                                                    gAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                                                 <  2:36958/149‑1 (MQ=255)
                                                                                            ggTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCAccc                                                         >  2:358966/1‑149 (MQ=255)
                                                                                                 ggCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGac                                                              >  3:182290/1‑139 (MQ=255)
                                                                                                 ggCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGac                                                              <  4:182290/139‑1 (MQ=255)
                                                                                                       tCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGc                                               <  4:239972/148‑1 (MQ=255)
                                                                                                         tGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATc                                            >  6:57451/1‑149 (MQ=255)
                                                                                                         tGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATc                                            >  7:27796/1‑149 (MQ=255)
                                                                                                          gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCAGTACTCTCAGCTGCCCTGACACGATCCCAAGa                                                                                     <  1:39486/108‑1 (MQ=255)
                                                                                                          gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGaacaac                                                                                >  6:246860/1‑112 (MQ=255)
                                                                                                          gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGaacaac                                                                                <  5:246860/112‑1 (MQ=255)
                                                                                                          gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATg                                                                 <  6:339690/127‑1 (MQ=255)
                                                                                                          gAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATcc                                           >  5:302744/1‑149 (MQ=255)
                                                                                                          gAAGATCTGGGGTGCCTCGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACAACCCCATACCCCGCATc                                            >  4:149791/1‑148 (MQ=255)
                                                                                                           aaGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCt                                          >  7:115443/1‑149 (MQ=255)
                                                                                                            aGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCt                                          >  6:19839/1‑148 (MQ=255)
                                                                                                             gATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTc                                         >  8:31017/1‑148 (MQ=255)
                                                                                                               tCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCCTCCTCAtc                                      >  4:175962/1‑149 (MQ=255)
                                                                                                                    gggcccTGCGAGACCGAGCGCCCCGCGGGCGGTGCTCGATCTCCCAGGCACCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCAtcgtcg                                  >  6:312170/5‑147 (MQ=21)
                                                                                                                    ggTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGa                                 <  1:358966/149‑1 (MQ=255)
                                                                                                                          tGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAgacga                           >  2:314589/1‑149 (MQ=255)
                                                                                                                            ggAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAgacgac                          <  7:31017/148‑1 (MQ=255)
                                                                                                                                   gAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCt                  <  8:115443/149‑1 (MQ=255)
                                                                                                                                    aGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCTc                 <  3:149791/149‑1 (MQ=255)
                                                                                                                                         cgcccgcGGGGCGCTCGGTCTCGCAGGGCCCCCCATCGACATTCATACTCTCCGCTGCCCTGACACGATCCCAAGACCAACAAGGGCACAGCCATGGACCCCCCCATACCCCGCCTCCTCACCGTCGAAGACGACCAGCGCCTCGCCg             >  7:270901/1‑148 (MQ=255)
                                                                                                                                                   gcgcTCGGTCTCGCAGGGCCCGCCATCTACATTCATACTCTCAGCAGCCCTGACACGAGACCAAGAACAACAAGGCCAACGCCATGGACCAGCCCAAACACCGCAGCCTCAACGGCGAAGACGACCACCGCCTCGCCGCACTCAccgcc  >  4:41268/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTTTGCAGTATTTTTCTGCACCGTGAAGCCCTTGGGCAGTTTTTTTTGGGCGCGAAGCCTGTCACAAAAAACGGTTCAGTCGAGGGGGCGGTAGGGCGTGTCTGAAGATCTGGGGTGCCTGGGAGATCGAGCACCGCCCGCGGGGCGCTCGGTCTCGCAGGGCCCCCCATCTACATTCATACTCTCAGCTGCCCTGACACGATCCCAAGAACAACAAGGCCAAAGCCATGGACCACCCCATACCCCGCATCCTCATCGTCGAAGACGACCAGCGCCTCGCCGAACTCACCGCC  >  NC_002947/1355895‑1356190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: