Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,851,692 T→A 5.4% intergenic (+84/+35) apt → / ← recR adenine phosphoribosyltransferase/recombination protein RecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,851,6920TA 5.4% 247.4 / 2.8 87intergenic (+84/+35)apt/recRadenine phosphoribosyltransferase/recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base A (3/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

GGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGCC  >  NC_002947/4851547‑4851840
                                                                                                                                                  |                                                                                                                                                      
ggTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTAcc                                                                                                                                                     >  1:284798/1‑149 (MQ=255)
ggTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTAcc                                                                                                                                                     >  8:318835/1‑149 (MQ=255)
ggTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCCGCTTGGCGGCCCATTCGCGGGTAAACCCGCGCCTACAAAACATACGCGGCTTTCGTAGGAGCGGGTTTAcc                                                                                                                                                     <  7:11784/149‑1 (MQ=17)
     cAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                               <  2:407758/149‑2 (MQ=255)
     cAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCcgcg                                                                                                                                                >  5:29712/1‑148 (MQ=255)
      aaGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTAccaccg                                                                                                                                                 >  1:155300/1‑143 (MQ=255)
      aaGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                              >  4:287470/1‑149 (MQ=255)
      aaGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGATAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                              >  2:361571/1‑149 (MQ=255)
        gCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTATTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAAt                                                                                                                                             >  1:341899/1‑149 (MQ=255)
         cGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGGATGTTTTGTAGGCGCGGGTTTACCCGCGAATgg                                                                                                                                           >  7:7348/1‑149 (MQ=255)
             cGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGgccg                                                                                                                                       <  6:30701/149‑1 (MQ=255)
             cGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCGTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTCCGAAACCCGCGTATGTTTTGTAGTCTAGGGTTTACCCGCGAATGGgccg                                                                                                                                       <  7:66375/149‑1 (MQ=39)
                     cTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAGAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCgg                                                                                                                               >  4:241771/1‑149 (MQ=255)
                      tACCTTCTGCCTGACCGAGTTTTCCCTTAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGc                                                                                                                              >  4:417680/1‑149 (MQ=255)
                      tACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGc                                                                                                                              >  8:308767/1‑149 (MQ=255)
                       aCCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGcc                                                                                                                             >  2:125278/1‑149 (MQ=255)
                        ccTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGc                                                                                                                              <  8:234049/147‑1 (MQ=255)
                        ccTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGGGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTTCGCAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGGCGCTTGGCGGcc                                                                                                                             <  7:29476/148‑1 (MQ=38)
                              gCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATg                                                                                                                       >  5:281180/1‑148 (MQ=255)
                               ccTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAAt                                                                                                                        <  6:104245/146‑1 (MQ=255)
                               ccTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                      <  6:233136/148‑1 (MQ=255)
                               ccTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGCTTGGCGGcccattc                                                                                                                        <  2:155300/146‑5 (MQ=38)
                                 tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGCTTGGCGGccc                                                                                                                            <  2:29182/140‑1 (MQ=255)
                                 tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGCTTGGCGGccc                                                                                                                            >  1:29182/1‑140 (MQ=255)
                                  gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                      >  5:319098/1‑145 (MQ=255)
                                  gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                      <  6:319098/145‑1 (MQ=255)
                                     cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCagag                                                                                                               >  3:128611/1‑149 (MQ=255)
                                     cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCaga                                                                                                                <  4:47620/148‑1 (MQ=255)
                                     cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCa                                                                                                                  >  3:47620/1‑146 (MQ=255)
                                          tttCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTATGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGaaa                                                                                                          >  7:164688/1‑149 (MQ=255)
                                           ttCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAAt                                                                                                         >  8:27079/1‑149 (MQ=255)
                                            tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATg                                                                                                        >  2:36908/1‑149 (MQ=255)
                                              ccTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATgg                                                                                                       <  4:128611/148‑1 (MQ=255)
                                              ccTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATgg                                                                                                       <  7:199392/148‑1 (MQ=255)
                                                      aGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGcc                                                                                              <  6:29712/149‑1 (MQ=255)
                                                       gTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGccc                                                                                             <  2:248408/149‑1 (MQ=255)
                                                         aTTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTg                                                                                           >  8:399434/1‑149 (MQ=255)
                                                          ttGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCTGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTg                                                                                           <  1:114561/149‑1 (MQ=255)
                                                          ttGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGc                                                                                          <  2:310670/149‑1 (MQ=255)
                                                             aGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCaaa                                                                                       <  3:417680/149‑1 (MQ=255)
                                                             aGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCGACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAAGCAGAGCGAAATGGGCCGGCGCCCTGccac                                                                                       >  8:62951/1‑146 (MQ=255)
                                                             aGTGGGCTTTTGGGGCGGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCaa                                                                                        <  4:410255/148‑1 (MQ=255)
                                                               tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                                                                                      <  2:341899/148‑1 (MQ=255)
                                                                          gggCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAg                                                                           <  7:318835/149‑1 (MQ=255)
                                                                           ggCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAgg                                                                         <  7:308767/149‑1 (MQ=255)
                                                                              tgctTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTa                                                                      >  7:388089/1‑149 (MQ=255)
                                                                               gctTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCag                                                                                                                 <  2:97879/106‑1 (MQ=255)
                                                                               gctTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCag                                                                                                                 >  1:97879/1‑106 (MQ=255)
                                                                                 ttttCAGCCGATTCGCTGATAAACCCTCTCCTACTAAAGCCGTGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccg                                                                   <  1:361571/149‑1 (MQ=38)
                                                                                   tGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgcc                                                                 >  4:71449/1‑149 (MQ=255)
                                                                                      aGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTc                                                              >  5:402891/1‑149 (MQ=255)
                                                                                       gCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc                                                             >  2:79992/1‑149 (MQ=255)
                                                                                        cccATTCGCGGGTACACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                            <  7:62951/149‑1 (MQ=255)
                                                                                        cccATTCGCGGGTAAACCCGCTCCTACGAGAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                            <  3:241771/149‑1 (MQ=255)
                                                                                        cccATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGcc                                                                                              >  6:161537/1‑116 (MQ=255)
                                                                                        cccATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGcc                                                                                              <  5:161537/116‑1 (MQ=255)
                                                                                        cccATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                            >  8:176401/1‑149 (MQ=255)
                                                                                          cATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccacc                                                          <  2:45470/149‑1 (MQ=255)
                                                                                          cATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccac                                                           >  4:144271/1‑148 (MQ=255)
                                                                                           aTTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccacca                                                         <  7:176401/149‑1 (MQ=255)
                                                                                             tCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtc                                                       >  7:229087/1‑149 (MQ=255)
                                                                                               gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtctc                                                     >  2:405989/1‑149 (MQ=255)
                                                                                               gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGCCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtctc                                                     <  1:36908/149‑1 (MQ=255)
                                                                                               gagggTAAACCCGCTCCTACGAATGCCGCGTATGTTTTGTAGGCGCGTGTTTACCCTCGAATGTGCCGCCAAGCGGCCCCAATGAATCAGAGAGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtctc                                                     <  8:7348/147‑1 (MQ=255)
                                                                                                cgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtctc                                                     >  1:421682/1‑148 (MQ=255)
                                                                                                cgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTcc                                                    <  7:27079/149‑1 (MQ=255)
                                                                                                 gggTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCa                                                   <  3:287470/149‑1 (MQ=255)
                                                                                                  ggTAAACCCGCTCCTACGAACGCCGCGTTTGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAg                                                  <  6:402891/149‑1 (MQ=255)
                                                                                                       aCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgc                                             >  1:205820/1‑149 (MQ=255)
                                                                                                            cTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGccc                                        <  2:284798/149‑1 (MQ=255)
                                                                                                             tCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCgccgcc                                         >  2:322527/1‑147 (MQ=255)
                                                                                                               cTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt                                     <  8:388089/149‑1 (MQ=255)
                                                                                                                  cGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                   <  2:421682/148‑1 (MQ=255)
                                                                                                                  cGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTg                                    <  2:375614/147‑1 (MQ=255)
                                                                                                                  cGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGc                                  <  8:164688/149‑1 (MQ=255)
                                                                                                                     aaGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGcaca                               >  2:388263/1‑149 (MQ=255)
                                                                                                                         cgcgTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCAt                           >  3:368230/1‑149 (MQ=255)
                                                                                                                                gTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAAt                    >  3:223259/1‑149 (MQ=255)
                                                                                                                                   ttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGc                  <  1:322527/148‑1 (MQ=255)
                                                                                                                                      tAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga              >  3:229244/1‑149 (MQ=255)
                                                                                                                                      tAGGCGCGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGCATCAGCGCGAAATGTGCCCGCGCCCTGCAAAGGCATGGGCCACGGAACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGCCACCCCATGGGCAAAGCgaga              >  3:316767/1‑149 (MQ=255)
                                                                                                                                       attcgcgGGTAAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt                                     <  3:350778/122‑1 (MQ=255)
                                                                                                                                       attcgcgGGTAAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGt                                     >  4:350778/4‑125 (MQ=255)
                                                                                                                                        ttcgcgGGTAAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACg            >  8:263583/3‑149 (MQ=255)
                                                                                                                                        ggCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGCAGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga              <  1:125278/148‑1 (MQ=255)
                                                                                                                                            cgGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGccac        <  1:79992/149‑1 (MQ=255)
                                                                                                                                               gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGccacc       >  7:168183/1‑147 (MQ=255)
                                                                                                                                                  tACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGcc  >  3:155258/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGCC  >  NC_002947/4851547‑4851840

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: