Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,704,636 T→A 11.5% M1M (TTG→ATG) †‡ PP_1499 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,704,6360TA11.5% 149.8 / 5.5 63M1M (TTG→ATG) †‡PP_1499hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base A (5/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CCGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAA  >  NC_002947/1704488‑1704784
                                                                                                                                                     |                                                                                                                                                    
ccGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCt                                                                                                                                                      >  8:100572/1‑149 (MQ=255)
     ggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGccc                                                                                                                                                 >  1:268436/1‑149 (MQ=255)
          acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGCCGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                            <  7:283/149‑1 (MQ=255)
          acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                            <  3:185171/149‑1 (MQ=255)
           cggcgacggcgTGCTGATCGACAGCCCGGCGGAGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGGTCCAGCGCTCGGGCAACTGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCTAag                                                                                                                                           <  5:158089/149‑1 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                            >  8:289078/1‑147 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                          <  2:28080/149‑1 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                          <  3:18769/149‑1 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                          <  5:85572/149‑1 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCCCGGGCTTGCCCCCGaaca                                                                                                                                          >  2:225850/1‑147 (MQ=255)
                        cTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGAGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAgg                                                                                                                              >  5:24828/1‑149 (MQ=255)
                         tGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGt                                                                                                                             >  8:240595/1‑149 (MQ=255)
                         tGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                             <  5:132830/133‑1 (MQ=255)
                         tGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                             >  6:132830/1‑133 (MQ=255)
                           aTCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGt                                                                                                                             <  1:200386/147‑1 (MQ=255)
                              gACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCg                                                                                                                         <  7:100572/148‑1 (MQ=255)
                                cAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGGTTCgtgt                                                                                                                      >  7:121985/1‑149 (MQ=255)
                                 aGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGAGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTa                                                                                                                     <  6:24828/149‑1 (MQ=255)
                                   cccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTGGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACt                                                                                                                   <  1:225850/149‑1 (MQ=255)
                                            gcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAg                                                                                                          >  8:47633/1‑149 (MQ=255)
                                            gcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAg                                                                                                          >  8:87483/1‑149 (MQ=255)
                                                    tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGt                                                                                                  >  1:36207/1‑149 (MQ=255)
                                                      gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTTCCGGTgg                                                                                                <  2:190750/149‑1 (MQ=255)
                                                      gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTgg                                                                                                >  8:187987/1‑149 (MQ=255)
                                                            cTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAAt                                                                                          <  7:247478/149‑1 (MQ=255)
                                                               gacgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATAcc                                                                                       >  3:310863/1‑149 (MQ=255)
                                                                acgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                      <  3:270626/149‑1 (MQ=255)
                                                                    cccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGCGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCa                                                                                   >  3:235113/1‑148 (MQ=255)
                                                                     ccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGGGCAGAGGCCGGGTCAGGTTTGGTGGACTAGAAAGGAGTTGCCGGTGGGCGAATAACAGCcaa                                                                                 >  4:313768/1‑149 (MQ=255)
                                                                       ggCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAc                                                                                >  7:70097/1‑148 (MQ=255)
                                                                          cAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                              >  8:192693/1‑131 (MQ=255)
                                                                          cAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                              <  7:192693/131‑1 (MQ=255)
                                                                           aaCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATcc                                                                           <  2:36207/149‑1 (MQ=255)
                                                                                   ccTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTg                                                                   >  1:67274/1‑149 (MQ=255)
                                                                                    cTGTTCGTCGTTGTGGCGCGGTTCCAGCGCTCGGGAAAATGAGCGTGTTGGCGTCTTGGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTg                                                                   <  8:70097/148‑1 (MQ=255)
                                                                                      gTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTcc                                                                >  3:216775/1‑149 (MQ=255)
                                                                                      gTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTCATGGCGTCTTAACGCGCATTACCCACAACAATGCCCTACAAGTTTCGAGTACTAGTAAGGAATTAGAGCTGGACTCATACCAGCAAACCAACAgcgcacagtac                                                                >  4:117268/1‑142 (MQ=255)
                                                                                           tcgAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCg                                                           >  6:147117/1‑149 (MQ=255)
                                                                                              aaGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAt                                                        <  7:187987/149‑1 (MQ=255)
                                                                                               aGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtc                                                       <  7:289078/149‑1 (MQ=255)
                                                                                                    cgcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTcgc                                                  <  4:310863/149‑1 (MQ=255)
                                                                                                     gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAACGTcgcg                                                 >  2:281057/1‑149 (MQ=255)
                                                                                                      cgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                                <  2:268436/149‑1 (MQ=255)
                                                                                                         gtCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCGTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGcac                                             <  7:26753/148‑1 (MQ=255)
                                                                                                           ccAGCGCTCCGGCAAATGACCGTGTTGGCGTCTTCCCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTACCCGTGGACGAATACCAGCAAACCATCCGCGCCCAGTCCGACCGCAACGTCGCGGca                                              >  8:82870/1‑146 (MQ=255)
                                                                                                            cAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGa                                          >  7:242007/1‑149 (MQ=255)
                                                                                                               cgcTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGc                                       <  2:117739/149‑1 (MQ=255)
                                                                                                                 cTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCATGCCCGCCAAGAGGCCCGTACAGGTTTCGAGTACTAGAAAGGAGTTCCCGGTGGACCAATACCAGCAAACCATACCCGCCCAGTCCGACCCCAACCGCGCGGCACACACGCCg                                     >  4:240089/1‑149 (MQ=255)
                                                                                                                      gCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATcc                                 <  1:122436/148‑1 (MQ=255)
                                                                                                                      gCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGACAGGTTTAGTGTACTTGAAAGGAGGTCCCGGTGGACGAATACCACCAAACCAAACGCGCCCAGTACGACCGCATCGTAGCGGCACAAACGCCGATcc                                 >  2:122436/1‑148 (MQ=255)
                                                                                                                       cAAATGTCCTTGTTGGCTTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTAGAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAGCATCCGCGCCCTGTCGGACCGAATCGTCGCGGCACAGACTCCGATCCgg                               <  3:117268/149‑1 (MQ=255)
                                                                                                                       cAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCgg                               <  8:121985/149‑1 (MQ=255)
                                                                                                                        aaaTGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATc                                  <  2:241425/145‑1 (MQ=255)
                                                                                                                        aaaTGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATc                                  >  1:241425/1‑145 (MQ=255)
                                                                                                                                   ttGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGAcgcg                   <  5:147117/149‑1 (MQ=255)
                                                                                                                                   ttGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGAcgc                    <  7:240595/148‑1 (MQ=255)
                                                                                                                                   ttGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTACGACCGCAACGTCGCGGCACAGACGCCGATCCGGGTCCATGAcgcg                   >  6:293553/1‑149 (MQ=255)
                                                                                                                                     ggCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAACGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCgg                  <  1:281057/148‑1 (MQ=255)
                                                                                                                                      gcctCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt                 >  3:98973/4‑148 (MQ=255)
                                                                                                                                          cTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt            >  7:144271/1‑149 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg           >  5:13564/1‑149 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac       <  7:47633/149‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacg      <  2:67274/149‑1 (MQ=255)
                                                                                                                                                  gCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacgac    <  8:242007/149‑1 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGacaa  >  8:245725/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAA  >  NC_002947/1704488‑1704784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: