Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,949 A→G 7.6% T82A (ACC→GCC)  PP_0188 → uroporphyrin‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,9490AG 7.6% 186.4 / 3.6 65T82A (ACC→GCC) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAA  >  NC_002947/239801‑240096
                                                                                                                                                    |                                                                                                                                                   
cccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                     >  4:306639/1‑149 (MQ=255)
  caccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAcc                                                                                                                                                   <  3:306639/149‑1 (MQ=255)
   accTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAccc                                                                                                                                                  >  5:27553/1‑149 (MQ=255)
    ccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCa                                                                                                                                                 <  2:183808/149‑1 (MQ=255)
          aaGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGGGTCAgcgcct                                                                                                                                           <  5:310640/149‑3 (MQ=255)
               cTCCGGCAGTGGCCTGGCAGCGCTTGCCTTGCTGCTTGGTGCGTCCGGGGTGGCGGTAGGTGGCTGTGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCACTGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAgg                                                                                                                                      <  1:61387/149‑1 (MQ=255)
                tCCGGCAGTGGCCTGGCAGCTCTGCCCTTGCTGCATGGTGCGGCCCGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGAAGTTAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAAgc                                                                                                                                     >  5:219965/1‑149 (MQ=255)
                     gaGTGGCCTGGCATCGGGGGGCTTGGTGCTTGGTGCGGCCGGGGTTGCTGTAGGTGGGTGGGGTGTCTGGCAGGTGCGTCATCTGCTGGGCAGTGAGTTCAACCTTGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGGGCTgc                                                                                                                                <  6:219965/148‑1 (MQ=255)
                      aGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgca                                                                                                                               <  2:312292/149‑1 (MQ=255)
                         ggCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgcagc                                                                                                                             >  4:281642/1‑148 (MQ=255)
                           ccTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAAc                                                                                                                          <  3:42104/149‑1 (MQ=255)
                           ccTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCCTCGGCGCGCTgggtcagc                                                                                                                          <  1:116570/149‑8 (MQ=255)
                            cTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGGTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACg                                                                                                                         <  5:94841/149‑1 (MQ=255)
                                cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                     <  1:140661/149‑1 (MQ=255)
                                 aGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                     >  3:34554/1‑148 (MQ=255)
                                      tGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAcagcag                                                                                                               <  7:190527/149‑1 (MQ=255)
                                          cTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCa                                                                                                           <  6:27553/149‑1 (MQ=255)
                                                gctTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTgcg                                                                                                      >  6:37981/1‑148 (MQ=255)
                                                   tGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAAc                                                                                                  >  2:337422/1‑149 (MQ=255)
                                                     gTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTg                                                                                                <  3:281642/149‑1 (MQ=255)
                                                      tGCGGCCGGGGTGGCGGTAGGTGGCTGGTGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTgg                                                                                               <  1:285712/149‑1 (MQ=255)
                                                         ggCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCa                                                                                            <  1:78704/149‑1 (MQ=255)
                                                         ggCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCa                                                                                            <  5:211028/149‑1 (MQ=255)
                                                              gggTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGt                                                                                                         <  7:200878/131‑1 (MQ=255)
                                                              gggTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGt                                                                                                         >  8:200878/1‑131 (MQ=255)
                                                                   gCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGc                                                                                      >  4:78257/1‑145 (MQ=255)
                                                                   gCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGc                                                                                      <  3:78257/145‑1 (MQ=255)
                                                                           tGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGc                                                                           >  4:315194/1‑148 (MQ=255)
                                                                            gggtGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGc                                                                           >  8:263894/4‑147 (MQ=255)
                                                                            ggCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGACTGCCAGCCGCCAGCGa                                                                         >  6:248154/1‑149 (MQ=255)
                                                                                    tgtCTGGCGGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTagaag                                                                 <  5:307054/149‑1 (MQ=255)
                                                                                    tgtCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTag                                                                    >  3:154905/1‑146 (MQ=255)
                                                                                            aGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAg                                                                                                                                       <  6:116748/71‑1 (MQ=255)
                                                                                            aGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAg                                                                                                                                       >  5:116748/1‑71 (MQ=255)
                                                                                            aGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCgg                                                         >  2:216588/1‑149 (MQ=255)
                                                                                              gTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAccc                                                                                                                                                  >  8:58656/1‑58 (MQ=255)
                                                                                              gTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGAccc                                                                                                                                                  <  7:58656/58‑1 (MQ=255)
                                                                                              gTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCg                                                       >  7:308350/1‑149 (MQ=255)
                                                                                                gCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGCGGCGCGGCAGCAACGGGCCCAGCAGCGCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCGAGAAGACCGCCGGCggt                                                     >  8:285770/1‑147 (MQ=255)
                                                                                                 cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGGAGCAACGTGATCAGCAGATCAGTGCGCAACTTGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                                    <  5:248154/149‑1 (MQ=255)
                                                                                                 cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                                    <  5:37981/149‑1 (MQ=255)
                                                                                                   tCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGt                                                                                                         >  5:235607/1‑94 (MQ=255)
                                                                                                   tCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGt                                                                                                         <  6:235607/94‑1 (MQ=255)
                                                                                                       cTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGc                                              >  3:158492/1‑149 (MQ=255)
                                                                                                       cTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGc                                              <  1:337422/149‑1 (MQ=255)
                                                                                                       cTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGc                                              >  4:273072/1‑149 (MQ=255)
                                                                                                        tGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGcc                                             >  8:168594/1‑149 (MQ=255)
                                                                                                          cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGccc                                            <  3:273072/148‑1 (MQ=255)
                                                                                                             ggCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGtt                                        <  3:315194/149‑1 (MQ=255)
                                                                                                               cAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTa                                                >  1:318834/1‑139 (MQ=255)
                                                                                                               cAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTa                                                <  2:318834/139‑1 (MQ=255)
                                                                                                                  tGAGTTCAACTAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                                   <  6:258024/149‑1 (MQ=255)
                                                                                                                   gAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAggg                                  >  2:35434/1‑149 (MQ=255)
                                                                                                                      ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCCGCGAGCTAGAAGCCCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                               >  6:92710/1‑149 (MQ=255)
                                                                                                                        cAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                              <  4:34554/148‑1 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCa                                           >  6:308992/1‑132 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCa                                           <  5:308992/132‑1 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagc                          <  7:168594/149‑1 (MQ=255)
                                                                                                                             tgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGATGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                                 <  1:221687/139‑1 (MQ=255)
                                                                                                                             tgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                                 >  2:221687/2‑140 (MQ=255)
                                                                                                                             tgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                        >  1:87768/2‑149 (MQ=255)
                                                                                                                              gggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCagcagc                       >  3:4041/1‑149 (MQ=255)
                                                                                                                               ggCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCg                      <  1:229220/149‑1 (MQ=255)
                                                                                                                                             ggCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGcc        >  7:27033/1‑149 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTg      >  7:11530/1‑148 (MQ=255)
                                                                                                                                                   gACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGaaa  >  4:17856/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAA  >  NC_002947/239801‑240096

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: ‑