Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,701,215 G→T 7.9% D60E (GAC→GAA)  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,2150GT 7.9% 201.2 / 6.5 63D60E (GAC→GAA) kdpCpotassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base T (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CTGCGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAAC  >  NC_002947/4701071‑4701356
                                                                                                                                                |                                                                                                                                             
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)
                gTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCg                                                                                                                            >  4:114200/1‑148 (MQ=255)
                   gCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcc                                                                                                                        <  5:265832/149‑1 (MQ=255)
                   gCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcc                                                                                                                        >  6:265832/1‑149 (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                                                                                                                      <  2:119640/127‑1 (MQ=255)
                                           gccagcGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAAcccc                                                                                                                      >  1:119640/1‑127 (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                                                                                               >  3:215484/1‑103 (MQ=255)
                                                                                          ggTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGt                                                                                               <  4:215484/103‑1 (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               >  6:168253/1‑147 (MQ=255)
                                                                                                                              ccAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAc               >  8:14141/1‑147 (MQ=255)
                                                                                                                                    aTCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtca       >  5:180336/1‑149 (MQ=255)
                                                                                                                                    aTCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtca       >  8:122756/1‑149 (MQ=255)
                                                                                                                                    aTCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtca       >  5:29329/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)
                                                                                                                                                |                                                                                                                                             
CTGCGCCTGATACAACGTGGCCGCATCGCCCTTTACCCGCTCAGCCAGCGCCGGGTTGCTCGGCGACAGGTTGCTGGCGCCGCTGGCCACGGTCGCGTAGGCACCTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAAC  >  NC_002947/4701071‑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: ‑