Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,944 C→T 8.0% A80V (GCG→GTG)  PP_0188 → uroporphyrin‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,9440CT 8.0% 175.3 / 3.9 62A80V (GCG→GTG) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base T (3/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

CCACCCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCT  >  NC_002947/239797‑240091
                                                                                                                                                   |                                                                                                                                                    
ccaccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAAGTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcg                                                                                                                                                     >  5:310405/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                                    >  1:300465/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                                    >  6:291852/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                                    >  6:270548/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGACAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGCGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                                    >  6:256354/1‑149 (MQ=255)
  accccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCt                                                                                                                                                   >  2:268115/1‑149 (MQ=255)
    cccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                 <  2:300465/149‑1 (MQ=255)
            ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                 <  5:178584/141‑1 (MQ=255)
            ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                                 >  6:178584/1‑141 (MQ=255)
             cAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑gcgcg                                                                                                                                       <  5:275962/149‑1 (MQ=255)
             cAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑gcgcg                                                                                                                                       >  6:312571/1‑149 (MQ=255)
              aaGCGCGCCAGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGGGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑gcgcgc                                                                                                                                      <  5:23738/149‑1 (MQ=255)
               aGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCAA‑GCGCGcc                                                                                                                                     >  2:245317/1‑149 (MQ=255)
                  gcTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAg                                                                                                                                  >  5:37251/1‑149 (MQ=255)
                      cGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGgcg                                                                                                                               <  2:225877/148‑1 (MQ=255)
                         cAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCCTCGGCGCGCTgg                                                                                                                           <  3:308483/149‑2 (MQ=255)
                          aGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGTGGCGCTgca                                                                                                                          <  1:150427/149‑1 (MQ=255)
                              gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTgca                                                                                                                          >  7:298041/1‑145 (MQ=255)
                              gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTgca                                                                                                                          <  8:298041/145‑1 (MQ=255)
                              gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCaa                                                                                                                      <  3:130037/149‑1 (MQ=255)
                                   gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                                    <  8:327880/115‑1 (MQ=255)
                                   gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                                    >  7:327880/1‑115 (MQ=255)
                                   gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGt                                                                                                                   <  7:10521/147‑1 (MQ=255)
                                    cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                <  6:37251/149‑1 (MQ=255)
                                    cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                <  5:270548/149‑1 (MQ=255)
                                      gcgcTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGCGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAca                                                                                                              <  1:250575/149‑1 (MQ=255)
                                        gcTGGCCTTGCTGCTTGTTGCGGCCGGGGTGGCGGTTGTTGGCTGGGGTGTCTGGCATGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCTCTGCATCAACGTGAAcagc                                                                                                            <  8:246534/149‑1 (MQ=255)
                                        gcTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAcagc                                                                                                            <  2:42635/149‑1 (MQ=255)
                                            gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                >  3:103472/1‑141 (MQ=255)
                                            gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                                <  4:103472/141‑1 (MQ=255)
                                            gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGa                                                                                                         <  1:245317/148‑1 (MQ=255)
                                               ttGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAg                                                                                                     >  1:317559/1‑149 (MQ=255)
                                                            cGGCCGGGGTTGCGGTAGGTGGCTGGGGTGTCTTGGAGGTGCGTCAACTTCAGGGCAGTGTGTTCAACCAGGGCCAGCACATCGAGGCGGTGACCCATCCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                         <  1:268115/149‑1 (MQ=255)
                                                               ccGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGt                                                                                     >  2:106663/1‑149 (MQ=255)
                                                                         ggTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccg                                                                           >  2:260671/1‑149 (MQ=255)
                                                                         ggTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccg                                                                           >  2:204946/1‑149 (MQ=255)
                                                                             ggTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGCGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCa                                                                        <  5:256354/148‑1 (MQ=255)
                                                                                      ggTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAc                                                                                                               >  4:223683/1‑100 (MQ=255)
                                                                                      ggTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAc                                                                                                               <  3:223683/100‑1 (MQ=255)
                                                                                        tgtCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTagaag                                                            <  2:188914/149‑1 (MQ=255)
                                                                                          tCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAc                                                          >  7:25932/1‑149 (MQ=255)
                                                                                              gCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAccgc                                                       <  5:312571/148‑1 (MQ=255)
                                                                                                aGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCgg                                                    >  1:176170/1‑149 (MQ=255)
                                                                                                     cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                               <  5:291852/149‑1 (MQ=255)
                                                                                                     cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                               >  1:227291/1‑149 (MQ=255)
                                                                                                              cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                                <  5:274942/107‑1 (MQ=255)
                                                                                                              cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                                >  6:274942/1‑107 (MQ=255)
                                                                                                                    cgacaGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCa                                <  3:281398/145‑1 (MQ=255)
                                                                                                                      tGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                              <  5:288992/149‑1 (MQ=255)
                                                                                                                         gTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                            >  7:110901/1‑148 (MQ=255)
                                                                                                                          ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGACAGCGACCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                          >  7:312588/1‑149 (MQ=255)
                                                                                                                            cAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                         <  2:317559/148‑1 (MQ=255)
                                                                                                                            cAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                         <  1:260671/148‑1 (MQ=255)
                                                                                                                              aCCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcag                      >  8:303941/1‑149 (MQ=255)
                                                                                                                               ccAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCATTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGc                                 <  7:221269/137‑1 (MQ=255)
                                                                                                                               ccAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGc                                 >  8:221269/1‑137 (MQ=255)
                                                                                                                                 tgggTCAGCGCCTCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                   >  8:116259/2‑149 (MQ=255)
                                                                                                                                   ggtcAGCGCCTCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGc                                                                  >  2:3058/4‑100 (MQ=255)
                                                                                                                                   ggtcAGCGCCTCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGc                                                                  <  1:3058/97‑1 (MQ=255)
                                                                                                                                    gtcAGCGCCTCGATGTGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGt                >  7:60654/3‑149 (MQ=255)
                                                                                                                                        gCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                             <  2:4419/116‑1 (MQ=255)
                                                                                                                                        gCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                             >  1:4419/1‑116 (MQ=255)
                                                                                                                                          gcctCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          <  7:116259/146‑1 (MQ=255)
                                                                                                                                           caTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          >  3:327343/1‑148 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                              <  3:328617/127‑1 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                              >  4:328617/1‑127 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                   >  4:224596/1‑138 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag                   <  3:224596/138‑1 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc          >  2:103482/1‑147 (MQ=255)
                                                                                                                                                  gcgcTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCt  <  1:204946/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCACCCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCT  >  NC_002947/239797‑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: ‑