Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,911 T→C 12.5% V105A (GTG→GCG) ‡ 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,9110TC12.5% 584.6 / 56.5 216V105A (GTG→GCG) ‡pydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base T (90/99);  new base C (6/21);  total (96/120)
Fisher's exact test for biased strand distribution p-value = 1.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 224 total aligned reads displayed.

TGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGG  >  NC_002947/4552764‑4553060
                                                                                                                                                   |                                                                                                                                                     
  tcgtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCg                                                                                                                                                   <  2:1431876/150‑1 (MQ=255)
   cgtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGc                                                                                                                                                  <  1:5418521/150‑1 (MQ=255)
   cgtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGc                                                                                                                                                  <  2:3232166/150‑1 (MQ=255)
    gtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCt                                                                                                                                                 >  1:3142096/1‑150 (MQ=255)
      cgcgCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGa                                                                                                                                               <  1:3019374/150‑1 (MQ=255)
           tACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTg                                                                                                                                          >  2:447088/1‑150 (MQ=255)
             cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGcc                                                                                                                                        <  1:485005/150‑1 (MQ=255)
              tCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCg                                                                                                                                       >  2:4283227/1‑150 (MQ=255)
                ggCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTg                                                                                                                                     >  1:83884/1‑150 (MQ=255)
                  cgcACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCa                                                                                                                                   <  2:363915/150‑1 (MQ=255)
                   gcACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCAt                                                                                                                                  >  2:537643/1‑150 (MQ=255)
                      cTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATaga                                                                                                                               <  2:3945429/150‑1 (MQ=255)
                       tATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                                                                              <  2:5256139/150‑1 (MQ=255)
                            ccccAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATc                                                                                                                         <  1:1774724/150‑1 (MQ=255)
                              ccAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGt                                                                                                                       >  1:2709953/1‑150 (MQ=255)
                                  ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                                                                   <  2:4218924/150‑1 (MQ=255)
                                  ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGANCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                                                                   <  2:4219044/150‑1 (MQ=255)
                                  ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGCACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                                                                   <  2:3067076/150‑1 (MQ=255)
                                   cGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAg                                                                                                                  >  1:4927021/1‑150 (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)
                                              tggtctCAACAATATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGcc                                                                                                       <  2:3780113/145‑1 (MQ=255)
                                                 gATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCt                                                                                                    >  2:2809816/1‑150 (MQ=255)
                                                  aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                   <  2:4940123/150‑1 (MQ=255)
                                                  aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                                                                   <  1:223006/150‑1 (MQ=255)
                                                   tCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtgg                                                                                                  <  1:1463068/150‑1 (MQ=255)
                                                    caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                                                                                                 <  2:2477013/150‑1 (MQ=255)
                                                    caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                                                                                                 <  1:372316/150‑1 (MQ=255)
                                                    caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                                                                                                 <  1:3101817/150‑1 (MQ=255)
                                                     aacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtg                                                                                                <  2:2599013/150‑1 (MQ=255)
                                                      acaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtgg                                                                                               >  1:3783897/1‑150 (MQ=255)
                                                       caaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGa                                                                                              >  1:4253096/1‑150 (MQ=255)
                                                        aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                                                             <  1:791928/150‑1 (MQ=255)
                                                        aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGACCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                                                             <  1:3375539/150‑1 (MQ=255)
                                                          tatCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAga                                                                                           >  1:1060288/1‑149 (MQ=255)
                                                            tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                                                         >  2:1195446/1‑150 (MQ=255)
                                                            tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                                                         >  2:3473710/1‑150 (MQ=255)
                                                            tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                                                         >  1:2385388/1‑150 (MQ=255)
                                                            tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                                                         >  2:1688722/1‑147 (MQ=255)
                                                            tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                                                         <  1:713479/150‑4 (MQ=255)
                                                              gAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACg                                                                                       >  2:4597525/1‑150 (MQ=255)
                                                                    cccACCGACCGCTCACTGGAAACCACCCTGCGTGAAATCACCCCGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                 <  1:3822131/148‑1 (MQ=255)
                                                                    aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                 <  2:3315803/150‑1 (MQ=255)
                                                                    aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                 <  2:5129264/150‑1 (MQ=255)
                                                                     tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg                                                                                >  1:3682316/1‑150 (MQ=255)
                                                                       accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTc                                                                                                               <  1:4713922/117‑1 (MQ=255)
                                                                       accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                              >  1:4367402/1‑150 (MQ=255)
                                                                       accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                              <  1:3415609/150‑1 (MQ=255)
                                                                        ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                                                             >  1:5334887/1‑150 (MQ=255)
                                                                        ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                                                             >  2:3283472/1‑150 (MQ=255)
                                                                         cgaccgGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCgaagaagaagc                                                                                                                         >  1:1571847/1‑103 (MQ=255)
                                                                         cgaccgGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCgaagaagaagc                                                                                                                         <  2:1571847/105‑3 (MQ=255)
                                                                          gaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGc                                                                           <  1:5193066/150‑1 (MQ=255)
                                                                           accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa                                                                          <  1:4708720/150‑1 (MQ=255)
                                                                             cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                                                                        >  1:705228/1‑150 (MQ=255)
                                                                               ctctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                              <  2:2958578/142‑1 (MQ=255)
                                                                               ctctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                              >  1:2958578/1‑142 (MQ=255)
                                                                                tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAg                                                                     <  1:2645188/150‑1 (MQ=255)
                                                                                  tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                                                                              <  2:1040066/91‑1 (MQ=255)
                                                                                  tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                                                                              >  1:1040066/1‑91 (MQ=255)
                                                                                  tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                                                                   <  1:2809816/150‑1 (MQ=255)
                                                                                   ctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTg                                                                  >  2:2335581/1‑150 (MQ=255)
                                                                                    tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                                                                   <  2:752346/100‑1 (MQ=255)
                                                                                     ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                <  1:2072282/150‑1 (MQ=255)
                                                                                     ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                <  1:5110894/150‑1 (MQ=255)
                                                                                     ggAAATCAACCTGCGTGAAATCACCCAGGAGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                <  2:4805776/150‑1 (MQ=255)
                                                                                         gtgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACttt                                                            <  1:4548384/147‑2 (MQ=255)
                                                                                          tgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACtttg                                                           <  1:1713874/148‑3 (MQ=255)
                                                                                          tCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCg                                                           <  2:490607/150‑1 (MQ=255)
                                                                                           cAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCgg                                                          >  2:838120/1‑150 (MQ=255)
                                                                                            aaCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCg                                                                                                                        >  2:3695506/1‑87 (MQ=255)
                                                                                             aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                        >  2:4887380/1‑150 (MQ=255)
                                                                                             aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                                                        <  1:5394919/150‑1 (MQ=255)
                                                                                                tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc                                                     >  2:295406/1‑150 (MQ=255)
                                                                                                    tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                 >  2:3509034/1‑150 (MQ=255)
                                                                                                    tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                 >  2:1443787/1‑150 (MQ=255)
                                                                                                     gAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCAc                                                >  1:1691202/1‑150 (MQ=255)
                                                                                                             ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                                                                   >  1:5038363/1‑75 (MQ=255)
                                                                                                             ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggggg                                                                                   <  2:4710457/107‑1 (MQ=255)
                                                                                                             ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGcc                                        <  2:1613108/150‑1 (MQ=255)
                                                                                                                 gTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAg                                    >  2:5117746/1‑150 (MQ=255)
                                                                                                                     agaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGcgcg                                <  2:2121746/150‑1 (MQ=255)
                                                                                                                     agaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGcgcg                                >  2:4104068/1‑150 (MQ=255)
                                                                                                                       aagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGc                              <  2:2637085/150‑1 (MQ=255)
                                                                                                                       aagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGAAGCCAGAACGGAATGCCGGAGCGCGGc                              >  2:1012923/1‑150 (MQ=255)
                                                                                                                              ggCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATggggg                       <  2:3359241/150‑1 (MQ=255)
                                                                                                                                cccGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcg                     <  1:2108947/150‑1 (MQ=255)
                                                                                                                                   gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                  <  1:547353/150‑1 (MQ=255)
                                                                                                                                   gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                  <  1:3768210/150‑1 (MQ=255)
                                                                                                                                     ccGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGt                <  1:4853889/150‑1 (MQ=255)
                                                                                                                                     ccGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGACGGAGCGCGGCATgg                          <  1:3403462/140‑1 (MQ=255)
                                                                                                                                      cgcgCACTGATCGTGTCGCTGATGGTGCCGTGCATAGTGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTc               <  2:586568/150‑1 (MQ=255)
                                                                                                                                       gcgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCg              >  1:939270/1‑150 (MQ=255)
                                                                                                                                       gcgcACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCg              >  1:1180649/1‑150 (MQ=255)
                                                                                                                                            cTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAg         <  1:2880031/150‑1 (MQ=255)
                                                                                                                                             tGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAgg        >  1:1798043/1‑150 (MQ=255)
                                                                                                                                                  gtgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCg   >  1:4324750/1‑150 (MQ=255)
                                                                                                                                                   tgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCgg  >  2:170100/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGG  >  NC_002947/4552764‑4553060

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: