Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,380,544 C→T 9.1% V72V (GTC→GTT)  tnpT‑I → cointegrate resolution protein T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,380,5440CT 9.1% 217.4 / 8.6 77V72V (GTC→GTT) tnpT‑Icointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base T (3/4);  total (38/39)
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.90e-01

ACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGG  >  NC_002947/3380399‑3380690
                                                                                                                                                 |                                                                                                                                                   
aCCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCgg                                                                                                                                                   <  6:46067/148‑1 (MQ=255)
aCCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGc                                                                                                                                                  >  1:126905/1‑149 (MQ=255)
   cAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACAGCCTCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCaa                                                                                                                                               >  6:218939/1‑149 (MQ=255)
      gCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGCCAACTg                                                                                                                                            <  3:195439/149‑1 (MQ=255)
       cATCGACGCGGTACGGGTCTAACTTGGCAATACCTGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTgg                                                                                                                                           <  5:439761/149‑1 (MQ=255)
       cATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTgg                                                                                                                                           <  3:280440/149‑1 (MQ=255)
        aTCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGc                                                                                                                                          <  3:348548/149‑1 (MQ=255)
           gACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCGACCAGGTTGGCCAACTgttcg                                                                                                                                        <  8:267658/148‑5 (MQ=255)
                     gggTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTagag                                                                                                                             >  5:298991/1‑149 (MQ=255)
                        tCGAACTTGGCAATCCCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACTGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAgaggag                                                                                                                          <  8:167898/149‑1 (MQ=255)
                           aaCTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAgaggag                                                                                                                          >  1:374051/1‑146 (MQ=255)
                            aCTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCGACTGGCAGATCAGCTAGAGGAGGATg                                                                                                                      >  8:27862/1‑149 (MQ=255)
                                    aTACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGcc                                                                                                              <  2:126905/149‑1 (MQ=255)
                                        cGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCg                                                                                                          >  1:403736/1‑149 (MQ=255)
                                         ggCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACAGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGt                                                                                                         >  1:131995/1‑149 (MQ=255)
                                           cTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTgg                                                                                                       >  2:306079/1‑149 (MQ=255)
                                              gAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGc                                                                                                      <  1:112370/147‑1 (MQ=255)
                                               aaaaCCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGcc                                                                                                     >  1:71781/1‑147 (MQ=255)
                                                  accaccATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGa                                                                                                <  6:298991/149‑1 (MQ=255)
                                                  accaccATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGAGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAg                                                                                                  >  6:192523/1‑147 (MQ=255)
                                                    caccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgc                                                                                              >  3:31623/1‑149 (MQ=255)
                                                    caccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgc                                                                                              >  4:212452/1‑149 (MQ=255)
                                                    caccaTTCACCGCTACCTGAAAGAACTCGAAGACGCAGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgc                                                                                              >  1:70498/1‑149 (MQ=255)
                                                      ccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCGACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgcgc                                                                                            <  7:27862/149‑1 (MQ=255)
                                                      ccaTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCCACTGGCAGATCCGCTAGAGGAGGCTGCGCAGGCCACCGTGGCCCAGGAgcgc                                                                                            >  6:428726/1‑149 (MQ=255)
                                                       caTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGcgcg                                                                                           <  1:2330/149‑1 (MQ=255)
                                                       caTTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTAGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGcgcg                                                                                           <  1:216335/149‑1 (MQ=255)
                                                           cACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaac                                                                                        <  3:69800/148‑1 (MQ=255)
                                                            aCCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaacaa                                                                                      <  2:374051/149‑1 (MQ=255)
                                                             ccGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGaacaa                                                                                      >  2:294608/1‑148 (MQ=255)
                                                              cGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACt                                                                                    <  4:31623/149‑1 (MQ=255)
                                                                 tACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑Ac                                                                                 >  3:226776/1‑148 (MQ=255)
                                                                   ccTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAAc                                                                              >  3:141703/1‑149 (MQ=255)
                                                                      gAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAAcgc                                                                            >  7:28224/1‑148 (MQ=255)
                                                                      gAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAAcgc                                                                            >  7:387549/1‑148 (MQ=255)
                                                                       aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAAcgcg                                                                           >  4:226082/1‑148 (MQ=255)
                                                                       aaaGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAAcgcg                                                                           <  2:71781/148‑1 (MQ=255)
                                                                           aaCTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACgg                                                                      <  2:403736/149‑1 (MQ=255)
                                                                             cTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgcgc                                                                                            >  6:192050/1‑126 (MQ=255)
                                                                             cTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAgcgc                                                                                            <  5:192050/126‑1 (MQ=255)
                                                                              tCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTc                                                                   >  2:47448/1‑149 (MQ=255)
                                                                               cGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTc                                                                   >  8:83190/1‑148 (MQ=255)
                                                                                gAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGa                                                                 >  3:16939/1‑149 (MQ=255)
                                                                                    aCGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTAc                                                             <  4:141703/149‑1 (MQ=255)
                                                                                      gCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCa                                                           <  7:231467/149‑1 (MQ=255)
                                                                                       ccGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTAc                                                             >  2:83231/1‑146 (MQ=255)
                                                                                             gcgGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCaa                                                    <  1:294608/149‑1 (MQ=255)
                                                                                                  cGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCg                                               >  5:258674/1‑149 (MQ=255)
                                                                                                    aaaTACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGc                                             >  6:58843/1‑149 (MQ=255)
                                                                                                         ccgccgCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                         <  2:70498/148‑1 (MQ=255)
                                                                                                         ccgccgCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                         <  3:226082/148‑1 (MQ=255)
                                                                                                          cgccgcGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTAACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGcc                                        <  1:83231/149‑1 (MQ=255)
                                                                                                             cgcgATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGa                                      >  4:251405/1‑147 (MQ=255)
                                                                                                                     cGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                         <  5:102453/136‑1 (MQ=255)
                                                                                                                     cGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGc                                         >  6:102453/1‑136 (MQ=255)
                                                                                                                      gCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATc                            >  7:24237/1‑148 (MQ=255)
                                                                                                                      gCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATc                            >  6:225147/1‑148 (MQ=255)
                                                                                                                      gCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATc                            >  4:429764/1‑148 (MQ=255)
                                                                                                                          aGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCATATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGTCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCcagca                       <  5:192523/149‑1 (MQ=255)
                                                                                                                          aGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCcag                         <  6:28402/147‑1 (MQ=255)
                                                                                                                             gAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGtt                    <  7:83190/149‑1 (MQ=255)
                                                                                                                             gAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTAGCCGAAAGCCGAATCCAGCAGtt                    <  4:413050/149‑1 (MQ=255)
                                                                                                                               aCAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTgg                  <  1:437921/149‑1 (MQ=255)
                                                                                                                                cAGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAg                                                   >  4:81624/1‑115 (MQ=255)
                                                                                                                                cAGTTGGCCGACCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAg                                                   <  3:81624/115‑1 (MQ=255)
                                                                                                                                cAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGa                 <  8:387549/149‑1 (MQ=255)
                                                                                                                                  gTTGGCCGACCAGGTTGGCCAACTGTCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGa                 >  3:148237/1‑147 (MQ=255)
                                                                                                                                  gTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGTGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGa                 >  1:64805/1‑147 (MQ=255)
                                                                                                                                  gTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGaaa               <  5:218939/149‑1 (MQ=255)
                                                                                                                                     ggCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGt             <  4:226776/148‑1 (MQ=255)
                                                                                                                                      gCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCa           <  1:47448/149‑1 (MQ=255)
                                                                                                                                       ccgaCCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGaaa                                    >  3:4561/4‑123 (MQ=255)
                                                                                                                                       ccgaCCAGGTTGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGaaa                                    <  4:4561/120‑1 (MQ=255)
                                                                                                                                         aaCCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCagag        >  7:83419/1‑149 (MQ=255)
                                                                                                                                         aaCCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCagag        <  1:306079/149‑1 (MQ=255)
                                                                                                                                         aaCCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCagag        <  8:28224/149‑1 (MQ=255)
                                                                                                                                               gTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAACGCCGAATCCAGCAGTTGGAATGGCAGAGCAGCgg  >  4:40526/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACCCAAGCATCGACGCGGTACGGGTCGAACTTGGCAATACCGGCTCGAAAACCACCATTCACCGCTACCTGAAAGAACTCGAAGACGCCGAGCGCGGCCGAAATACCGCCGCGATACCGCTCAGCGAACAGTTGGCCAACCTGGTCGGCCAACTGGCAGATCAGCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACT‑ACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGG  >  NC_002947/3380399‑3380690

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