Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,844 T→C 9.3% S83P (TCT→CCT) ‡ 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,8440TC 9.3% 537.1 / 31.4 183S83P (TCT→CCT) ‡pydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base T (76/90);  new base C (5/12);  total (81/102)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 190 total aligned reads displayed.

CCTTCGAGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAG  >  NC_002947/4552698‑4552993
                                                                                                                                                  |                                                                                                                                                     
ccTTCGAGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGctctc                                                                                                                                                    <  2:1032714/150‑1 (MQ=255)
 cTTCGAGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGctctct                                                                                                                                                   <  2:4956122/150‑1 (MQ=255)
     gAGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGaa                                                                                                                                               >  1:1074405/1‑150 (MQ=255)
     gAGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGaa                                                                                                                                               <  2:4542385/150‑1 (MQ=255)
      aGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGata                                                                                                                                              <  1:3658128/150‑3 (MQ=255)
        gCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATc                                                                                                                                            >  1:3359241/1‑150 (MQ=255)
             cTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCt                                                                                                                                       <  1:2961151/150‑1 (MQ=255)
                 ggcggcGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGt                                                                                                                                   <  1:1178819/150‑1 (MQ=255)
                 ggcggcGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGt                                                                                                                                   <  1:751071/150‑1 (MQ=255)
                            cTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCa                                                                                                                        <  2:831006/150‑1 (MQ=255)
                            cTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCa                                                                                                                        <  2:3723894/150‑1 (MQ=255)
                             tGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAg                                                                                                                       <  1:1884870/150‑1 (MQ=255)
                               gAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGt                                                                                                                     >  1:3579270/1‑150 (MQ=255)
                                  gACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaa                                                                                                                  <  1:2564869/150‑1 (MQ=255)
                                     ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                               <  2:5354265/150‑1 (MQ=255)
                                     ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                               <  1:3710833/150‑1 (MQ=255)
                                     ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                               >  2:2896559/1‑150 (MQ=255)
                                         ggCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGAc                                                                                                           <  1:982598/150‑1 (MQ=255)
                                            gAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAAc                                                                                                                                         <  1:5215974/117‑1 (MQ=255)
                                            gAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTgg                                                                                                        >  1:4062715/1‑150 (MQ=255)
                                            gAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTgg                                                                                                        >  2:4116750/1‑150 (MQ=255)
                                             aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGc                                                                                                       <  1:2901437/150‑1 (MQ=255)
                                                    tGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGTGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACcg                                                                                                >  2:5212771/1‑150 (MQ=255)
                                                           gTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTg                                                                                         >  2:2108947/1‑150 (MQ=255)
                                                           gTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTg                                                                                         >  2:4845544/1‑150 (MQ=255)
                                                              aaCGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATc                                                                                      <  2:4310323/150‑1 (MQ=255)
                                                              aaCGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATc                                                                                      >  2:4826681/1‑150 (MQ=255)
                                                                  tgtCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCgtgt                                                                                  <  2:1732166/150‑1 (MQ=255)
                                                                    tcgtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCg                                                                                <  2:1431876/150‑1 (MQ=255)
                                                                     cgtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGc                                                                               <  2:3232166/150‑1 (MQ=255)
                                                                      gtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGaa                                                                                                                                >  2:3157076/1‑100 (MQ=255)
                                                                      gtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGaa                                                                                                                                <  1:3157076/100‑1 (MQ=255)
                                                                       tcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGc                                                                                                                                     >  1:1013875/1‑94 (MQ=255)
                                                                               cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTg                                                                                                                    <  1:2604656/103‑1 (MQ=255)
                                                                               cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGcc                                                                     >  1:381183/1‑150 (MQ=255)
                                                                               cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGcc                                                                     <  1:485005/150‑1 (MQ=255)
                                                                                 cGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGt                                                                   <  1:2082465/150‑1 (MQ=255)
                                                                                  ggCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTg                                                                  >  1:83884/1‑150 (MQ=255)
                                                                                        cTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATaga                                                            <  2:3945429/150‑1 (MQ=255)
                                                                                         tATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                           >  1:4862212/1‑150 (MQ=255)
                                                                                         tATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                           <  2:5256139/150‑1 (MQ=255)
                                                                                              ccccAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATc                                                      <  1:1774724/150‑1 (MQ=255)
                                                                                                  aaCCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTgg                                                  <  1:4088234/150‑1 (MQ=255)
                                                                                                    ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                <  1:1596841/150‑1 (MQ=255)
                                                                                                    ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                <  2:4219018/150‑1 (MQ=255)
                                                                                                    ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                <  1:1463946/150‑1 (MQ=255)
                                                                                                    ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                <  2:4218924/150‑1 (MQ=255)
                                                                                                    ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGCACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                <  2:3067076/150‑1 (MQ=255)
                                                                                                     cGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAg                                               >  1:4927021/1‑150 (MQ=255)
                                                                                                        cAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTc                                            <  2:5379538/150‑1 (MQ=255)
                                                                                                             aCAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCt                                       >  1:350278/1‑150 (MQ=255)
                                                                                                                tggtctCAACAATATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGcc                                    <  2:3780113/145‑1 (MQ=255)
                                                                                                                 gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCg                                   >  1:338798/1‑150 (MQ=255)
                                                                                                                    aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                <  2:4797145/150‑1 (MQ=255)
                                                                                                                    aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                <  2:4939755/150‑1 (MQ=255)
                                                                                                                    aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                <  1:223006/150‑1 (MQ=255)
                                                                                                                    aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                <  1:5400961/150‑1 (MQ=255)
                                                                                                                    aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                <  2:5314896/150‑1 (MQ=255)
                                                                                                                    aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                >  1:1787459/1‑150 (MQ=255)
                                                                                                                     tCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtgg                               <  1:1463068/150‑1 (MQ=255)
                                                                                                                      caacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                              >  2:2601202/1‑150 (MQ=255)
                                                                                                                      caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                              <  1:372316/150‑1 (MQ=255)
                                                                                                                      caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                              <  1:4538041/150‑1 (MQ=255)
                                                                                                                      caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                              <  2:2477013/150‑1 (MQ=255)
                                                                                                                      caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                              <  1:3101817/150‑1 (MQ=255)
                                                                                                                       aacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGTAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtg                             >  2:1593466/1‑150 (MQ=255)
                                                                                                                        acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtgg                            >  1:1696462/1‑150 (MQ=255)
                                                                                                                        acaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtgg                            >  1:3783834/1‑150 (MQ=255)
                                                                                                                         caaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGa                           >  1:4253096/1‑150 (MQ=255)
                                                                                                                          aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                          <  2:4886366/150‑1 (MQ=255)
                                                                                                                          aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                          <  2:3345090/150‑1 (MQ=255)
                                                                                                                          aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                          <  1:791928/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)
                                                                                                                              tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                      >  2:3473710/1‑150 (MQ=255)
                                                                                                                              tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                      >  2:1195446/1‑150 (MQ=255)
                                                                                                                              tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                      <  1:713479/150‑4 (MQ=255)
                                                                                                                               cGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCAc                     >  1:5385803/1‑150 (MQ=255)
                                                                                                                                gAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACg                    >  2:4597525/1‑150 (MQ=255)
                                                                                                                                 aGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACgg                   <  1:1359411/150‑1 (MQ=255)
                                                                                                                                   ttGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACgggg                 <  2:2583942/150‑1 (MQ=255)
                                                                                                                                    tGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggggg                <  2:273320/150‑1 (MQ=255)
                                                                                                                                     gATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGc               >  2:28724/1‑150 (MQ=255)
                                                                                                                                      aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc              >  2:1653659/1‑150 (MQ=255)
                                                                                                                                      aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc              <  2:5129264/150‑1 (MQ=255)
                                                                                                                                      aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc              <  2:3315803/150‑1 (MQ=255)
                                                                                                                                       tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg             >  1:2979342/1‑150 (MQ=255)
                                                                                                                                         accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTc                                            >  2:4713922/1‑117 (MQ=255)
                                                                                                                                         accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTc                                            <  1:4713922/117‑1 (MQ=255)
                                                                                                                                         accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt           <  1:3415609/150‑1 (MQ=255)
                                                                                                                                          ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg          <  2:2304102/150‑1 (MQ=255)
                                                                                                                                          ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg          >  2:3283472/1‑150 (MQ=255)
                                                                                                                                           cgaccgGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCgaagaagaagc                                                      >  1:1571847/1‑103 (MQ=255)
                                                                                                                                           cgaccgGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCgaagaagaagc                                                      <  2:1571847/105‑3 (MQ=255)
                                                                                                                                             accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa       <  1:5055433/150‑1 (MQ=255)
                                                                                                                                             accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa       <  1:4708474/150‑1 (MQ=255)
                                                                                                                                               cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc     >  1:705228/1‑150 (MQ=255)
                                                                                                                                               cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc     >  1:545254/1‑150 (MQ=255)
                                                                                                                                                  tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAg  <  1:2645188/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCTTCGAGGCCGGCTGGGGCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAG  >  NC_002947/4552698‑4552993

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: