Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,324,162 C→A 6.0% R266R (CGG→AGG)  PP_4684 → kinase‑like domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,324,1620CA 6.0% 214.5 / 3.9 67R266R (CGG→AGG) PP_4684kinase‑like domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base A (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

CCTGGGCAACGCCACCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGATGGCGTTCAGCGTGC  >  NC_002947/5324017‑5324307
                                                                                                                                                  |                                                                                                                                                 
ccTGGGCAACGCCACCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGgc                                                                                                                                                <  4:51969/149‑1 (MQ=255)
       aaCGCCACCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGCGCCGGGCCAGGCACTtcag                                                                                                                                         <  7:73715/149‑2 (MQ=255)
         cGCCACCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtcatca                                                                                                                                       <  5:288635/149‑1 (MQ=255)
             aCCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGc                                                                                                                                     >  8:149944/1‑147 (MQ=255)
             aCCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCa                                                                                                                                   <  3:345157/149‑1 (MQ=255)
                    tCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGCGCCGGGCCAGGCACTtcagg                                                                                                                                        <  3:4389/137‑3 (MQ=255)
                    tCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGCGCCGGGCCAGGCACTtcagg                                                                                                                                        >  4:4389/1‑135 (MQ=255)
                                tggtgATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACgggg                                                                                                                <  7:207120/149‑1 (MQ=255)
                                  gtgATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGAGCCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACgggg                                                                                                                >  1:67706/1‑148 (MQ=255)
                                   tgATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGAt                                                                                                             <  6:261326/149‑1 (MQ=255)
                                   tgATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGAt                                                                                                             >  1:193826/1‑149 (MQ=255)
                                    gATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTTGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGAt                                                                                                             >  4:163816/1‑148 (MQ=255)
                                      tCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTa                                                                                                           <  3:163816/148‑1 (MQ=255)
                                        ttCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCTCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTAc                                                                                                          <  1:62176/147‑1 (MQ=255)
                                        ttCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTAc                                                                                                          >  2:62176/1‑147 (MQ=255)
                                          cGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTATGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTAc                                                                                                          <  1:146049/145‑1 (MQ=255)
                                           gATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTTGAACTGACGGGGGATTACGAAg                                                                                                      <  7:220414/148‑1 (MQ=255)
                                             ttGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGACCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGaa                                                                                                       >  7:19326/1‑145 (MQ=255)
                                             ttGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGACCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGaa                                                                                                       <  8:19326/145‑1 (MQ=255)
                                              tGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCt                                                                                                  >  4:315882/1‑149 (MQ=255)
                                                 aTCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           <  3:173275/105‑1 (MQ=255)
                                                 aTCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           >  4:173275/1‑105 (MQ=255)
                                                 aTCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTgg                                                                                                >  8:85745/1‑148 (MQ=255)
                                                  tCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGaa                                                                                              >  6:204165/1‑149 (MQ=255)
                                                  tCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGaa                                                                                              >  6:333052/1‑149 (MQ=255)
                                                  tCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGa                                                                                               >  8:61147/1‑148 (MQ=255)
                                                        tCAACGAGCCTTTACGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCt                                                                                         >  3:179109/1‑148 (MQ=255)
                                                         cAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTc                                                                                        <  4:179109/148‑1 (MQ=255)
                                                         cAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGAATACGAAGGCCTGGAACTGCTc                                                                                        <  6:219287/148‑1 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           <  4:338806/96‑1 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           >  5:248183/1‑96 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           >  5:252588/1‑96 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           <  6:252588/96‑1 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           >  3:338806/1‑96 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTtc                                                                                                                                           <  6:248183/96‑1 (MQ=255)
                                                          aaCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAAc                                                                                     >  5:148173/1‑149 (MQ=255)
                                                           aCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAg                                                                                                      >  2:234653/1‑131 (MQ=255)
                                                           aCGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAg                                                                                                      <  1:234653/131‑1 (MQ=255)
                                                             gAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACtt                                                                                   <  2:67706/149‑1 (MQ=255)
                                                                ccTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTa                                                                                <  7:61147/149‑1 (MQ=255)
                                                                ccTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTa                                                                                <  3:315882/149‑1 (MQ=255)
                                                                 cTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTAc                                                                               <  7:149944/149‑1 (MQ=255)
                                                                      cGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGc                                                                          <  3:211496/149‑1 (MQ=255)
                                                                         cTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTa                                                                       <  5:204165/149‑1 (MQ=255)
                                                                            aCCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTAc                                                                      <  2:193826/147‑1 (MQ=255)
                                                                                          ccGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACgg                                                       <  5:333052/149‑1 (MQ=255)
                                                                                          ccGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGAAGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCTTACgg                                                       <  5:249699/149‑1 (MQ=255)
                                                                                                 cGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGACGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAgg                                                >  8:156783/1‑149 (MQ=255)
                                                                                                 cGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAAc                                                                                             >  8:310433/1‑104 (MQ=255)
                                                                                                 cGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAAc                                                                                             <  7:310433/104‑1 (MQ=255)
                                                                                                 cGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAgg                                                >  7:173506/1‑149 (MQ=255)
                                                                                                           gCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGtt                                      <  7:156783/149‑1 (MQ=255)
                                                                                                           gCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGt                                       <  6:76187/148‑1 (MQ=255)
                                                                                                            ccATGGACCTGGAAGACCGCGGCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCa                                    <  6:148173/149‑1 (MQ=255)
                                                                                                               tGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGc                                  <  8:173506/149‑1 (MQ=255)
                                                                                                                       gAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCAtt                                         >  1:153704/1‑134 (MQ=255)
                                                                                                                       gAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCAtt                                         <  2:153704/134‑1 (MQ=255)
                                                                                                                       gAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATgccg                            <  7:85745/147‑1 (MQ=255)
                                                                                                                          gACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCa                                    <  7:60212/136‑1 (MQ=255)
                                                                                                                          gACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCa                                    >  8:60212/1‑136 (MQ=255)
                                                                                                                                  gatGAAGCGCCGGGCCAGGCACTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCGTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGATg               >  3:14924/3‑149 (MQ=255)
                                                                                                                                      aaGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTAAACTTCTACAAAGCCTACCGCGCGCTCGTACTCGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGAAGGCGt           >  4:223539/1‑149 (MQ=255)
                                                                                                                                           ccTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGc                                                                          <  5:15917/81‑1 (MQ=255)
                                                                                                                                           ccTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGc                                                                          >  6:15917/1‑81 (MQ=255)
                                                                                                                                             tGGCCCGGCGCTTCGTCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGATGGCGTTCAGCGt    >  2:88271/1‑149 (MQ=255)
                                                                                                                                               gCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGATGGCGTTCAGCGt    <  3:223539/147‑1 (MQ=255)
                                                                                                                                               gCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGATGGCGTTCAGCGTGc  >  5:179129/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCTGGGCAACGCCACCTTGATCGATGGCAAGGTGGTGATCTTCGATTGCATCGAGTTCAACGAGCCTTTCCGACTGACCGACGTCTACGCCGACACCGCCTTCCTGGCCATGGACCTGGAAGACCGCGGCCTGAAGTGCCTGGCCCGGCGCTTCATCAGCCAGTACCTGGAACTGACGGGGGATTACGAAGGCCTGGAACTGCTCAACTTCTACAAAGCCTACCGCGCGCTCGTACGGGCCAAGGTCGCATTGTTCAGCATGCCGGCCGACGCTGATGGCGTTCAGCGTGC  >  NC_002947/5324017‑5324307

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: