Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,633,579 T→G 5.0% L126R (CTG→CGG)  PP_3199 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,633,5790TG 5.0% 1130.2 / 32.2 328L126R (CTG→CGG) PP_3199monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (150/161);  new base G (8/9);  total (158/170)
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.48e-01

Only 100 of 330 total aligned reads displayed.

CATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATCAGC  >  NC_002947/3633430‑3633728
                                                                                                                                                     |                                                                                                                                                     
cATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCt                                                                                                                                                       >  1:2845095/1‑150 (MQ=255)
  tCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTgg                                                                                                                                                     >  1:5791091/1‑150 (MQ=255)
  tCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTgg                                                                                                                                                     <  2:5243139/150‑1 (MQ=255)
  tCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTgg                                                                                                                                                     >  1:3163106/1‑150 (MQ=255)
   cTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTggg                                                                                                                                                    <  2:260447/150‑1 (MQ=255)
       accaccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCAcc                                                                                                                                                >  1:5068829/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:1195630/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:5242258/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:5302692/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              <  2:4197973/150‑1 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:2292755/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:1580479/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:1335008/1‑150 (MQ=255)
         caccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTc                                                                                                                                              >  1:2032315/1‑150 (MQ=255)
          accGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCg                                                                                                                                             <  1:1072760/150‑1 (MQ=255)
          accGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCCCAGGCGCACCtgg                                                                                                                                             <  2:4871718/150‑3 (MQ=255)
            cGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCGAGGTGCCCAGGCGCACCTGGGt                                                                                                                                           <  2:5939445/150‑1 (MQ=255)
            cGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGt                                                                                                                                           <  2:561370/150‑1 (MQ=255)
            cGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGt                                                                                                                                           <  2:2869561/150‑1 (MQ=255)
            cGAAGGGCAGCAACTGGCGCGCAGCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGt                                                                                                                                           <  1:4207791/150‑1 (MQ=255)
              aaGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTgg                                                                                                                                         >  1:1392426/1‑150 (MQ=255)
                gggCAGCAACTGGCGCGCAGCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGaa                                                                                                                                       <  1:6540669/150‑1 (MQ=255)
                   cagcaACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGaa                                                                                                                                       >  1:916972/1‑147 (MQ=255)
                    agcaACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTCTCCGGCCCATGTGGGTATCTCGCGCCGCGCCCTGCACACGGTGTTGAGCGAAACTGCCATCGAACGGGGCCCCCAGGTGCGCCTGGGCACCTCGGGGGAAAGCt                                                                                                                                   >  2:2437689/1‑150 (MQ=255)
                      caACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                                 <  2:6682029/150‑1 (MQ=255)
                      caACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                                 <  1:1660832/150‑1 (MQ=255)
                      caACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                                 >  2:3452315/1‑150 (MQ=255)
                        aCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCAGAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGa                                                                                                                               <  2:354215/150‑1 (MQ=255)
                        aCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGTGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGa                                                                                                                               <  2:3229633/150‑1 (MQ=255)
                        aCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGa                                                                                                                               <  2:3920930/150‑1 (MQ=255)
                        aCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGa                                                                                                                               <  1:6054790/150‑1 (MQ=255)
                        aCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGa                                                                                                                               >  1:5606386/1‑150 (MQ=255)
                              gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                                 <  1:332701/142‑1 (MQ=255)
                              gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                                 >  2:332522/1‑142 (MQ=255)
                              gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAga                                                                                                                         >  2:4747469/1‑150 (MQ=255)
                               cgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAgac                                                                                                                        <  1:191868/150‑1 (MQ=255)
                                  aTCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAgacgac                                                                                                                     >  2:3991989/1‑150 (MQ=255)
                                           cAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTg                                                                                                            <  1:2991146/150‑1 (MQ=255)
                                            aGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTgg                                                                                                           >  2:4536794/1‑150 (MQ=255)
                                                   ggCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCTCAGTGGACGTGCt                                                                                                    <  2:4514944/150‑1 (MQ=255)
                                                   ggCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCt                                                                                                    <  1:4849804/150‑1 (MQ=255)
                                                    gCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTg                                                                                                   >  2:5589877/1‑150 (MQ=255)
                                                    gCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTg                                                                                                   <  1:6722964/150‑1 (MQ=255)
                                                     ccGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGt                                                                                                  >  2:2788837/1‑150 (MQ=255)
                                                      cGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGtt                                                                                                 >  2:2667279/1‑150 (MQ=255)
                                                           cGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCg                                                                                            <  1:1809857/150‑1 (MQ=255)
                                                             cAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGAt                                                                                          <  2:3776926/150‑1 (MQ=255)
                                                              aGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGTAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATg                                                                                         <  1:5102154/150‑1 (MQ=255)
                                                               gTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATgg                                                                                        <  2:3580711/150‑1 (MQ=255)
                                                                 aTCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCa                                                                                      <  1:2618607/150‑1 (MQ=255)
                                                                   ccGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGc                                                                                    >  1:2308033/1‑150 (MQ=255)
                                                                     ggCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCa                                                                                  >  2:6514975/1‑150 (MQ=255)
                                                                      gCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAg                                                                                 >  1:254540/1‑150 (MQ=255)
                                                                       ccAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAgg                                                                                <  2:5109397/150‑1 (MQ=255)
                                                                         aaTGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAgggg                                                                              <  1:2471260/150‑1 (MQ=255)
                                                                          aTGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGgc                                                                             <  2:5302692/150‑1 (MQ=255)
                                                                            gtgGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGgcgc                                                                           >  1:5427667/1‑150 (MQ=255)
                                                                              gggTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTa                                                                         >  1:4980618/1‑150 (MQ=255)
                                                                              gggTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTa                                                                         >  1:2412809/1‑150 (MQ=255)
                                                                                   tctcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTc                                                                    >  2:1328279/1‑150 (MQ=255)
                                                                                    ctcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCt                                                                   <  1:3814276/150‑1 (MQ=255)
                                                                                     tcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTg                                                                  >  2:1815645/1‑150 (MQ=255)
                                                                                     tcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTg                                                                  >  2:1424977/1‑150 (MQ=255)
                                                                                     tcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTg                                                                  >  1:4789752/1‑150 (MQ=255)
                                                                                     tcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTg                                                                  >  2:4341037/1‑150 (MQ=255)
                                                                                        cgccgcgccCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATg                                                               <  2:685520/150‑1 (MQ=255)
                                                                                        cgccgcgccCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATg                                                               <  1:3252025/150‑1 (MQ=255)
                                                                                              gccCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGgc                                                         <  2:1683075/150‑1 (MQ=255)
                                                                                                  tGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCg                                                     >  1:6211327/1‑150 (MQ=255)
                                                                                                  tGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCg                                                     >  1:4475076/1‑150 (MQ=255)
                                                                                                   gCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGAc                                                                                                         <  2:6608068/97‑1 (MQ=255)
                                                                                                      caAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACATCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACgg                                                 <  2:2304402/150‑1 (MQ=255)
                                                                                                         ggTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGgtt                                              <  2:3058653/150‑1 (MQ=255)
                                                                                                            gtTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGgttgtt                                           <  1:4033829/150‑1 (MQ=255)
                                                                                                                   gAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGaaa                                    <  2:2002743/150‑1 (MQ=255)
                                                                                                                     aaCTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGt                                  <  1:843514/150‑1 (MQ=255)
                                                                                                                        tGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTGGTTCTCGAAAGTgcg                               >  2:5023484/1‑150 (MQ=255)
                                                                                                                          ccATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCa                             <  1:225819/150‑1 (MQ=255)
                                                                                                                            aTCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTCCGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGc                           >  1:1688887/1‑150 (MQ=255)
                                                                                                                            aTCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGc                           >  2:6245211/1‑150 (MQ=255)
                                                                                                                            aTCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGc                           >  1:2659360/1‑150 (MQ=255)
                                                                                                                            aTCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGc                           >  2:1006125/1‑150 (MQ=255)
                                                                                                                             tCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTCCGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCg                          >  2:1201529/1‑150 (MQ=255)
                                                                                                                             tCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCg                          >  2:2746603/1‑150 (MQ=255)
                                                                                                                                 aCTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCt                      >  2:2082470/1‑150 (MQ=255)
                                                                                                                                    gggCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGtt                   <  2:5791091/150‑1 (MQ=255)
                                                                                                                                      gCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCGAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCg                 <  2:6047394/150‑1 (MQ=255)
                                                                                                                                        aCCGAGGTGCCCAGGCGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGt               >  1:4168293/1‑150 (MQ=255)
                                                                                                                                        aCCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGt               >  1:6031931/1‑150 (MQ=255)
                                                                                                                                         cccAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTg              >  2:5039362/1‑150 (MQ=255)
                                                                                                                                          ccAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGGCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGa             >  1:3132875/1‑150 (MQ=255)
                                                                                                                                           cAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGAt            <  2:1590503/150‑1 (MQ=255)
                                                                                                                                            aGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATc           <  2:254540/150‑1 (MQ=255)
                                                                                                                                            aGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATc           <  2:4693248/150‑1 (MQ=255)
                                                                                                                                            aGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATc           <  2:4605465/150‑1 (MQ=255)
                                                                                                                                            aGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATc           >  2:2250703/1‑150 (MQ=255)
                                                                                                                                              gTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGt         <  1:4341037/150‑1 (MQ=255)
                                                                                                                                                 cgcCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTAt      <  1:6382186/150‑1 (MQ=255)
                                                                                                                                                  gcCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATc     <  1:4907075/150‑1 (MQ=255)
                                                                                                                                                     tGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATCAGc  <  1:3607536/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATCAGC  >  NC_002947/3633430‑3633728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑