Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,704,636 T→A 11.8% 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.8% 169.3 / 13.2 68M1M (TTG→ATG) †‡PP_1499hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (5/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

CCGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGAC  >  NC_002947/1704488‑1704779
                                                                                                                                                    |                                                                                                                                                
ccGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCt                                                                                                                                                  >  3:128427/1‑149 (MQ=255)
   gaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                               <  2:104872/149‑1 (MQ=255)
   gaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTAGCGGGCTTgg                                                                                                                                               >  7:304600/1‑148 (MQ=255)
      gaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCcg                                                                                                                                            >  7:178027/1‑149 (MQ=255)
         aaCGTCGACGGCGTGCTGATCGACAGACAGGCGGCTCGGGTGCTGGCCTACCTGGACGACCCGGCCAAATCAGACCTGTTCGTCGAAGTGGCGCGTTTCCAGCGCTCGGGAAAATGAACGTGTTGGCGTCTTCTCGGGCTTTCAcgcg                                                                                                                                          >  2:231952/1‑148 (MQ=255)
         aaCGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCcgcg                                                                                                                                          >  6:88776/1‑148 (MQ=255)
          acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                        <  4:86016/149‑1 (MQ=255)
          acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                        <  4:243819/149‑1 (MQ=255)
          acggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                        <  7:50521/149‑1 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                        >  2:58178/1‑147 (MQ=255)
            ggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                      <  4:128427/149‑1 (MQ=255)
               gacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                   >  7:256412/1‑149 (MQ=255)
                     gTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTAc                                                                                                                             <  5:170296/149‑1 (MQ=255)
                          gATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                               >  2:118983/1‑126 (MQ=255)
                          gATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                               <  1:118983/126‑1 (MQ=255)
                               aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgtg                                                                                                                   <  8:178027/149‑1 (MQ=255)
                               aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgtg                                                                                                                   >  1:24286/1‑149 (MQ=255)
                               aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTAGCGGGCTTGACCGAGAAGAGGCCGGTACACGTATCgtg                                                                                                                   >  6:277579/1‑149 (MQ=255)
                                  gCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTGGCGGGCGTGGCCGCGCAGGGGCCGGGTCAGGGTTGgtggt                                                                                                                 >  3:192052/1‑146 (MQ=255)
                                  ccccGGCGGCGGGGGTGCTGGCCTACCTGGACCACCCGGCAAACTCCGACCTGTTCGGCGATGGGGCGCGGTTCCAGCGCTCGGGAAAAGGCCCGTGTTGGGGTCTTCGGGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTa                                                                                                                 <  4:192052/147‑1 (MQ=255)
                                    ccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTa                                                                                                              <  2:292084/149‑1 (MQ=255)
                                        cggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGaaa                                                                                                          <  2:27865/149‑1 (MQ=255)
                                         ggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAg                                                                                                         <  7:94804/149‑1 (MQ=255)
                                           cgcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                       <  2:197393/149‑1 (MQ=255)
                                           cgcgGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                       <  5:175845/149‑1 (MQ=255)
                                                  gCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                >  3:314312/1‑149 (MQ=255)
                                                     ggCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                             >  3:149304/1‑149 (MQ=255)
                                                      gCCTACCTGGACGACCCGTCCAACTCCGACGTGTTCGTCGAAGTGGCGCGGTTTCAGCGCTCGGGCAAATTACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCTAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                             <  8:154798/148‑1 (MQ=255)
                                                      gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTGGGGAAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGa                                                                                                       <  5:245270/138‑1 (MQ=255)
                                                      gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGGTCCCGGTgg                                                                                            >  8:104593/1‑149 (MQ=255)
                                                        cTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                          >  4:215984/1‑149 (MQ=255)
                                                        cTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                          <  4:7772/149‑1 (MQ=255)
                                                         tACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGAc                                                                                          >  7:90499/1‑148 (MQ=255)
                                                           ccTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGAGGAGCGGCCCGTTCCCGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGGCgc                                                                                        >  3:346118/1‑147 (MQ=255)
                                                             tGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGACAg                                                                                                                           >  5:220481/1‑111 (MQ=255)
                                                                acgacCGGGCCAACTCGGAGAGGTTCGTCTAAGTGGCGCGGTTCCAGCGCTCGGTCAAATTACGGTGTTGGGTTCTTCGCGGGGTTGCCCGCTATGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                  <  2:159190/149‑1 (MQ=255)
                                                                acgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                  <  3:141849/149‑1 (MQ=255)
                                                                  gacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGc                                                                                <  3:215984/149‑1 (MQ=255)
                                                                        gCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑Ca                                                                         <  8:256412/149‑1 (MQ=255)
                                                                           aaCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATcc                                                                      <  1:58178/149‑1 (MQ=255)
                                                                                cGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCaaa                                                                             >  7:345433/1‑138 (MQ=255)
                                                                                cGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGACCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCaaa                                                                             <  8:345433/138‑1 (MQ=255)
                                                                                  aCCTGTTGGTCGATGGGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCt                                                               <  8:304600/149‑1 (MQ=255)
                                                                                  aCCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCt                                                               <  5:88776/149‑1 (MQ=255)
                                                                                          gtcgAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATACGCGCCCTGTCCGAcc                                                       >  1:8121/1‑149 (MQ=255)
                                                                                                     gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTcgcg                                            >  6:364985/1‑149 (MQ=255)
                                                                                                     gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTcgcg                                            <  8:90499/149‑1 (MQ=255)
                                                                                                       ggTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACGCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                           <  1:231952/149‑1 (MQ=255)
                                                                                                             aGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGATGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGAc                                    <  4:149304/149‑1 (MQ=255)
                                                                                                             aGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGAc                                    >  1:307030/1‑149 (MQ=255)
                                                                                                              gcgcTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                           <  4:41359/141‑1 (MQ=255)
                                                                                                              gcgcTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                           >  3:41359/1‑141 (MQ=255)
                                                                                                                 cTCGGGCAAATGACCGTGTTGGCGTCTTCCCCGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCCCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGcc                                 >  8:31367/1‑148 (MQ=255)
                                                                                                                                    cggCCTCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCgg             >  2:309182/2‑149 (MQ=255)
                                                                                                                                      gcctCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt            >  4:150229/4‑148 (MQ=255)
                                                                                                                                      gCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt            <  4:314312/148‑1 (MQ=255)
                                                                                                                                      gCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGt            <  2:24286/148‑1 (MQ=255)
                                                                                                                                        gTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGaa         <  5:364985/149‑1 (MQ=255)
                                                                                                                                         tCTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAg        >  4:151509/1‑149 (MQ=255)
                                                                                                                                         tCTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAg        >  3:215392/1‑149 (MQ=255)
                                                                                                                                          cTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt       >  6:326453/1‑149 (MQ=255)
                                                                                                                                          cTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGGGGACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CCTCCGCGCCCTGTCCGACCGCATCGGCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt       <  7:255985/149‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg      >  1:144352/1‑149 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg      >  2:234079/1‑149 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCtt                    >  3:290615/1‑135 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCtt                    <  4:290615/135‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGga   >  8:11628/1‑149 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac  <  4:346118/149‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgac  >  2:108397/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAAC‑CATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGAC  >  NC_002947/1704488‑1704779

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