Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,249,446 G→A 6.9% R126C (CGT→TGT)  PP_4624 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,249,4460GA 6.9% 172.2 / 2.2 58R126C (CGT→TGT) PP_4624alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base A (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATTACGCACCAGGTCGTCACGCCCGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCAGGCTGTAAGTGTGAT  >  NC_002947/5249301‑5249589
                                                                                                                                                 |                                                                                                                                                
acaTTACGCACCAGGTCGTCACGCCCGAACGACCGCCCGACCAGCCAGCCAAAGCGGCTGCGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACTCGCAACTGGCGATGTCGGCCCgctg                                                                                                                                               >  3:255228/1‑149 (MQ=255)
  aTTACGCACCAGGTCGTCACGCCCGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACgctgct                                                                                                                                             <  4:196249/149‑1 (MQ=255)
            aGGTCGTCACGCCCGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACACTGCTCATGGTGTcg                                                                                                                                   >  6:190671/1‑149 (MQ=255)
             ggTCGTCACGCCCGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTcgc                                                                                                                                  <  5:174207/149‑1 (MQ=255)
               tcgtcACGCCCGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCgcgc                                                                                                                                <  2:22501/149‑1 (MQ=255)
                        ccGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGCGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTc                                                                                                                       >  7:126400/1‑149 (MQ=255)
                          gAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCt                                                                                                                     <  4:313890/149‑1 (MQ=255)
                                gcgcCCGACCAGCCAGCC‑AAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACACTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCAc                                                                                                              >  2:38343/1‑149 (MQ=255)
                                     ctaCCAGCCAGCCAACGCGGCTGAGCAGCAGCTTCTGGATCTGTAGCAAGCGGCAGCTTCCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTg                                                                                                          <  5:152490/147‑1 (MQ=255)
                                     cGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCAt                                                                                                                                          <  1:154833/117‑1 (MQ=255)
                                     cGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCAt                                                                                                                                          >  2:154833/1‑117 (MQ=255)
                                       aCCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGTTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTgcc                                                                                                        >  1:158914/1‑149 (MQ=255)
                                        ccagccagCC‑AAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACACTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTgccgc                                                                                                      <  1:38343/149‑1 (MQ=255)
                                               gccaAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGt                                                                                               >  2:241795/1‑149 (MQ=255)
                                               gccaAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGt                                                                                               >  7:343939/1‑149 (MQ=255)
                                                  aaaGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGt                                                                                            >  5:156610/1‑149 (MQ=255)
                                                   aaGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTg                                                                                           >  3:345412/1‑149 (MQ=255)
                                                     gCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTggg                                                                                         >  5:244705/1‑149 (MQ=255)
                                                       ggCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGc                                                                                        >  8:321368/1‑148 (MQ=255)
                                                          tGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCcag                                                                                     >  8:4781/1‑148 (MQ=255)
                                                          tGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAACCACGCTGCCGCC‑GTAGTCGTGGGCcag                                                                                     >  4:350300/1‑148 (MQ=255)
                                                            agcagcagcTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCcagc                                                                                    <  4:245183/147‑1 (MQ=255)
                                                            agcagcagcTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCcagc                                                                                    <  7:321368/147‑1 (MQ=255)
                                                            agcagcagcTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCcagc                                                                                    <  4:250037/147‑1 (MQ=255)
                                                                    cTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACgg                                                                          >  1:119480/1‑149 (MQ=255)
                                                                     ttCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGTTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACggg                                                                         <  2:158914/149‑1 (MQ=255)
                                                                          gATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACACTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCt                                                                       <  5:190671/146‑1 (MQ=255)
                                                                              aGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGa                                                                 <  6:156610/148‑1 (MQ=255)
                                                                              aGTAGCATGCGGCAGCTTGCGGGGCACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGAc                                                                <  8:126400/149‑1 (MQ=255)
                                                                               gTAGCATGCGGCAGCTTTCGGGGAACCCCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCCCCAGCAGTTCCTGACCCACGCTGCCGCCTGTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGAc                                                                <  1:241795/149‑1 (MQ=255)
                                                                                  gCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGc                                                            <  5:221557/149‑1 (MQ=255)
                                                                                   cATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCa                                                           <  7:177489/149‑1 (MQ=255)
                                                                                         gCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGa                                                     >  8:131132/1‑149 (MQ=255)
                                                                                           aGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGAc                                                    >  3:86557/1‑148 (MQ=255)
                                                                                              tttCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTg                                                                                           <  6:209618/106‑1 (MQ=255)
                                                                                              tttCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTg                                                                                           >  5:209618/1‑106 (MQ=255)
                                                                                                  ggggAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATg                                            >  4:6224/1‑149 (MQ=255)
                                                                                                  ggggAACAGCCCGCTGTCGAGACACACGCAACAGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTGCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAAt                                             <  6:218436/148‑1 (MQ=255)
                                                                                                   gggAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGc                                           <  6:244705/149‑1 (MQ=255)
                                                                                                    ggAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGcc                                          >  2:139680/1‑149 (MQ=255)
                                                                                                       aCAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCctgc                                       >  5:44435/1‑149 (MQ=255)
                                                                                                       aCAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCctgc                                       >  5:39225/1‑149 (MQ=255)
                                                                                                          gCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCtgctgc                                    <  4:345412/149‑1 (MQ=255)
                                                                                                             cGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAgg                                 <  2:119480/149‑1 (MQ=255)
                                                                                                             cGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAgg                                 >  6:57113/1‑149 (MQ=255)
                                                                                                                 gTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTctg                             >  7:90549/1‑149 (MQ=255)
                                                                                                                 gTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTctg                             >  8:225151/1‑149 (MQ=255)
                                                                                                                  ttGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTctgc                            <  7:4781/149‑1 (MQ=255)
                                                                                                                  ttGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTctgc                            >  4:351207/1‑149 (MQ=255)
                                                                                                                  ttGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAACCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTctgc                            <  3:350300/149‑1 (MQ=255)
                                                                                                                      gaAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGGCGCGCCAGCAGTTCCTGCGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTctgcctgc                        <  7:231762/149‑1 (MQ=255)
                                                                                                                           acGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCAt                   >  4:107485/1‑149 (MQ=255)
                                                                                                                             gCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCAGGCTGCCGCC‑TTAGTCGTGGGCCAGCAGGTGCACTGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCa                 <  8:170790/149‑1 (MQ=255)
                                                                                                                             gCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCa                 <  8:343939/149‑1 (MQ=255)
                                                                                                                                     cGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCAGGCTGt           >  5:219315/1‑147 (MQ=255)
                                                                                                                                     cGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCAGGCTGTaa         <  4:86557/149‑1 (MQ=255)
                                                                                                                                         gTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCctg                                        <  6:86263/114‑1 (MQ=255)
                                                                                                                                         gTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCctg                                        >  5:86263/1‑114 (MQ=255)
                                                                                                                                            ggCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCAGGCTGTAAGTGTGa   <  7:131132/148‑1 (MQ=255)
                                                                                                                                            ggCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCAGGCTGTAAGTGTGAt  >  5:330109/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACATTACGCACCAGGTCGTCACGCCCGAACGAGCGCCCGACCAGCCAGCCAAAGCGGCTGAGCAGCAGCTTCTGGATCAGTAGCATGCGGCAGCTTTCGGGGAACAGCCCGCTGTTGAGAAACACGCAACTGGCGATGTCGGCACGCTGCTCATGGTGTCGCGCCAGCAGTTCCTGAGCCACGCTGCCGCC‑GTAGTCGTGGGCCAGCAGGTGCACGGGCTGCGCGACCTGCAGATGGACAAGCAATGCCTGCTGCAGGTCTGCCTGCTCCATCAGGCTGTAAGTGTGAT  >  NC_002947/5249301‑5249589

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