Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,960 G→C 9.9% P121P (CCG→CCC 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,9600GC 9.9% 775.0 / 56.8 242P121P (CCG→CCCpydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base G (116/102);  new base C (5/19);  total (121/121)
Fisher's exact test for biased strand distribution p-value = 4.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 244 total aligned reads displayed.

GGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGC  >  NC_002947/4552811‑4553108
                                                                                                                                                     |                                                                                                                                                    
  gATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCt                                                                                                                                                    >  2:2809816/1‑150 (MQ=255)
   aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                                   <  1:5400961/150‑1 (MQ=255)
   aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                                   >  1:1787459/1‑150 (MQ=255)
   aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                                   <  2:5314896/150‑1 (MQ=255)
   aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                                   <  1:223006/150‑1 (MQ=255)
   aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                                                                   <  2:4940123/150‑1 (MQ=255)
    tCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtgg                                                                                                                                                  <  1:1463068/150‑1 (MQ=255)
     caacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                                                                                                                                 >  2:2601202/1‑150 (MQ=255)
       acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtgg                                                                                                                                               >  1:1696462/1‑150 (MQ=255)
       acaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtgg                                                                                                                                               >  1:3783897/1‑150 (MQ=255)
        caaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGa                                                                                                                                              >  1:4253096/1‑150 (MQ=255)
         aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                                                                                                             <  1:791928/150‑1 (MQ=255)
         aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                                                                                                             <  2:3726776/150‑1 (MQ=255)
         aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGACCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                                                                                                             <  1:3375539/150‑1 (MQ=255)
           tatCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAga                                                                                                                                           >  1:1060288/1‑149 (MQ=255)
           tatCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAga                                                                                                                                           <  2:4282461/150‑2 (MQ=255)
             tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                                                                                                         >  1:2385388/1‑150 (MQ=255)
             tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                                                                                                         <  1:713211/150‑4 (MQ=255)
             tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                                                                                                         >  2:1688722/1‑147 (MQ=255)
               gAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAACTGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaact                                                                                                                                       <  2:794715/150‑6 (MQ=255)
                aGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACgg                                                                                                                                      <  1:1359411/150‑1 (MQ=255)
                     aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                                 >  1:1877954/1‑150 (MQ=255)
                     aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                                 <  2:3315803/150‑1 (MQ=255)
                        accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                                                                              >  1:4367402/1‑150 (MQ=255)
                         ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                                                                                                             >  1:5334887/1‑150 (MQ=255)
                         ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                                                                                                             >  2:3283472/1‑150 (MQ=255)
                            accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa                                                                                                                          <  1:4708474/150‑1 (MQ=255)
                                ctctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                                                                              >  1:2958578/1‑142 (MQ=255)
                                 tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAg                                                                                                                     <  1:2645160/150‑1 (MQ=255)
                                      ggAAATCAACCTGCGTGAAATCACCCAGGAGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                                                                <  2:4805776/150‑1 (MQ=255)
                                        aaaTCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACt                                                                                                              >  2:4753903/1‑150 (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)
                                            gaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTgg                                                                                                          <  2:2222128/149‑1 (MQ=255)
                                              aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                                                                        <  1:2670924/150‑1 (MQ=255)
                                              aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                                                                        <  1:5394919/150‑1 (MQ=255)
                                              aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                                                                        >  2:4887380/1‑150 (MQ=255)
                                                cTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGc                                                                                                      <  2:3579270/150‑1 (MQ=255)
                                                 tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc                                                                                                     >  2:295681/1‑150 (MQ=255)
                                                     tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                                 >  2:3509034/1‑150 (MQ=255)
                                                     tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                                 >  1:5004626/1‑150 (MQ=255)
                                                       aaaTCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACg                                                                                               >  2:4648624/1‑150 (MQ=255)
                                                          tCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGAAGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTa                                                                                            <  1:5212771/150‑1 (MQ=255)
                                                             cccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGc                                                                                         >  2:2973288/1‑150 (MQ=255)
                                                                 ggTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGa                                                                                     <  2:3404814/150‑1 (MQ=255)
                                                                  gTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAg                                                                                    >  2:5117926/1‑150 (MQ=255)
                                                                        aagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGc                                                                              >  1:5163191/1‑150 (MQ=255)
                                                                            gcTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTGGCTGTCCGCACGGCATgagcgag                                                                                    <  2:4376028/139‑7 (MQ=255)
                                                                            gcTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTGGCTGTCCGCACGGCATGAGCGAGCGCGGCATgg                                                                          <  2:2460713/149‑1 (MQ=255)
                                                                            aCTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgg                                                                          <  1:2860088/150‑1 (MQ=255)
                                                                              tGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgggg                                                                        <  1:3118860/150‑1 (MQ=255)
                                                                               ggCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATggggg                                                                       <  2:3359241/150‑1 (MQ=255)
                                                                                 cccGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcg                                                                     <  1:4845544/150‑1 (MQ=255)
                                                                                 cccGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcg                                                                     >  1:4311976/1‑150 (MQ=255)
                                                                                    gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                                                                  <  1:547641/150‑1 (MQ=255)
                                                                                    gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                                                                  <  1:3768210/150‑1 (MQ=255)
                                                                                      ccGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGACGGAGCGCGGCATgg                                                                          <  1:3403462/140‑1 (MQ=255)
                                                                                       cgcgCACTGCTCGTGTCGCGGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTc                                                               <  2:4062715/150‑1 (MQ=255)
                                                                                       cgcgCACTGATCGTGTCGCTGATGGTGCCGTGCATAGTGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTc                                                               <  2:586568/150‑1 (MQ=255)
                                                                                        gcgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCg                                                              >  1:1180607/1‑150 (MQ=255)
                                                                                        gcgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCg                                                              >  1:939270/1‑150 (MQ=255)
                                                                                         cgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGCACCGCGGCATGGGGGCGGCCGTCgg                                                             >  1:3064694/1‑150 (MQ=255)
                                                                                            aCTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCa                                                          <  2:3142096/150‑1 (MQ=255)
                                                                                             cTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAg                                                         <  1:2880031/150‑1 (MQ=255)
                                                                                                   gtgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCg                                                   >  1:3602727/1‑150 (MQ=255)
                                                                                                      tCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGaa                                                >  1:3218737/1‑150 (MQ=255)
                                                                                                      tCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGaa                                                <  1:295406/150‑1 (MQ=255)
                                                                                                       cGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCgactc                                               <  2:2709953/150‑5 (MQ=255)
                                                                                                         cTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAatat                                             >  2:4175110/1‑150 (MQ=255)
                                                                                                         cTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAatat                                             >  1:2804716/1‑150 (MQ=255)
                                                                                                            aTGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgtg                                          <  2:83884/150‑1 (MQ=255)
                                                                                                             tGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTgg                                         <  2:1463529/150‑1 (MQ=255)
                                                                                                               gTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                                 >  2:2802030/1‑92 (MQ=255)
                                                                                                                 gCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAAt                                     <  1:4116750/150‑1 (MQ=255)
                                                                                                                 gCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAAt                                     <  2:1358278/150‑1 (MQ=255)
                                                                                                                  ccGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATg                                    >  1:2087581/1‑150 (MQ=255)
                                                                                                                   cGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATgg                                   <  1:4826681/150‑1 (MQ=255)
                                                                                                                         tAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCAccc                             <  2:4862212/150‑1 (MQ=255)
                                                                                                                          agagGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCg                            >  1:2535105/1‑150 (MQ=255)
                                                                                                                            agGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAGCTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCt                          >  1:490312/1‑150 (MQ=255)
                                                                                                                            agGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAatat                                             >  2:2366535/1‑131 (MQ=255)
                                                                                                                            agGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCt                          >  2:1033536/1‑150 (MQ=255)
                                                                                                                             ggAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTg                         >  1:4408271/1‑150 (MQ=255)
                                                                                                                              gAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTgg                        <  1:1035475/150‑1 (MQ=255)
                                                                                                                               aaTCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGt                       >  2:1124284/1‑150 (MQ=255)
                                                                                                                                   gTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCaa                   >  1:3400922/1‑150 (MQ=255)
                                                                                                                                   gTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCaa                   <  2:1099814/150‑1 (MQ=255)
                                                                                                                                   gTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCaa                   >  2:1894908/1‑150 (MQ=255)
                                                                                                                                     ggAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGa                 <  2:1066941/150‑1 (MQ=255)
                                                                                                                                      gAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGAc                >  1:3540350/1‑150 (MQ=255)
                                                                                                                                           tCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACt           >  1:491053/1‑150 (MQ=255)
                                                                                                                                             aTCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgc         >  1:3245203/1‑150 (MQ=255)
                                                                                                                                              tCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgcg        >  1:5062116/1‑150 (MQ=255)
                                                                                                                                              tCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgcg        >  1:2334304/1‑150 (MQ=255)
                                                                                                                                              tCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgc         >  2:1541151/1‑149 (MQ=255)
                                                                                                                                                cTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCAGGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgcgcg      >  2:2510967/1‑150 (MQ=255)
                                                                                                                                                 tGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgcgcgc     >  1:724854/1‑150 (MQ=255)
                                                                                                                                                 tGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgcgcgc     >  2:1635933/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGC  >  NC_002947/4552811‑4553108

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: