Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,549 G→C 5.5% R125P (CGC→CCC)  PP_3199 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,5490GC 5.5% 968.6 / 33.3 294R125P (CGC→CCC) PP_3199monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (119/155);  new base C (14/2);  total (133/161)
Fisher's exact test for biased strand distribution p-value = 5.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 305 total aligned reads displayed.

GGCCATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATC  >  NC_002947/3630400‑3630698
                                                                                                                                                     |                                                                                                                                                     
  ccATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGcc                                                                                                                                                     >  2:76955/1‑150 (MQ=255)
   cATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCt                                                                                                                                                    >  1:2344103/1‑150 (MQ=255)
    aTCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTg                                                                                                                                                   >  2:6900269/1‑150 (MQ=255)
     tCTATACCCCCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTgg                                                                                                                                                  <  2:6505737/150‑1 (MQ=255)
       tataCCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGCGCACCTgggt                                                                                                                                                <  1:4791305/150‑2 (MQ=255)
         taCCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTgggtgc                                                                                                                                              <  2:6077245/150‑4 (MQ=255)
         taCCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCAc                                                                                                                                              >  2:3003691/1‑150 (MQ=255)
         taCCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCAc                                                                                                                                              >  2:290448/1‑150 (MQ=255)
          accaccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCAcc                                                                                                                                             >  2:7175069/1‑150 (MQ=255)
           ccaccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCt                                                                                                                                            >  1:960860/1‑150 (MQ=255)
             accGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCg                                                                                                                                          >  2:5712890/1‑150 (MQ=255)
              ccGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCgg                                                                                                                                         <  2:1839883/150‑1 (MQ=255)
               cGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGt                                                                                                                                        <  1:4205227/150‑1 (MQ=255)
               cGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATACGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGt                                                                                                                                        >  1:1214177/1‑150 (MQ=255)
                gAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTg                                                                                                                                       <  1:845311/150‑1 (MQ=255)
                gAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTg                                                                                                                                       <  1:3208131/150‑1 (MQ=255)
                  aGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGa                                                                                                                                     <  1:4568992/150‑1 (MQ=255)
                   gggCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGaa                                                                                                                                    >  1:4892489/1‑150 (MQ=255)
                   gggCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGaa                                                                                                                                    >  1:8686235/1‑150 (MQ=255)
                     gcagcaACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAg                                                                                                                                  <  2:4773962/150‑1 (MQ=255)
                      cagcaACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGc                                                                                                                                 >  2:421923/1‑150 (MQ=255)
                      cagcaACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGACCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGc                                                                                                                                 >  1:9065727/1‑150 (MQ=255)
                       aggcaCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCCCAGGCGCACCtgg                                                                                                                                          >  2:4908759/5‑138 (MQ=255)
                         caACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                              <  1:7543589/150‑1 (MQ=255)
                         caACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTc                                                                                                                              >  2:3600353/1‑150 (MQ=255)
                          aaCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCg                                                                                                                             <  1:5330909/150‑1 (MQ=255)
                           aCTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCACGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGa                                                                                                                            <  2:7392206/150‑1 (MQ=255)
                                cgcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAg                                                                                                                       <  2:5269613/150‑1 (MQ=255)
                                 gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAga                                                                                                                      >  1:2273511/1‑150 (MQ=255)
                                 gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAga                                                                                                                      >  2:8707785/1‑150 (MQ=255)
                                 gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAga                                                                                                                      >  2:1981414/1‑150 (MQ=255)
                                 gcgcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAga                                                                                                                      >  2:2442428/1‑150 (MQ=255)
                                   gcATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAATTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAgacg                                                                                                                    >  2:6389897/1‑150 (MQ=255)
                                     aTCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAgacgac                                                                                                                  >  1:8967381/1‑150 (MQ=255)
                                     aTCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAgacgac                                                                                                                  >  1:3798082/1‑150 (MQ=255)
                                     aTCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAAGGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAACGCTTCGAGCACgacgac                                                                                                                  >  2:1649708/1‑150 (MQ=255)
                                      tCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACa                                                                                                                 <  1:8223371/150‑1 (MQ=255)
                                      tCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACa                                                                                                                 <  1:3249098/150‑1 (MQ=255)
                                         cAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAgcg                                                                                                              <  2:8541423/150‑1 (MQ=255)
                                         cAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAgcg                                                                                                              >  2:7844011/1‑150 (MQ=255)
                                           ggCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCa                                                                                                            <  2:877184/150‑1 (MQ=255)
                                           ggCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCa                                                                                                            <  1:1863241/150‑1 (MQ=255)
                                            gCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAg                                                                                                           >  1:6324234/1‑150 (MQ=255)
                                              cAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTg                                                                                                         <  2:8920432/150‑1 (MQ=255)
                                                gcgcCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTTGa                                                                                                       <  2:1884397/150‑1 (MQ=255)
                                                gcgcCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGa                                                                                                       <  1:3613639/150‑1 (MQ=255)
                                                 cgcCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGAc                                                                                                      <  1:7844011/150‑1 (MQ=255)
                                                 cgcCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGAc                                                                                                      >  2:3739467/1‑150 (MQ=255)
                                                  gcCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACg                                                                                                     >  1:6458800/1‑150 (MQ=255)
                                                  gcCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACg                                                                                                     >  1:3716723/1‑150 (MQ=255)
                                                       gCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTg                                                                                                <  1:8579427/150‑1 (MQ=255)
                                                        ccGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTg                                                                                                                                       <  1:5155178/110‑1 (MQ=255)
                                                            tCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCAc                                                                                           <  2:1405695/150‑1 (MQ=255)
                                                             ccGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCAcc                                                                                          <  1:3073188/150‑1 (MQ=255)
                                                              cGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCg                                                                                         >  2:2047710/1‑150 (MQ=255)
                                                               gCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGa                                                                                        >  1:5661928/1‑150 (MQ=255)
                                                               gCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGa                                                                                        >  1:5471453/1‑150 (MQ=255)
                                                                cAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGAt                                                                                       <  1:4673886/150‑1 (MQ=255)
                                                                    aTCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCa                                                                                   <  2:5898640/150‑1 (MQ=255)
                                                                             aTGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGgc                                                                          <  1:2013586/150‑1 (MQ=255)
                                                                                   gTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACg                                                                    <  1:6569633/150‑1 (MQ=255)
                                                                                      tctcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTc                                                                 >  1:1481331/1‑150 (MQ=255)
                                                                                       ctcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCt                                                                <  2:3475537/150‑1 (MQ=255)
                                                                                        tcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTg                                                               >  2:8593975/1‑150 (MQ=255)
                                                                                        tcGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTg                                                               >  1:6199414/1‑150 (MQ=255)
                                                                                              cgcgccCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTc                                                         <  1:1248171/150‑1 (MQ=255)
                                                                                                      gCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGa                                                 <  1:1002478/150‑1 (MQ=255)
                                                                                                       cacaAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGAc                                                >  1:4069490/1‑150 (MQ=255)
                                                                                                          aaGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAgcg                                                                                                              <  1:2359414/85‑1 (MQ=255)
                                                                                                          aaGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACggg                                             >  2:1302436/1‑150 (MQ=255)
                                                                                                          aaGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACggg                                             <  1:7268281/150‑1 (MQ=255)
                                                                                                           aGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGgt                                            <  1:269718/150‑1 (MQ=255)
                                                                                                            ggTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGgtt                                           <  1:7959678/150‑1 (MQ=255)
                                                                                                            ggTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGgtt                                           <  1:7434099/150‑1 (MQ=255)
                                                                                                              tgtTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGgttgt                                         >  1:4990249/1‑150 (MQ=255)
                                                                                                                  gAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTtctc                                     >  2:5584594/1‑150 (MQ=255)
                                                                                                                  gAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCCGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTtctc                                     <  1:6900269/150‑1 (MQ=255)
                                                                                                                    gCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGa                                   >  1:5247772/1‑150 (MQ=255)
                                                                                                                     cGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGaa                                  >  1:2553695/1‑150 (MQ=255)
                                                                                                                      gAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGaaa                                 <  1:521971/150‑1 (MQ=255)
                                                                                                                      gAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGaaa                                 <  2:9065727/150‑1 (MQ=255)
                                                                                                                      gAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGaaa                                 <  1:290448/150‑1 (MQ=255)
                                                                                                                      gAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGaaa                                 <  1:8288434/150‑1 (MQ=255)
                                                                                                                      gAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGATGTTCTCGaaa                                 >  1:2216115/1‑150 (MQ=255)
                                                                                                                              cATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGTCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAg                         >  1:2764240/1‑150 (MQ=255)
                                                                                                                              cATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAg                         >  1:5972104/1‑150 (MQ=255)
                                                                                                                               aTCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGc                        >  1:3937031/1‑150 (MQ=255)
                                                                                                                                tCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCg                       <  2:3920074/150‑1 (MQ=255)
                                                                                                                                    aCTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCt                   <  1:4344477/150‑1 (MQ=255)
                                                                                                                                    aCTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCt                   <  2:1019390/150‑1 (MQ=255)
                                                                                                                                       gggCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGtt                <  2:6199414/150‑1 (MQ=255)
                                                                                                                                           aCCGAGGTGCCCAGGCGCACCTGGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGa                                   <  2:8313428/127‑1 (MQ=255)
                                                                                                                                           aCCGAGGTGCCCAGGCGCACCTGGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGt            >  2:2605371/1‑150 (MQ=255)
                                                                                                                                           aCCGAGGTGCCCAGGCGCACCTGGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGt            >  2:4570754/1‑150 (MQ=255)
                                                                                                                                             ccAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGa          >  1:472943/1‑150 (MQ=255)
                                                                                                                                             ccAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACAGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGa          >  1:4012686/1‑150 (MQ=255)
                                                                                                                                                   gcgcCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTa    >  1:7712402/1‑150 (MQ=255)
                                                                                                                                                   gcgcCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTa    <  2:1214177/150‑1 (MQ=255)
                                                                                                                                                    cgcCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTAt   <  1:520304/150‑1 (MQ=255)
                                                                                                                                                     gcCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATc  >  1:3516367/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCATCTATACCACCGAAGGGCAGCAACTGGCGCGCATCCCAGGCCAGCGCCTGGCCGGTCCGCAGTATCCGGCCAATGTGGGTATCTCGCGCCGCGCCCTGCACAAGGTGTTGAGCGAAACTGCCATCGAACTGGGCACCCAGGTGCGCCTGGGCACCTCGGTGGAAAGCTTCGAGCAAGACGACAGCGCAGTGGACGTGCTGTTCACCGATGGCAGCCAGGGGCGCTACGGTCTGATGGTCGGCGCCGACGGGTTGTTCTCGAAAGTGCGCAGCGTGCTGTTCGGTGATCGTTATC  >  NC_002947/3630400‑3630698

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