Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,701,210 T→A 7.4% Q62L (CAG→CTG)  kdpC ← potassium‑transporting ATPase subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,701,2100TA 7.4% 198.4 / 4.1 67Q62L (CAG→CTG) kdpCpotassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base A (2/3);  total (32/35)
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.23e-01

GCCCCTGCTGCGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAAC  >  NC_002947/4701064‑4701356
                                                                                                                                                  |                                                                                                                                                  
gCCCCTGCTGCGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTg                                                                                                                                                   >  4:129056/1‑148 (MQ=255)
gCCCCTGCTGCGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTg                                                                                                                                                   <  2:64839/148‑1 (MQ=255)
       ctgcGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCGCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCt                                                                                                                                           <  4:127503/149‑1 (MQ=255)
          cgcCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTggg                                                                                                                                        <  7:146533/149‑1 (MQ=255)
          cgcCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTggg                                                                                                                                        <  3:32377/149‑1 (MQ=255)
          cgcCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTgg                                                                                                                                         <  5:324343/148‑1 (MQ=255)
             cTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCg                                                                                                                                      >  6:272812/1‑148 (MQ=255)
              tGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGAt                                                                                                                                    >  7:65984/1‑149 (MQ=255)
                aTACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATc                                                                                                                                   <  5:136054/148‑1 (MQ=255)
                 tACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCAGGACTTCCAGGGCGAtggc                                                                                                                                 <  7:297941/149‑4 (MQ=255)
                    aaCGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGc                                                                                                                              >  8:239477/1‑149 (MQ=255)
                       gTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCTGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTg                                                                                                                                                   <  6:101931/125‑1 (MQ=255)
                       gTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCTGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTg                                                                                                                                                   >  5:101931/1‑125 (MQ=255)
                       gTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCg                                                                                                                            >  4:114200/1‑148 (MQ=255)
                          gCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcc                                                                                                                        >  6:265832/1‑149 (MQ=255)
                          gCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcc                                                                                                                        <  5:265832/149‑1 (MQ=255)
                                    ccTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTgg                                                                                                                                         >  2:312879/1‑122 (MQ=255)
                                    ccTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTgg                                                                                                                                         <  1:312879/122‑1 (MQ=255)
                                          cccGCTCAGCCAGCGCCGGGTTGCTCGGCGACCGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTgc                                                                                                        <  8:179583/149‑1 (MQ=255)
                                           ccGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCtt                                                                                                          <  3:189364/146‑1 (MQ=255)
                                           ccGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCtt                                                                                                          >  4:189364/1‑146 (MQ=255)
                                           ccGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCAGGACTTCCAGGGCGAtg                                                                                                                                   >  2:35284/1‑120 (MQ=255)
                                           ccGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCAGGACTTCCAGGGCGAtg                                                                                                                                   <  1:35284/121‑2 (MQ=255)
                                               tCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATc                                                                                                   >  7:97094/1‑149 (MQ=255)
                                                 agccagcGCAGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCGATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATc                                                                                                   <  4:19793/147‑1 (MQ=255)
                                                  gccagcGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcccc                                                                                                                      >  1:119640/1‑127 (MQ=255)
                                                  gccagcGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcccc                                                                                                                      <  2:119640/127‑1 (MQ=255)
                                                   ccagcGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATcgcg                                                                                                >  8:270280/1‑148 (MQ=255)
                                                    cagcGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCAGGACTTCCAGGGCGAt                                                                                                                                    >  5:205194/1‑111 (MQ=255)
                                                    cagcGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCAGGACTTCCAGGGCGAt                                                                                                                                    <  6:205194/111‑1 (MQ=255)
                                                      gcgcCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGCACCCCGCACCTGGCCTTg                                                                                                         >  7:247088/1‑136 (MQ=255)
                                                      gcgcCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTg                                                                                                         <  8:247088/136‑1 (MQ=255)
                                                          cGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGc                                                                                        <  4:192643/149‑1 (MQ=255)
                                                                 cTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGtt                                                                                 >  6:146545/1‑149 (MQ=255)
                                                                    ggCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTg                                                                                <  5:297701/147‑1 (MQ=255)
                                                                    ggCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGCCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTg                                                                                <  5:272812/147‑1 (MQ=255)
                                                                      cGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCt                                                                            >  1:296505/1‑149 (MQ=255)
                                                                             ttGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCgg                                                                     <  5:146545/149‑1 (MQ=255)
                                                                              tGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCgg                                                                     >  6:73340/1‑148 (MQ=255)
                                                                                   gcgcCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGaaa                                                                 <  5:73340/147‑1 (MQ=255)
                                                                                   gcgcCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGaaa                                                                 <  7:270280/147‑1 (MQ=255)
                                                                                   gcgcCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGc                                                               <  7:239477/149‑1 (MQ=255)
                                                                                       cgcTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCggg                                                                    >  5:312566/1‑140 (MQ=255)
                                                                                       cgcTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCggg                                                                    <  6:312566/140‑1 (MQ=255)
                                                                                       cgcTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCg                                                                      >  6:322971/1‑138 (MQ=255)
                                                                                       cgcTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCg                                                                      <  5:322971/138‑1 (MQ=255)
                                                                                          tGGCCACGGGCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTgg                                                        <  6:169035/149‑1 (MQ=255)
                                                                                            gCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTggg                                                       <  3:129056/148‑1 (MQ=255)
                                                                                            gCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGc                                                      <  2:296505/149‑1 (MQ=255)
                                                                                                 ggTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGt                                                                                               <  4:215484/103‑1 (MQ=255)
                                                                                                 ggTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGt                                                                                               >  3:215484/1‑103 (MQ=255)
                                                                                                             aCCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGc                                       >  6:258220/1‑147 (MQ=255)
                                                                                                                tGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCg                                  >  6:211199/1‑149 (MQ=255)
                                                                                                                  ccGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCgg                                 <  5:285275/148‑1 (MQ=255)
                                                                                                                     aGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATa                              >  8:19530/1‑148 (MQ=255)
                                                                                                                       ggCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAg                           <  3:114200/149‑1 (MQ=255)
                                                                                                                           gcgAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGCACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGgcg                       >  5:36832/1‑149 (MQ=255)
                                                                                                                                   aaCCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGgcgc                      >  5:143135/1‑142 (MQ=255)
                                                                                                                                   aaCCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGgcgc                      <  6:143135/142‑1 (MQ=255)
                                                                                                                                     ccAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAc               >  8:14141/1‑147 (MQ=255)
                                                                                                                                     ccAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAc               >  6:168253/1‑147 (MQ=255)
                                                                                                                                           aTCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtca       >  5:180336/1‑149 (MQ=255)
                                                                                                                                           aTCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtca       >  5:29329/1‑149 (MQ=255)
                                                                                                                                           aTCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtca       >  8:122756/1‑149 (MQ=255)
                                                                                                                                            tCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtcat      <  7:122756/149‑1 (MQ=255)
                                                                                                                                                ccTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATcaac  <  6:36832/149‑1 (MQ=255)
                                                                                                                                                ccTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATcaa   >  7:293497/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCCCTGCTGCGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAAC  >  NC_002947/4701064‑4701356

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: