Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,419,630 G→C 6.5% E106Q (GAG→CAG)  PP_4759 → transcriptional regulator, GntR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,419,6300GC 6.5% 200.0 / 3.9 61E106Q (GAG→CAG) PP_4759transcriptional regulator, GntR family
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base C (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GAAGCTGCAGGCCTCTGGGCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACCTGCTGATTGCC  >  NC_002947/5419488‑5419776
                                                                                                                                              |                                                                                                                                                  
gAAGCTGCAGGCCTCTGGGCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGc                                                                                                                                               >  5:31430/1‑148 (MQ=255)
      gCAGGCCTCTGGGCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTg                                                                                                                                        >  6:80903/1‑149 (MQ=255)
       cAGGCCTCTGGGCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGc                                                                                                                                       <  3:152595/149‑1 (MQ=255)
                 ggCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGc                                                                                                                              <  6:97482/148‑1 (MQ=255)
                 ggCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTCAGGGCCAGGCTGCCGCGCTGGc                                                                                                                              <  7:60402/148‑1 (MQ=255)
                    tggtggAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCCTCAAGGCCaat                                                                                                                                           >  6:35692/1‑130 (MQ=255)
                    tggtggAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCCTCAAGGCCaat                                                                                                                                           <  5:35692/132‑3 (MQ=255)
                         gAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAg                                                                                                                                            >  5:310049/1‑126 (MQ=255)
                         gAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAg                                                                                                                                            <  6:310049/126‑1 (MQ=255)
                         gAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcg                                                                                                                     >  2:380642/1‑149 (MQ=255)
                          aGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcgt                                                                                                                    <  8:234527/149‑1 (MQ=255)
                              gcgcCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcgtcgt                                                                                                                 >  7:75938/1‑148 (MQ=255)
                                gcCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGAc                                                                                                              >  7:187405/1‑149 (MQ=255)
                                  cACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCTGCGCTGGCCGCACTGCGTCGTGACGa                                                                                                            <  5:80903/149‑1 (MQ=255)
                                   aCGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcgtcgt                                                                                                                 >  8:223579/1‑143 (MQ=255)
                                   aCGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcgtcgt                                                                                                                 <  7:223579/143‑1 (MQ=255)
                                       gATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAgggg                                                                                                        >  2:68257/1‑148 (MQ=255)
                                       gATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAgggg                                                                                                        >  8:366840/1‑148 (MQ=255)
                                            gTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAggg                                                                                                         >  2:306994/1‑142 (MQ=255)
                                            gTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAggg                                                                                                         <  1:306994/142‑1 (MQ=255)
                                            gTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTg                                                                                                  >  1:58621/1‑149 (MQ=255)
                                             tACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTgg                                                                                                 >  4:138509/1‑149 (MQ=255)
                                                    gTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCAt                                                                                          >  6:48133/1‑149 (MQ=255)
                                                         ggAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTAGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGcc                                                                                     >  1:375386/1‑149 (MQ=255)
                                                          gAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGcc                                                                                     <  6:31430/148‑1 (MQ=255)
                                                           aaCGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAg                                                                                   >  5:280191/1‑149 (MQ=255)
                                                             ggCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGgc                                                                                 <  3:315365/148‑1 (MQ=255)
                                                             cGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGgc                                                                                 >  3:275123/1‑149 (MQ=255)
                                                             cGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCAGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGgc                                                                                 >  7:167838/1‑149 (MQ=255)
                                                                   gcccaggccGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAggtcc                                                                               <  5:333661/145‑4 (MQ=255)
                                                                   gcccaggccGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAggtcc                                                                               >  6:333661/1‑142 (MQ=255)
                                                                    cccaggccGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcgtcg                                                                                                                  >  7:272695/1‑109 (MQ=255)
                                                                    cccaggccGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGcgtcg                                                                                                                  <  8:272695/109‑1 (MQ=255)
                                                                      caggccGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGa                                                                        >  5:338868/1‑149 (MQ=255)
                                                                           cggccTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATc                                                                   <  1:68257/149‑1 (MQ=255)
                                                                           cggccTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATc                                                                   <  2:58621/149‑1 (MQ=255)
                                                                           cggccTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATc                                                                   <  5:48133/149‑1 (MQ=255)
                                                                                tgagcgTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGAc                                                              >  6:193768/4‑149 (MQ=255)
                                                                                   gcgTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTAGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTg                                                           <  2:375386/149‑1 (MQ=255)
                                                                                      ttGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGcc                                                        <  6:338868/149‑1 (MQ=255)
                                                                                        gcggcgGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGc                                                      >  2:404423/1‑149 (MQ=255)
                                                                                          ggcggAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTg                                                    >  1:173526/1‑149 (MQ=255)
                                                                                           gcggAGAGCGCAGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGc                                                   <  8:167838/149‑1 (MQ=255)
                                                                                              gagaGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCgg                                                <  3:138509/149‑1 (MQ=255)
                                                                                              gagaGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCgg                                                <  1:380642/149‑1 (MQ=255)
                                                                                                 agCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGc                                                   <  3:5715/143‑1 (MQ=255)
                                                                                                 agCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGc                                                   >  4:5715/1‑143 (MQ=255)
                                                                                                      ggcgAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGctgcctgc                                               <  2:292539/142‑1 (MQ=255)
                                                                                                      ggcgAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGc                                               >  1:292539/1‑142 (MQ=255)
                                                                                                              tgCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCg                                  >  8:205861/1‑147 (MQ=255)
                                                                                                              tgCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGaa                                <  1:404423/149‑1 (MQ=255)
                                                                                                               gcgcGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCAt                                    <  4:226463/144‑1 (MQ=255)
                                                                                                               gcgcGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCAt                                    >  3:226463/1‑144 (MQ=255)
                                                                                                                          tgGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGgcgc                    <  5:208650/149‑1 (MQ=255)
                                                                                                                                 gCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCAcc             >  4:18227/1‑149 (MQ=255)
                                                                                                                                   agcctGGCCCTCAAGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACctg           >  8:250942/5‑149 (MQ=255)
                                                                                                                                     aTTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACctgct         >  2:75321/1‑149 (MQ=255)
                                                                                                                                         gCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACCTGCTGAt      <  1:75321/148‑1 (MQ=255)
                                                                                                                                           cTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACCTGCTGATTGc   >  8:25456/1‑149 (MQ=255)
                                                                                                                                           cTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGAACTGGCCGCTGCCGGCGACAGCTGCATCGCAGCCGACCGGCGCTTTCACCTGCTCATTGc   >  3:286247/1‑149 (MQ=255)
                                                                                                                                            ttGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACCTGCTGATTGcc  >  4:337496/1‑149 (MQ=255)
                                                                                                                                              gAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGaa                                >  1:403406/1‑117 (MQ=255)
                                                                                                                                              gAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGc                              >  2:294062/1‑119 (MQ=255)
                                                                                                                                              gAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGc                              <  1:294062/119‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GAAGCTGCAGGCCTCTGGGCTGGTGGAGCCGCGCCACGGGATCGGTACGTTCGTGATGGAACGCCAGGCCCAGGCCGGCCTGCGCGTTGCGGCGGAGAGCGCGGCGAATGTGCGCGATCTGCTGGAGCTGCGTATTGGCCTTGAGGGCCAGGCTGCCGCGCTGGCCGCACTGCGTCGTGACGAGGGGCATCTGGCGCGCATGCGCCAGGCGCTGGATGACTATCAGGACCTGGCCGCTGCCGGCGACAGCTGCATCGAAGCCGACCGGCGCTTTCACCTGCTGATTGCC  >  NC_002947/5419488‑5419776

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