Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,625,957 G→C 5.0% Y515* (TAC→TAG)  PP_4944 ← carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,625,9570GC 5.0% 261.2 / 2.6 79Y515* (TAC→TAG) PP_4944carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

TCCGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACA  >  NC_002947/5625817‑5626101
                                                                                                                                            |                                                                                                                                                
tCCGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTcgc                                                                                                                                           <  2:95506/148‑1 (MQ=255)
 ccGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCgcgc                                                                                                                                         <  6:237173/149‑1 (MQ=255)
    aGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTg                                                                                                                                      <  2:303823/149‑1 (MQ=255)
     gCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGa                                                                                                                                     <  8:321202/149‑1 (MQ=255)
     gCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGa                                                                                                                                     <  5:220913/149‑1 (MQ=255)
       cGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAg                                                                                                                                    >  8:61376/1‑148 (MQ=255)
           gaaCATGTTCAGGGCATCTGTCGGCGAGCAGATCATCGGTTCAGCACTACGTTTCAGAGAGGTGTTCATCGACACACCGTTGCCAGTCAGGTCTTACAGTGCCTTCATGATGTCGTATTAGCGCGGGTTGTAGTcgctc                                                                                                                                         <  6:234706/139‑3 (MQ=255)
           gaaCATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCgcgc                                                                                                                                         >  5:234706/1‑139 (MQ=255)
                tGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCa                                                                                                                          >  5:326863/1‑149 (MQ=255)
                  ttCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCAcg                                                                                                                        >  3:210232/1‑149 (MQ=255)
                    cAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCAcgc                                                                                                                       <  7:61376/148‑1 (MQ=255)
                       ggCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAg                                                                                                                   >  6:129101/1‑149 (MQ=255)
                          aTCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGAGTACAACCCGCGCtac                                                                                                                                      <  1:159815/127‑3 (MQ=255)
                          aTCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGAGTACAACCCGCGCtac                                                                                                                                      >  2:165346/1‑125 (MQ=255)
                          aTCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGAGTACAACCCGCGCtac                                                                                                                                      <  1:165346/127‑3 (MQ=255)
                          aTCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGAGTACAACCCGCGCtac                                                                                                                                      >  2:159815/1‑125 (MQ=255)
                               tcggCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTcc                                                                                                           >  2:170767/1‑149 (MQ=255)
                                    gAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTg                                                                                                      <  6:187917/149‑1 (MQ=255)
                                      gCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAATCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGca                                                                                                    >  1:176342/1‑149 (MQ=255)
                                            catcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACtt                                                                                              >  3:205711/1‑149 (MQ=255)
                                             atcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTc                                                                                             >  4:130642/1‑149 (MQ=255)
                                              tcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGaa                                                                                                                    >  1:63539/1‑125 (MQ=255)
                                              tcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGaa                                                                                                                    <  2:63539/125‑1 (MQ=255)
                                              tcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTcc                                                                                            >  5:281962/1‑149 (MQ=255)
                                              tcGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATCTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTcc                                                                                            >  4:140686/1‑149 (MQ=255)
                                                ggTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCgg                                                                                          <  1:170767/149‑1 (MQ=255)
                                                 gTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCgg                                                                                          >  2:312255/1‑148 (MQ=255)
                                                 gTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCAGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCgg                                                                                          >  7:138956/1‑148 (MQ=255)
                                                        acgacgGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCggg                                                                                  >  6:270321/1‑149 (MQ=255)
                                                        acgacgGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCggg                                                                                  >  1:221223/1‑149 (MQ=255)
                                                            cgGGTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGtttt                                                                              <  4:8790/149‑1 (MQ=255)
                                                             ggTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTc                                                                             <  3:140686/149‑1 (MQ=255)
                                                             ggTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTc                                                                             <  4:205711/149‑1 (MQ=255)
                                                              gTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTcc                                                                            <  1:312255/149‑1 (MQ=255)
                                                                 cAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCAc                                                                          <  6:281962/148‑1 (MQ=255)
                                                                    cGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCtt                                                                      >  5:32184/1‑149 (MQ=255)
                                                                        gtgtTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGc                                                                  >  5:66795/1‑149 (MQ=255)
                                                                          gtTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGa                                                                >  1:81840/1‑149 (MQ=255)
                                                                                 gACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaa                                                         <  5:270321/149‑1 (MQ=255)
                                                                                  acacacCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaag                                                        >  6:155718/1‑149 (MQ=255)
                                                                                  acacacCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaag                                                        >  4:227486/1‑149 (MQ=255)
                                                                                  acacacCGTTGCCAATCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaag                                                        <  2:176342/149‑1 (MQ=255)
                                                                                      acCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtg                                                     >  1:158907/1‑148 (MQ=255)
                                                                                      acCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtg                                                     >  3:147922/1‑148 (MQ=255)
                                                                                       ccGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtg                                                     >  3:278514/1‑147 (MQ=255)
                                                                                        cGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtgaag                                                  <  5:129101/149‑1 (MQ=255)
                                                                                        cGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtgaag                                                  <  8:138956/149‑1 (MQ=255)
                                                                                        cGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtgaag                                                  >  6:313965/1‑149 (MQ=255)
                                                                                         gTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCgaaggtgaagg                                                 >  6:58910/1‑149 (MQ=255)
                                                                                             ccAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCAt                                             <  4:210232/149‑1 (MQ=255)
                                                                                             ccAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCAt                                             <  3:227486/149‑1 (MQ=255)
                                                                                              cAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCag                                             >  5:317408/1‑147 (MQ=255)
                                                                                                 tCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAAc                                         <  3:130642/149‑1 (MQ=255)
                                                                                                     gTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTc                                                                                                         <  8:336905/81‑1 (MQ=255)
                                                                                                      tCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGgcgc                                    <  6:66795/149‑1 (MQ=255)
                                                                                                       cTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGcgcg                                   >  3:347168/1‑149 (MQ=255)
                                                                                                          ccAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGcgcg                                   <  6:32184/146‑1 (MQ=255)
                                                                                                          ccAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGa                                >  6:31252/1‑149 (MQ=255)
                                                                                                            aGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATg                              <  4:278514/149‑1 (MQ=255)
                                                                                                             gTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATgg                             >  3:132203/1‑149 (MQ=255)
                                                                                                                ccTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTc                           >  3:208391/1‑148 (MQ=255)
                                                                                                                ccTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTc                           >  1:252516/1‑148 (MQ=255)
                                                                                                                 cTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGa                         >  2:310685/1‑149 (MQ=255)
                                                                                                                 cTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGa                         >  8:264715/1‑149 (MQ=255)
                                                                                                                  ttCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGAc                        >  8:345089/1‑149 (MQ=255)
                                                                                                                   tcatcaTGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGAcc                       >  8:258519/1‑149 (MQ=255)
                                                                                                                      tcatGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTg                    <  2:81840/149‑1 (MQ=255)
                                                                                                                      tcatGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTg                    <  2:158907/149‑1 (MQ=255)
                                                                                                                         tGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGatc                 <  5:155718/149‑1 (MQ=255)
                                                                                                                                gtagCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGc           <  4:147922/148‑1 (MQ=255)
                                                                                                                                  agCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGAACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCggg        >  3:193883/1‑149 (MQ=255)
                                                                                                                                   gcgcGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCgggg       >  4:263524/1‑149 (MQ=255)
                                                                                                                                   gcgcGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCgggg       <  3:263524/149‑1 (MQ=255)
                                                                                                                                   gcgcGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCggg        <  7:258519/148‑1 (MQ=255)
                                                                                                                                   gcgcGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCggg        <  2:221223/148‑1 (MQ=255)
                                                                                                                                    cgcgGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGc      <  5:31252/149‑1 (MQ=255)
                                                                                                                                      cgGGTTGTACTCGCGTTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGa    >  3:112655/1‑149 (MQ=255)
                                                                                                                                        ggTTGTACTCGCGCTTGAGCACCTTGGCAGGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACa  <  3:50206/149‑1 (MQ=255)
                                                                                                                                        ggTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTTCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACa  <  6:326863/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
TCCGAGCCGAAGAACATGTTCAGGGCATCGGTCGGCGAGCAGATCATCGGTTCACCACGACGGTTCAGCGAGGTGTTCAGCGACACACCGTTGCCAGTCAGGTCTTCCAGTGCCTTCATCATGTCGTAGTAGCGCGGGTTGTACTCGCGCTTGAGCACCTGGGCACGCGAAGTGCCGTCCTCGTGCACCACTTCCGGCACGCGGGTTTTCCACTCTTCAGCGACTTCGAAGGTGAAGGTCATGAACGGCGCGGGATGGTCGACCTTGATCATCTGCGGGGCGACA  >  NC_002947/5625817‑5626101

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