Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,949 A→G 6.3% 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 6.3% 181.7 / 2.7 63T82A (ACC→GCC) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/33);  new base G (2/2);  total (28/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 = 9.25e-01

CCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCT  >  NC_002947/239801‑240091
                                                                                                                                                    |                                                                                                                                               
cccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                 <  2:300465/149‑1 (MQ=255)
        ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                 >  6:178584/1‑141 (MQ=255)
        ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                 <  5:178584/141‑1 (MQ=255)
         cAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgcg                                                                                                                                       <  5:275962/149‑1 (MQ=255)
         cAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgcg                                                                                                                                       >  6:312571/1‑149 (MQ=255)
          aaGCGCGCCAGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGGGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAgcgcgc                                                                                                                                      <  5:23738/149‑1 (MQ=255)
           aGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCAAGCGCGcc                                                                                                                                     >  2:245317/1‑149 (MQ=255)
              gcTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAg                                                                                                                                  >  5:37251/1‑149 (MQ=255)
                  cGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGgcg                                                                                                                               <  2:225877/148‑1 (MQ=255)
                     cAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCCTCGGCGCGCTgg                                                                                                                           <  3:308483/149‑2 (MQ=255)
                      aGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGTGGCGCTgca                                                                                                                          <  1:150427/149‑1 (MQ=255)
                          gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgca                                                                                                                          <  8:298041/145‑1 (MQ=255)
                          gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTgca                                                                                                                          >  7:298041/1‑145 (MQ=255)
                          gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCaa                                                                                                                      <  3:130037/149‑1 (MQ=255)
                               gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGt                                                                                                                   <  7:10521/147‑1 (MQ=255)
                                cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                <  6:37251/149‑1 (MQ=255)
                                cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                <  5:270548/149‑1 (MQ=255)
                                  gcgcTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGCGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAca                                                                                                              <  1:250575/149‑1 (MQ=255)
                                    gcTGGCCTTGCTGCTTGTTGCGGCCGGGGTGGCGGTTGTTGGCTGGGGTGTCTGGCATGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCTCTGCATCAACGTGAAcagc                                                                                                            <  8:246534/149‑1 (MQ=255)
                                    gcTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAcagc                                                                                                            <  2:42635/149‑1 (MQ=255)
                                        gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                >  3:103472/1‑141 (MQ=255)
                                        gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                <  4:103472/141‑1 (MQ=255)
                                        gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGa                                                                                                         <  1:245317/148‑1 (MQ=255)
                                           ttGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAg                                                                                                     >  1:317559/1‑149 (MQ=255)
                                                        cGGCCGGGGTTGCGGTAGGTGGCTGGGGTGTCTTGGAGGTGCGTCAACTTCAGGGCAGTGTGTTCAACCAGGGCCAGCACATCGAGGCGGTGACCCATCCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                         <  1:268115/149‑1 (MQ=255)
                                                           ccGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGt                                                                                     >  2:106663/1‑149 (MQ=255)
                                                                     ggTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccg                                                                           >  2:260671/1‑149 (MQ=255)
                                                                     ggTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccg                                                                           >  2:204946/1‑149 (MQ=255)
                                                                         ggTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGCGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCa                                                                        <  5:256354/148‑1 (MQ=255)
                                                                                  ggTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAc                                                                                                               <  3:223683/100‑1 (MQ=255)
                                                                                  ggTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAAc                                                                                                               >  4:223683/1‑100 (MQ=255)
                                                                                    tgtCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTagaag                                                            <  2:188914/149‑1 (MQ=255)
                                                                                      tCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAc                                                          >  7:25932/1‑149 (MQ=255)
                                                                                          gCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAccgc                                                       <  5:312571/148‑1 (MQ=255)
                                                                                            aGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCgg                                                    >  1:176170/1‑149 (MQ=255)
                                                                                                 cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                               >  1:227291/1‑149 (MQ=255)
                                                                                                 cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                               <  5:291852/149‑1 (MQ=255)
                                                                                                          cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                                >  6:274942/1‑107 (MQ=255)
                                                                                                          cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                                <  5:274942/107‑1 (MQ=255)
                                                                                                                cgacaGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCa                                <  3:281398/145‑1 (MQ=255)
                                                                                                                  tGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                              <  5:288992/149‑1 (MQ=255)
                                                                                                                     gTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                            >  7:110901/1‑148 (MQ=255)
                                                                                                                      ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGACAGCGACCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                          >  7:312588/1‑149 (MQ=255)
                                                                                                                        cAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                         <  2:317559/148‑1 (MQ=255)
                                                                                                                        cAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                         <  1:260671/148‑1 (MQ=255)
                                                                                                                          aCCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcag                      >  8:303941/1‑149 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCATTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGc                                 <  7:221269/137‑1 (MQ=255)
                                                                                                                           ccAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGc                                 >  8:221269/1‑137 (MQ=255)
                                                                                                                             tgggTCAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                   >  8:116259/2‑149 (MQ=255)
                                                                                                                               ggtcAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGc                                                                  <  1:3058/97‑1 (MQ=255)
                                                                                                                               ggtcAGCGCCTCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGc                                                                  >  2:3058/4‑100 (MQ=255)
                                                                                                                                gtcAGCGCCTCGATGTGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGt                >  7:60654/3‑149 (MQ=255)
                                                                                                                                    gCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                             >  1:4419/1‑116 (MQ=255)
                                                                                                                                    gCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                             <  2:4419/116‑1 (MQ=255)
                                                                                                                                      gcctCGATGTGCTGGCCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          <  7:116259/146‑1 (MQ=255)
                                                                                                                                       caTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          >  3:327343/1‑148 (MQ=255)
                                                                                                                                        aTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                              >  4:328617/1‑127 (MQ=255)
                                                                                                                                        aTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                              <  3:328617/127‑1 (MQ=255)
                                                                                                                                        aTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                   >  4:224596/1‑138 (MQ=255)
                                                                                                                                        aTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                   <  3:224596/138‑1 (MQ=255)
                                                                                                                                        aTCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          >  2:103482/1‑147 (MQ=255)
                                                                                                                                              gcgcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCt  <  1:204946/149‑1 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGTAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGTCAGCCGCCAGCGAGCTAGAAGACCGCCGGc                                                   <  6:72584/98‑1 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGTAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGTCAGCCGCCAGCGAGCTAGAAGACCGCCGGc                                                   >  5:72584/1‑98 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          <  3:176844/139‑1 (MQ=255)
                                                                                                                                                gcTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          >  4:176844/1‑139 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCAGCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCT  >  NC_002947/239801‑240091

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: