Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,475,226 T→G 5.8% D496A (GAC→GCC)  PP_3087 ← excinuclease ABC subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,475,2260TG 5.8% 216.4 / 3.5 69D496A (GAC→GCC) PP_3087excinuclease ABC subunit A
Reads supporting (aligned to +/- strand):  ref base T (28/37);  new base G (2/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TTCCAGCAGGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTAC  >  NC_002947/3475090‑3475371
                                                                                                                                        |                                                                                                                                                 
ttCCAGCAGGGCGTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGc                                                                                                                                       <  5:195024/149‑1 (MQ=255)
 tCCAGCTGGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGGCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGc                                                                                                                                       <  5:189444/148‑1 (MQ=255)
   cagcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCATCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCggg                                                                                                                                    <  2:64489/149‑1 (MQ=255)
    agcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCgggg                                                                                                                                   <  6:351732/149‑1 (MQ=255)
    agcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCgggg                                                                                                                                   >  2:16836/1‑149 (MQ=255)
     gcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCACCGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCgggg                                                                                                                                   >  3:94934/1‑148 (MQ=255)
     gcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCgggg                                                                                                                                   <  4:94934/148‑1 (MQ=255)
       agGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGgc                                                                                                                                >  4:62929/1‑149 (MQ=255)
       agGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGgc                                                                                                                                <  7:135597/149‑1 (MQ=255)
         ggCCTGGCTGACCAGGCTGGAATTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGgcgc                                                                                                                              >  6:161110/1‑149 (MQ=255)
               gCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCgtgt                                                                                                                          <  3:62929/147‑1 (MQ=255)
               gCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTgg                                                                                                                        <  1:58257/149‑1 (MQ=255)
               gCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTgg                                                                                                                        >  4:255366/1‑149 (MQ=255)
                cTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTACCGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGCGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGc                                                                                                                       <  5:288136/149‑1 (MQ=255)
                cTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGc                                                                                                                       <  6:151511/149‑1 (MQ=255)
                  gACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCt                                                                                                                     >  6:184101/1‑149 (MQ=255)
                  gACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCt                                                                                                                     >  5:205918/1‑149 (MQ=255)
                   aCCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGcc                                                                                                                      <  1:73844/147‑1 (MQ=255)
                      aGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTAc                                                                                                                 >  4:201284/1‑149 (MQ=255)
                        gCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTGATCGAGGTTAGTGCGCGTGATGCCGTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTCCCg                                                                                                               >  7:130065/1‑149 (MQ=255)
                            gATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                             <  1:16836/147‑1 (MQ=255)
                            gATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                             <  5:184101/147‑1 (MQ=255)
                              tttGCCAGAGCCGGAAATACCGCTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCAt                                                                                                         <  5:134931/149‑1 (MQ=255)
                                  cccGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGaa                                                                                                     >  1:81875/1‑149 (MQ=255)
                                  cccGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGCGGCGCGTGTGGCCTGTACCGGGGCATCGaa                                                                                                     <  1:31151/149‑1 (MQ=255)
                                   ccGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAACTGCGGCGCTCAGGTTATCGAGGGTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAAc                                                                                                    <  7:264635/149‑1 (MQ=255)
                                       gCCGGAAATTCCGGTATCGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCTTGGTTATTGCGCGTGATGGCTTCAAGCTTCAGCCAGTCGCGGGCGTGGCGGGGGGGGCGTGTGGCCTGTACCGGGGCATCGAACaggt                                                                                                <  4:75948/149‑1 (MQ=255)
                                       gCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACaggt                                                                                                <  2:81875/149‑1 (MQ=255)
                                           tcaaTACCGGGAACGGCGGTAAAGCAGCCCCACGGAAATGCGGCGCGCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACaggtaggt                                                                                            <  8:130065/147‑1 (MQ=255)
                                               tACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCgg                                                                                        <  6:205918/149‑1 (MQ=255)
                                                 ccGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGt                                                                                      <  3:255366/149‑1 (MQ=255)
                                                   ggTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTg                                                                                     >  2:48168/1‑148 (MQ=255)
                                                    gTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCa                                                                                   >  4:439/1‑149 (MQ=255)
                                                          gcggTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGTGGGTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTg                                                                             <  3:201284/149‑1 (MQ=255)
                                                                  gcagcCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCt                                                                     >  6:175985/1‑149 (MQ=255)
                                                                   cagcCCAACGGAAATGCGGCGCTCAGGTTAACGAGGTTAATGCCCGTGATGCCTTCAACCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTCTACCGGGGCATCGAACAGCTAGGTGCGGGTGCATGACTGATGCACCTc                                                                    >  2:79035/1‑149 (MQ=255)
                                                                    agcCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCg                                                                   <  5:161110/149‑1 (MQ=255)
                                                                         aaCGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGGGTGATGCCTTCAAGCTTCAGGCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAg                                                              <  1:79035/149‑1 (MQ=255)
                                                                          aCGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAgg                                                             <  3:439/149‑1 (MQ=255)
                                                                            ggAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGcc                                                           >  5:280013/1‑149 (MQ=255)
                                                                              aaaTGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGtgtg                                                                                                                         >  3:352025/1‑85 (MQ=255)
                                                                              aaaTGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGtgtg                                                                                                                         <  4:352025/85‑1 (MQ=255)
                                                                                         tCAGGTTATCGAGGTTATTGCGCGTGAGGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGcc                                                <  2:302063/147‑1 (MQ=255)
                                                                                            ggTTATCGAGGTTATTGCGCGGGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTg                                           <  6:77119/149‑1 (MQ=255)
                                                                                                   gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACt                                                                              <  6:341328/107‑1 (MQ=255)
                                                                                                   gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACt                                                                              >  5:341328/1‑107 (MQ=255)
                                                                                                   gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCggcgggc                                                    >  7:35798/1‑133 (MQ=255)
                                                                                                   gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCggcgggc                                                    <  8:35798/133‑1 (MQ=255)
                                                                                                   gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGa                                    >  5:67201/1‑149 (MQ=255)
                                                                                                         aTTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTgc                              >  6:192202/1‑149 (MQ=255)
                                                                                                          ttGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTgcc                             >  2:226209/1‑149 (MQ=255)
                                                                                                          ttGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTgc                              <  5:192202/148‑1 (MQ=255)
                                                                                                             cgcgTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTgccgc                           <  6:280013/148‑1 (MQ=255)
                                                                                                               cgTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCgt                        <  1:234571/149‑1 (MQ=255)
                                                                                                                       cTTCAAGCTTCAGACAGTGGGGGGCCTGTAGGTGGGCGCGTGTGTCCTGTTCCGGCGAATCTAAATGGTAGGTGCGGGTTCATGACTGTTGCACCTCTGCCAGTCCGGCTGGCGTCCAGCTGTAGAGGATGGTGCCGCCGTGTTCGCCg                <  8:305068/149‑1 (MQ=255)
                                                                                                                        ttCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCgg               >  6:296240/1‑149 (MQ=255)
                                                                                                                          cAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCgg               >  2:168202/1‑147 (MQ=255)
                                                                                                                           aaGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCaa            >  4:64185/1‑148 (MQ=255)
                                                                                                                             gCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTc                    >  6:245496/1‑139 (MQ=255)
                                                                                                                             gCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTc                    <  5:245496/139‑1 (MQ=255)
                                                                                                                              cTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACATCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCgg         >  4:284517/1‑149 (MQ=255)
                                                                                                                                 ccgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCg          <  5:101306/143‑1 (MQ=255)
                                                                                                                                 ccgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCg          >  6:101306/3‑145 (MQ=255)
                                                                                                                                 ccgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGcc      >  7:104136/3‑149 (MQ=255)
                                                                                                                                  cgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGccc     <  8:104136/148‑1 (MQ=255)
                                                                                                                                  aGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCAGCCGTGTTAGCCGGCAGCGGGcc      >  3:20856/1‑148 (MQ=255)
                                                                                                                                   gCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTATACCGGGGCATCGAACAGGTAGGTACGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGccc     >  1:258818/1‑148 (MQ=255)
                                                                                                                                    ccAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTa   <  1:48168/149‑1 (MQ=255)
                                                                                                                                     cAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTAc  <  5:204271/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
TTCCAGCAGGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTAC  >  NC_002947/3475090‑3475371

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