Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,915,826 A→T 7.2% intergenic (+254/‑35) fklB‑II → / → PP_1715 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,915,8260AT 7.2% 183.8 / 3.4 65intergenic (+254/‑35)fklB‑II/PP_1715FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base T (3/2);  total (35/30)
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.48e-01

GAAGTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCGG  >  NC_002947/1915678‑1915973
                                                                                                                                                    |                                                                                                                                                   
gaagTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACaa                                                                                                                                                     <  8:319297/149‑1 (MQ=255)
  agTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAAt                                                                                                                                                    >  4:272802/1‑148 (MQ=255)
   gTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATa                                                                                                                                                   >  1:325219/1‑148 (MQ=255)
   gTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGGCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATaa                                                                                                                                                  <  6:268919/149‑1 (MQ=255)
    tAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATaaa                                                                                                                                                 >  5:36603/1‑149 (MQ=255)
    tAGGATGTTTCTGACGCAAATCGTTTGCGACGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAAAGCCCGATCTGACCGGGCCCCTTGCGTGGCAGGCACATGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATaaa                                                                                                                                                 >  1:209395/1‑149 (MQ=255)
      ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtg                                                                                                                                               >  1:290708/1‑149 (MQ=255)
       gATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                              <  7:117481/149‑1 (MQ=255)
       gATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAt                                                                                                                                                <  7:176730/147‑1 (MQ=255)
       gATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAt                                                                                                                                                >  8:176730/1‑147 (MQ=255)
                 aCGAAAATCGTTTGCGCCGCAGGTAGTCATTGCGTGAAACGAATGTGTAAAAAACGCCCGATCTGACCGGGCCCATTGCCTGGCAGGAACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGATATCCACAATAAATGTGGATAACAtt                                                                                                                                    >  4:110215/1‑149 (MQ=255)
                  cGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCttt                                                                                                                                   <  6:36603/149‑1 (MQ=255)
                   gCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCtttt                                                                                                                                  >  3:50825/1‑149 (MQ=255)
                   gCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCttt                                                                                                                                   >  6:332106/1‑148 (MQ=255)
                   gCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCttt                                                                                                                                   <  5:332106/148‑1 (MQ=255)
                    cAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTc                                                                                                                                 <  1:204994/149‑1 (MQ=255)
                     aaaTCGTTTGCGCCGCAGGTACTCATCGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTAGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATCTGGATAACGATTTcc                                                                                                                                >  3:68957/1‑149 (MQ=255)
                             tGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGa                                                                                                                         <  1:43236/148‑1 (MQ=255)
                               cgccgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGc                                                                                                                      <  2:290708/149‑1 (MQ=255)
                                  cgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTcc                                                                                                                   <  8:114669/149‑1 (MQ=255)
                                       gTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAg                                                                                                              >  3:201279/1‑149 (MQ=255)
                                                    gAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAt                                                                                                 >  8:111460/1‑149 (MQ=255)
                                                     aaaCGACTGTGTAAAAAACGCCCGATCTGCCCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtt                                                                                                <  2:280336/149‑1 (MQ=255)
                                                       aCGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCACTGTTCACAAGGTTATCCACAATAAATGTgg                                                                                                                                            >  7:265838/1‑103 (MQ=255)
                                                       aCGACTGTCTAAAAAACGCCCGATCTGAGCTGGCCCCTTGCCTGGCAGGCACCTGGGGGCGTAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTgg                                                                                                                                            <  8:265838/103‑1 (MQ=255)
                                                           cTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGccc                                                                                          <  2:156413/149‑1 (MQ=255)
                                                              tgtAAAAAACGCCCGATCTGCCCGGGCCCCTTGCCTTGCAGTCACCTGGGGTCGGAAACGTAGCACTTTTCACAACGGTATCCACAATAATTGTTGATAAACTTTTCCTCTCTGAGGCTCCATGAGCGCACCCTGGAATTTTTCTCGtt                                                                                       <  2:209395/149‑1 (MQ=255)
                                                              tgtAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtttt                                                                                              >  6:61243/1‑142 (MQ=255)
                                                              tgtAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtttt                                                                                              <  5:61243/142‑1 (MQ=255)
                                                               gtAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGttt                                                                                      >  7:201967/1‑149 (MQ=255)
                                                                tAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTc                                                                                     >  6:161117/1‑149 (MQ=255)
                                                                  aaataCGCCCGTTCTGACCGGGTCCCTTGGCTGGCAGGGTCCTGGGGGCGGTTACGTGGCCCTGTTGACAAGGGTATCCACATTAAATGTGGCTAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCt                                                                                   <  4:68957/149‑1 (MQ=255)
                                                                    aaaCGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGc                                                                                 >  4:44006/1‑149 (MQ=255)
                                                                                 gACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                             <  1:77948/108‑1 (MQ=255)
                                                                                 gACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                             >  2:77948/1‑108 (MQ=255)
                                                                                 gACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAgcg                                                                    >  3:16705/1‑149 (MQ=255)
                                                                                  aCCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAgcgc                                                                   >  4:26571/1‑149 (MQ=255)
                                                                                      ggCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCGGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGc                                                               >  7:295988/1‑149 (MQ=255)
                                                                                       gCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTg                                                                <  7:87054/147‑1 (MQ=255)
                                                                                       gCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTg                                                                >  8:87054/1‑147 (MQ=255)
                                                                                              gCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCGTGCCATTGGCCGAGCGCCTGCTCAGcc                                                         <  6:359222/147‑1 (MQ=255)
                                                                                                cTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGAGCCTGCTCAGCCGCgg                                                     >  8:139961/1‑149 (MQ=255)
                                                                                                 tGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCt                                                                                   <  7:213130/118‑1 (MQ=255)
                                                                                                 tGGCAGGCACCTCGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCACGGAGGCTCCATGAGCGCACCATGGAATTTTGCCCGTTTACt                                                                                   >  8:213130/1‑118 (MQ=255)
                                                                                                         aCCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGcc                                            <  4:50825/149‑1 (MQ=255)
                                                                                                           cTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCgg                                          >  4:4073/1‑149 (MQ=255)
                                                                                                              ggggCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTc                                                                                                                                 <  3:337094/59‑1 (MQ=255)
                                                                                                              ggggCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTc                                                                                                                                 >  4:337094/1‑59 (MQ=255)
                                                                                                                 gCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAACTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgt                                    <  8:42347/149‑1 (MQ=255)
                                                                                                                  cGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgtt                                   <  4:201279/149‑1 (MQ=255)
                                                                                                                   ggAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgttg                                  >  7:113886/1‑149 (MQ=255)
                                                                                                                   ggAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgttg                                  >  6:104347/1‑149 (MQ=255)
                                                                                                                       aCGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTcg                              <  8:201967/149‑1 (MQ=255)
                                                                                                                         gTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTcgcg                            <  8:295988/149‑1 (MQ=255)
                                                                                                                            gCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTg                         >  7:75781/1‑149 (MQ=255)
                                                                                                                                  ttCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGgcgcg                    <  3:26571/148‑1 (MQ=255)
                                                                                                                                       aaGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGcc                                            <  5:369883/119‑1 (MQ=255)
                                                                                                                                       aaGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGcc                                            >  6:369883/1‑119 (MQ=255)
                                                                                                                                         ggTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGcc            >  1:334487/1‑149 (MQ=255)
                                                                                                                                           ttATCCACATTTATTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgtt                                   >  4:215264/1‑124 (MQ=255)
                                                                                                                                           ttATCCACATTTATTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgtt                                   <  3:215264/124‑1 (MQ=255)
                                                                                                                                            tATCCACATTTATTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCcgc         >  7:253899/1‑149 (MQ=255)
                                                                                                                                             aTCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGGGGCGCGCAAAGGCCcgc         <  7:254484/148‑1 (MQ=255)
                                                                                                                                              tCCACATTTATTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCgcgc       >  4:107076/1‑149 (MQ=255)
                                                                                                                                                 acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTc    <  3:44006/149‑1 (MQ=255)
                                                                                                                                                 acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTc    >  6:312781/1‑149 (MQ=255)
                                                                                                                                                   aaTAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCAACCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCgg  <  3:110215/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGTAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCGG  >  NC_002947/1915678‑1915973

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