Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,949 T→G 10.7% F118V (TTC→GTC) ‡ pydA → NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,552,9490TG10.7% 678.4 / 56.0 231F118V (TTC→GTC) ‡pydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base T (107/99);  new base G (5/20);  total (112/119)
Fisher's exact test for biased strand distribution p-value = 2.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 231 total aligned reads displayed.

CAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGAC  >  NC_002947/4552802‑4553094
                                                                                                                                                   |                                                                                                                                                 
cAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTc                                                                                                                                                 <  2:5379538/150‑1 (MQ=255)
 aGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCa                                                                                                                                                <  1:3414651/150‑1 (MQ=255)
 aGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCa                                                                                                                                                <  1:4098502/150‑1 (MQ=255)
  ggTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCAt                                                                                                                                               <  2:4780305/150‑1 (MQ=255)
    ggtcggtgATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATcc                                                                                                                                             >  1:4653401/8‑150 (MQ=255)
     aCAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCt                                                                                                                                            >  1:350278/1‑150 (MQ=255)
      cAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTg                                                                                                                                           <  1:815811/150‑1 (MQ=255)
        tggtctCAACAATATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGcc                                                                                                                                         <  2:3780113/145‑1 (MQ=255)
         gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCg                                                                                                                                        >  1:338798/1‑150 (MQ=255)
            aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                     <  2:4939755/150‑1 (MQ=255)
            aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                     <  1:5400961/150‑1 (MQ=255)
              caacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                                                                                                                   >  2:5193066/1‑150 (MQ=255)
              caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                                                                                                                                   <  2:2477013/150‑1 (MQ=255)
              caacaaGATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                                                                                                                   <  1:74727/150‑1 (MQ=255)
               aacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGTAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtg                                                                                                                                  >  2:1593466/1‑150 (MQ=255)
               aacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtg                                                                                                                                  <  2:2599013/150‑1 (MQ=255)
                acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtgg                                                                                                                                 >  1:1696462/1‑150 (MQ=255)
                  aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                                                                                               <  2:3345090/150‑1 (MQ=255)
                    tatCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAga                                                                                                                             <  2:4282461/150‑2 (MQ=255)
                    tatCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAga                                                                                                                             >  1:1060288/1‑149 (MQ=255)
                      tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                                                                                           <  1:713211/150‑4 (MQ=255)
                        gAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAACTGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaact                                                                                                                         <  2:794715/150‑6 (MQ=255)
                         aGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACgg                                                                                                                        <  1:1359411/150‑1 (MQ=255)
                              cccACCGACCGCTCACTGGAAACCACCCTGCGTGAAATCACCCCGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                   <  1:3822131/148‑1 (MQ=255)
                              aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                   >  1:1877954/1‑150 (MQ=255)
                              aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                   >  2:1653659/1‑150 (MQ=255)
                              aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                   <  2:3315803/150‑1 (MQ=255)
                               tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg                                                                                                                  >  1:2979342/1‑150 (MQ=255)
                               tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg                                                                                                                  >  1:3682316/1‑150 (MQ=255)
                                 accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTc                                                                                                                                                 >  2:4713922/1‑117 (MQ=255)
                                 accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                                                                >  1:4367402/1‑150 (MQ=255)
                                  ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                                                                                               <  2:2304102/150‑1 (MQ=255)
                                  ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                                                                                               >  1:5334887/1‑150 (MQ=255)
                                    gaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGc                                                                                                             <  1:5193066/150‑1 (MQ=255)
                                     accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa                                                                                                            <  1:5055433/150‑1 (MQ=255)
                                     accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa                                                                                                            <  1:4708720/150‑1 (MQ=255)
                                         ctctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                                                                <  2:2958578/142‑1 (MQ=255)
                                            tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                                                                                                     >  1:2223093/1‑150 (MQ=255)
                                            tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                                                                                                     <  1:4204836/150‑1 (MQ=255)
                                            tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                                                                                                     <  1:2809816/150‑1 (MQ=255)
                                               ggAAATCAACCTGCGTGAAATCACCCAGGAGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                                                  <  2:4805776/150‑1 (MQ=255)
                                                  agtgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACtt                                                                                               <  1:1688722/146‑1 (MQ=255)
                                                   gtgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACttt                                                                                              <  1:4548384/147‑2 (MQ=255)
                                                    tgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACtttg                                                                                             <  2:3783897/148‑3 (MQ=255)
                                                    tCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCg                                                                                             <  2:490607/150‑1 (MQ=255)
                                                       aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                                                          <  1:5394919/150‑1 (MQ=255)
                                                       aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                                                          >  2:4887380/1‑150 (MQ=255)
                                                       aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                                                          <  1:2670924/150‑1 (MQ=255)
                                                          tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc                                                                                       >  2:295406/1‑150 (MQ=255)
                                                           gCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGccc                                                                                      >  2:3082282/1‑150 (MQ=255)
                                                              tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                   >  1:5004626/1‑150 (MQ=255)
                                                              tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                   >  2:3508516/1‑150 (MQ=255)
                                                                aaaTCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACg                                                                                 >  2:4648624/1‑150 (MQ=255)
                                                                       ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGcc                                                                          >  2:1085646/1‑150 (MQ=255)
                                                                           gTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAg                                                                      >  2:5117926/1‑150 (MQ=255)
                                                                               agaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGcgcg                                                                  >  1:1725236/1‑150 (MQ=255)
                                                                               agaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGcgcg                                                                  <  2:2121746/150‑1 (MQ=255)
                                                                                 aagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGc                                                                >  1:5163191/1‑150 (MQ=255)
                                                                                 aagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGc                                                                >  1:2434760/1‑150 (MQ=255)
                                                                                 aagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGAAGCCAGAACGGAATGCCGGAGCGCGGc                                                                >  2:1012923/1‑150 (MQ=255)
                                                                                  tgcgcTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTGGCTGTCCGCACGGCATGAGCGAGCGCGGCa                                                               <  1:2934831/146‑1 (MQ=255)
                                                                                   ggACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCAt                                                              >  2:736596/1‑150 (MQ=255)
                                                                                     gcTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTGGCTGTCCGCACGGCATgagcgag                                                                      <  2:4376028/139‑7 (MQ=255)
                                                                                       tGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgggg                                                          <  1:3118860/150‑1 (MQ=255)
                                                                                          cccGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcg                                                       <  1:4845544/150‑1 (MQ=255)
                                                                                             gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                                                    <  1:3768210/150‑1 (MQ=255)
                                                                                             gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                                                    <  1:547641/150‑1 (MQ=255)
                                                                                               ccGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGt                                                  <  1:4853889/150‑1 (MQ=255)
                                                                                               ccGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGACGGAGCGCGGCATgg                                                            >  2:3403462/1‑140 (MQ=255)
                                                                                                cgcgCACTGCTCGTGTCGCGGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTc                                                 <  2:4062715/150‑1 (MQ=255)
                                                                                                 gcgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCg                                                >  1:939270/1‑150 (MQ=255)
                                                                                                 gcgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCg                                                >  1:1180649/1‑150 (MQ=255)
                                                                                                     aCTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCa                                            <  2:3142096/150‑1 (MQ=255)
                                                                                                       tGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAgg                                          >  1:3199149/1‑150 (MQ=255)
                                                                                                            gtgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCg                                     >  1:4324750/1‑150 (MQ=255)
                                                                                                             tgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcggc                                                     <  2:276035/133‑1 (MQ=255)
                                                                                                             tgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCgg                                    >  2:170100/1‑150 (MQ=255)
                                                                                                               tCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGaa                                  <  1:295681/150‑1 (MQ=255)
                                                                                                               tCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGaa                                  >  1:3218870/1‑150 (MQ=255)
                                                                                                                cGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCgactc                                 <  2:2709953/150‑5 (MQ=255)
                                                                                                                  cTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAatat                               >  2:4175110/1‑150 (MQ=255)
                                                                                                                   tGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtg                              <  2:3649871/150‑1 (MQ=255)
                                                                                                                    gATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgt                             >  1:3401979/1‑150 (MQ=255)
                                                                                                                     aTGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgtg                            <  2:83884/150‑1 (MQ=255)
                                                                                                                      tGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTgg                           <  2:1463529/150‑1 (MQ=255)
                                                                                                                      tGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTgg                           <  1:28724/150‑1 (MQ=255)
                                                                                                                        gTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                   <  1:2802030/92‑1 (MQ=255)
                                                                                                                          gCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAAt                       <  2:1358278/150‑1 (MQ=255)
                                                                                                                           ccGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATg                      >  1:2087581/1‑150 (MQ=255)
                                                                                                                            cGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATgg                     <  1:4826681/150‑1 (MQ=255)
                                                                                                                                     agGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAatat                               <  1:2366535/131‑1 (MQ=255)
                                                                                                                                      ggAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTg           <  2:2385388/150‑1 (MQ=255)
                                                                                                                                      ggAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTg           >  1:4408271/1‑150 (MQ=255)
                                                                                                                                       gAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTgg          <  1:1035475/150‑1 (MQ=255)
                                                                                                                                        aaTCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGt         <  2:2804716/150‑1 (MQ=255)
                                                                                                                                        aaTCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGt         >  2:1124284/1‑150 (MQ=255)
                                                                                                                                            gTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCaa     >  1:2564084/1‑150 (MQ=255)
                                                                                                                                              ggAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGa   <  2:3566597/150‑1 (MQ=255)
                                                                                                                                              ggAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGa   <  2:1066941/150‑1 (MQ=255)
                                                                                                                                               gAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGAc  >  1:3540615/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGAC  >  NC_002947/4552802‑4553094

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: