Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,475,226 T→G 5.7% 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.7% 253.6 / 4.5 83D496A (GAC→GCC) PP_3087excinuclease ABC subunit A
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base G (3/2);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

CCCACCAGTTCCAGCAGGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTACA  >  NC_002947/3475082‑3475372
                                                                                                                                                 |                                                                                                                                                    
cccACCAGTTCCAGCAGGGCCTGGCTGACCAGGCTGTATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTcgcg                                                                                                                                                  <  6:104791/149‑1 (MQ=255)
            agcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCgggg                                                                                                                                     >  8:53433/1‑149 (MQ=255)
             ggagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCggggg                                                                                                                                    <  2:2016/147‑1 (MQ=255)
             gcagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCggggg                                                                                                                                    >  7:94498/1‑149 (MQ=255)
              cagGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGGTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCgggggg                                                                                                                                   >  6:5292/1‑149 (MQ=255)
                            ccAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTg                                                                                                                      >  5:393693/1‑148 (MQ=255)
                             cAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTa                                                                                                                    >  7:109727/1‑149 (MQ=255)
                             cAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTa                                                                                                                    >  7:266907/1‑149 (MQ=255)
                             cAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCATGTGGCCTGTa                                                                                                                    <  3:62794/149‑1 (MQ=255)
                              aGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGACGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTAc                                                                                                                   >  5:114085/1‑149 (MQ=255)
                               ggCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCt                                                                                                                       <  8:258167/144‑1 (MQ=255)
                               ggCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCt                                                                                                                       >  7:258167/1‑144 (MQ=255)
                               ggCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTAcc                                                                                                                  >  8:4622/1‑149 (MQ=255)
                                  tGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                               >  5:37244/1‑149 (MQ=255)
                                  tGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                               >  2:313392/1‑149 (MQ=255)
                                  tGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCAGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                               >  6:210634/1‑149 (MQ=255)
                                    gATTTGCCCGAGCCGGAACTACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGGTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGGCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                               <  3:45848/147‑1 (MQ=255)
                                    gATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAACTGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                               <  1:6930/147‑1 (MQ=255)
                                    gATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCggg                                                                                                               <  5:163804/147‑1 (MQ=255)
                                      tgtGCCCGAGCCGGAAATACCGGTTACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCCAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCAt                                                                                                          <  5:314837/147‑1 (MQ=255)
                                           ccGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCTAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGaa                                                                                                      <  1:309930/149‑1 (MQ=255)
                                           ccGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAAc                                                                                                     <  5:82003/149‑1 (MQ=255)
                                              aGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACagg                                                                                                  <  4:73283/149‑1 (MQ=255)
                                                  ggAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACagg                                                                                                  >  8:90634/1‑145 (MQ=255)
                                                  ggAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACagg                                                                                                  <  7:90634/145‑1 (MQ=255)
                                                         ccGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGt                                                                                       <  5:184707/149‑1 (MQ=255)
                                                          cGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACa                                                                                                    <  1:247765/135‑1 (MQ=255)
                                                          cGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACa                                                                                                    >  2:247765/1‑135 (MQ=255)
                                                          cGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTg                                                                                      >  1:377328/1‑149 (MQ=255)
                                                           ggTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACag                                                                                                   >  7:144818/1‑135 (MQ=255)
                                                           ggTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACag                                                                                                   <  8:144818/135‑1 (MQ=255)
                                                            gTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTGGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCa                                                                                    <  6:37244/149‑1 (MQ=255)
                                                             tAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCgg                                                                                                                >  5:12680/1‑121 (MQ=255)
                                                             tAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCgg                                                                                                                <  6:12680/121‑1 (MQ=255)
                                                             tAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGGGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCAt                                                                                   <  6:141163/149‑1 (MQ=255)
                                                               aCGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGGTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGa                                                                                 <  5:5292/149‑1 (MQ=255)
                                                                 ggcggTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACt                                                                               >  7:204705/1‑149 (MQ=255)
                                                                   cggTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGt                                                                             <  2:91960/149‑1 (MQ=255)
                                                                     gTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTg                                                                           >  7:90620/1‑149 (MQ=255)
                                                                            agcCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCg                                                                    <  1:230221/149‑1 (MQ=255)
                                                                               ccAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGc                                                                  >  5:23348/1‑148 (MQ=255)
                                                                                cAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGc                                                                  >  7:6287/1‑147 (MQ=255)
                                                                                 aaCGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAg                                                               <  8:94498/149‑1 (MQ=255)
                                                                                 aaCGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAg                                                               >  3:124665/1‑149 (MQ=255)
                                                                                 aaCGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGCGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAg                                                               <  7:53433/149‑1 (MQ=255)
                                                                                   cGGAAATGCGGCGCTGAGGTTATCGATGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGC                                                             <  7:4622‑M1/149‑1 (MQ=255)
                                                                                         tGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCcggcgg                                                       <  2:237836/149‑1 (MQ=255)
                                                                                                 tCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGccc                                                >  3:218679/1‑149 (MQ=255)
                                                                                                    ggTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTg                                            <  2:377328/149‑1 (MQ=255)
                                                                                                        aTCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTa                                          >  2:351000/1‑147 (MQ=255)
                                                                                                         tCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTagag                                       <  1:313392/149‑1 (MQ=255)
                                                                                                         tCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTaga                                        <  8:109727/148‑1 (MQ=255)
                                                                                                           gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAgg                                      <  8:266907/148‑1 (MQ=255)
                                                                                                           gAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGa                                     <  6:114085/149‑1 (MQ=255)
                                                                                                             ggTTCTTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATg                                   <  8:255966/149‑1 (MQ=255)
                                                                                                                 aTTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTg                                <  8:204705/148‑1 (MQ=255)
                                                                                                                   tGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTgccg                             >  4:301312/1‑149 (MQ=255)
                                                                                                                     cgcgTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTgccgcc                           >  7:282559/1‑149 (MQ=255)
                                                                                                                       cgTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCgt                         <  6:23348/149‑1 (MQ=255)
                                                                                                                         tGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCgtgt                       <  8:240763/149‑1 (MQ=255)
                                                                                                                          gATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGtt                      >  6:297742/1‑149 (MQ=255)
                                                                                                                           aTGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTc                     >  1:94568‑M2/1‑148 (MQ=255)
                                                                                                                             gCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGTGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCCCCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTcg                    <  1:112101‑M2/149‑3 (MQ=255)
                                                                                                                               cTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTcgccg                 <  2:77019‑M2/149‑6 (MQ=255)
                                                                                                                               cTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACATCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTcgccg                 <  3:166559‑M2/149‑6 (MQ=255)
                                                                                                                                ttCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTcgccgg                <  7:240854‑M2/149‑7 (MQ=255)
                                                                                                                                ttCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTcgcc                  >  6:94528‑M2/1‑143 (MQ=255)
                                                                                                                                ttCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTcgcc                  <  5:94528‑M2/147‑5 (MQ=255)
                                                                                                                                     aGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGc            <  8:90620/149‑1 (MQ=255)
                                                                                                                                      gCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGaacacggcggcgccatcctctacagcgggccgcccgccggcctgg           >  1:193230‑M1/1‑104 (MQ=255)
                                                                                                                                       cTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGaacacggcggcgccatcctctacagcgggccgcccgccggcctggc          <  2:193230‑M1/149‑47 (MQ=255)
                                                                                                                                        ccccgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTg                                                                           >  1:162089/5‑83 (MQ=38)
                                                                                                                                        ccccgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTg                                                                           <  2:162089/79‑1 (MQ=38)
                                                                                                                                        ccccgCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCggg         >  6:216246/5‑149 (MQ=255)
                                                                                                                                           aGCCAGTCGCGGGCCTGGCGGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGc        <  4:218679/148‑1 (MQ=255)
                                                                                                                                           aGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGccc      <  6:393693/149‑1 (MQ=255)
                                                                                                                                            gCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCAGTAGAGGATGGTGCCGCCGTCTTCGCCGGCAGCGGGccc      >  7:137686/1‑148 (MQ=255)
                                                                                                                                            gCCAGGCCCGCGACTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACATGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGccc      >  2:78138/1‑148 (MQ=255)
                                                                                                                                             ccAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTa    >  4:99873/1‑149 (MQ=255)
                                                                                                                                             ccAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTa    >  6:268289/1‑149 (MQ=255)
                                                                                                                                              cAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCGTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTa    >  8:6197/1‑148 (MQ=255)
                                                                                                                                              cAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTAc   <  5:210634/149‑1 (MQ=255)
                                                                                                                                               aGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTACa  >  4:337534/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCACCAGTTCCAGCAGGGCCTGGCTGACCAGGCTGGATTTGCCCGAGCCGGAAATACCGGTAACGGCGGTAAAGCAGCCCAACGGAAATGCGGCGCTCAGGTTATCGAGGTTATTGCGCGTGATGCCTTCAAGCTTCAGCCAGTCGCGGGCCTGGCGGGGGGCGCGTGTGGCCTGTACCGGGGCATCGAACAGGTAGGTGCGGGTGCATGACTGTTGCACCTCGGCCAGGCCGGCGGGCGGCCCGCTGTAGAGGATGGTGCCGCCGTGTTCGCCGGCAGCGGGCCCTACA  >  NC_002947/3475082‑3475372

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