Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,873 C→A 11.1% T92T (ACC→ACA) ‡ 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,8730CA11.1% 621.6 / 54.0 209T92T (ACC→ACA) ‡pydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base C (88/97);  new base A (5/18);  total (94/115)
Fisher's exact test for biased strand distribution p-value = 2.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 216 total aligned reads displayed.

CTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCC  >  NC_002947/4552726‑4553022
                                                                                                                                                   |                                                                                                                                                     
cTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCa                                                                                                                                                     <  2:831006/150‑1 (MQ=255)
 tGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAg                                                                                                                                                    <  1:1884870/150‑1 (MQ=255)
    aaGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTg                                                                                                                                                 >  1:1463529/1‑150 (MQ=255)
      gACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaa                                                                                                                                               <  1:2564869/150‑1 (MQ=255)
        cccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaaga                                                                                                                                             <  2:3343701/150‑1 (MQ=255)
         ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                                                            <  1:3710833/150‑1 (MQ=255)
         ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                                                            <  2:1492536/150‑1 (MQ=255)
         ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                                                            >  2:2896559/1‑150 (MQ=255)
         ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaagaa                                                                                                                                            <  2:5354245/150‑1 (MQ=255)
          cTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGaagaag                                                                                                                                           <  1:4730511/150‑1 (MQ=255)
             ggCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGAc                                                                                                                                        <  1:982598/150‑1 (MQ=255)
                gAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTgg                                                                                                                                     >  2:4116750/1‑150 (MQ=255)
                gAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTgg                                                                                                                                     >  1:1358278/1‑150 (MQ=255)
                 aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGc                                                                                                                                    <  1:2901437/150‑1 (MQ=255)
                 aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGc                                                                                                                                    <  1:1360303/150‑1 (MQ=255)
                   gACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGccc                                                                                                                                  <  1:1748369/150‑1 (MQ=255)
                   gACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGccc                                                                                                                                  <  2:4605866/150‑1 (MQ=255)
                    aCCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCg                                                                                                                                 >  2:883567/1‑150 (MQ=255)
                     cccTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGa                                                                                                                                >  2:4204836/1‑150 (MQ=255)
                       cTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGAcc                                                                                                                              <  1:1483449/150‑1 (MQ=255)
                               gTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTg                                                                                                                      >  2:4845544/1‑150 (MQ=255)
                                  aaCGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATc                                                                                                                   >  1:1066941/1‑150 (MQ=255)
                                  aaCGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATc                                                                                                                   >  1:3566597/1‑150 (MQ=255)
                                         cgtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGc                                                                                                            <  2:3232166/150‑1 (MQ=255)
                                          gtcgCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCt                                                                                                           >  1:3142096/1‑150 (MQ=255)
                                            cgcgCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGa                                                                                                         <  1:3019374/150‑1 (MQ=255)
                                                cTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGt                                                                                                     <  2:5306133/150‑1 (MQ=255)
                                                   cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTg                                                                                                                                                 >  2:2604656/1‑103 (MQ=255)
                                                    tCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCg                                                                                                 >  2:4283227/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                                                                                        >  1:3726776/1‑150 (MQ=255)
                                                                       aCCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGa                                                                              <  2:4546539/150‑1 (MQ=255)
                                                                       aCCGTCAGGTACAGGGGATCAACAATATCCAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGa                                                                              <  2:1369742/150‑1 (MQ=255)
                                                                        ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                             <  1:1596841/150‑1 (MQ=255)
                                                                        ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                             <  2:4219018/150‑1 (MQ=255)
                                                                        ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                             <  2:1577851/150‑1 (MQ=255)
                                                                             aGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCa                                                                        <  1:3414651/150‑1 (MQ=255)
                                                                             aGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCa                                                                        <  1:4098502/150‑1 (MQ=255)
                                                                                  cAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTg                                                                   <  1:815811/150‑1 (MQ=255)
                                                                                     gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCg                                                                >  1:5039289/1‑150 (MQ=255)
                                                                                       gATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCt                                                              >  2:2809816/1‑150 (MQ=255)
                                                                                        aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                             <  2:4797145/150‑1 (MQ=255)
                                                                                        aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                             >  1:1787459/1‑150 (MQ=255)
                                                                                        aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                             <  2:4939755/150‑1 (MQ=255)
                                                                                        aTCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtg                                                             <  1:5400914/150‑1 (MQ=255)
                                                                                          caacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                                           >  2:5193066/1‑150 (MQ=255)
                                                                                          caacaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                                           <  1:350250/150‑1 (MQ=255)
                                                                                          caacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggt                                                           <  2:2477013/150‑1 (MQ=255)
                                                                                           aacaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtg                                                          <  2:2599013/150‑1 (MQ=255)
                                                                                            acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtgg                                                         >  1:1696462/1‑150 (MQ=255)
                                                                                            acaaTATCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCtggtgg                                                         >  1:3783897/1‑150 (MQ=255)
                                                                                              aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                       <  2:3726776/150‑1 (MQ=255)
                                                                                              aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                       <  2:4886366/150‑1 (MQ=255)
                                                                                              aaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGACCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGaa                                                       <  1:3375539/150‑1 (MQ=255)
                                                                                                  tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                   >  2:3473710/1‑150 (MQ=255)
                                                                                                  tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                   >  1:4087045/1‑150 (MQ=255)
                                                                                                  tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                                   >  1:2385388/1‑150 (MQ=255)
                                                                                                  tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                   <  1:713479/150‑4 (MQ=255)
                                                                                                  tCGAGCTGATCACCGACCGGCCGCTCGAAGTGAACCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAgaaa                                                   >  2:1688722/1‑147 (MQ=255)
                                                                                                   cGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCAc                                                  >  1:5385803/1‑150 (MQ=255)
                                                                                                          cccACCGACCGCTCACTGGAAACCACCCTGCGTGAAATCACCCCGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                           <  1:3822131/148‑1 (MQ=255)
                                                                                                           tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg                                          >  1:3682316/1‑150 (MQ=255)
                                                                                                             accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                        >  1:4367402/1‑150 (MQ=255)
                                                                                                              ccgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                                       >  2:3283472/1‑150 (MQ=255)
                                                                                                                 accgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCa                                    <  1:4708474/150‑1 (MQ=255)
                                                                                                                   cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                                  >  1:545254/1‑150 (MQ=255)
                                                                                                                   cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                                  >  1:705228/1‑150 (MQ=255)
                                                                                                                      tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAg                               <  1:2645160/150‑1 (MQ=255)
                                                                                                                        tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                                        >  1:1040066/1‑91 (MQ=255)
                                                                                                                        tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                             <  1:2809816/150‑1 (MQ=255)
                                                                                                                        tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                             >  1:2223093/1‑150 (MQ=255)
                                                                                                                         ctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTg                            >  2:2335581/1‑150 (MQ=255)
                                                                                                                          tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                             <  2:752346/100‑1 (MQ=255)
                                                                                                                          tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                             >  1:752346/1‑100 (MQ=255)
                                                                                                                           ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                          <  1:5110894/150‑1 (MQ=255)
                                                                                                                           ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                          <  1:2072282/150‑1 (MQ=255)
                                                                                                                           ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                          >  2:4239309/1‑150 (MQ=255)
                                                                                                                           ggAAATCAACCTGCGTGAAATCACCCAGGAGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                          <  2:4805776/150‑1 (MQ=255)
                                                                                                                             aaaTCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACt                        >  2:4753903/1‑150 (MQ=255)
                                                                                                                               gtgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACttt                      <  1:4548384/147‑2 (MQ=255)
                                                                                                                                tgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACtttg                     <  2:3783897/148‑3 (MQ=255)
                                                                                                                                tCAACCTGCGTGATATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCg                     <  2:490475/150‑1 (MQ=255)
                                                                                                                                 gaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTgg                    <  2:2222128/149‑1 (MQ=255)
                                                                                                                                  aaCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCg                                                                                  >  2:3695506/1‑87 (MQ=255)
                                                                                                                                   aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt                  <  1:2670924/150‑1 (MQ=255)
                                                                                                                                      tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc               >  2:295681/1‑150 (MQ=255)
                                                                                                                                       gCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGccc              >  2:3082282/1‑150 (MQ=255)
                                                                                                                                          tGAAATCACCCAGGTGAAGATGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa           >  1:4853057/1‑150 (MQ=255)
                                                                                                                                          tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa           >  2:3509034/1‑150 (MQ=255)
                                                                                                                                          tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa           >  2:1443782/1‑150 (MQ=255)
                                                                                                                                          tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa           >  1:5004626/1‑150 (MQ=255)
                                                                                                                                               tCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGAAGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTa      <  1:5212771/150‑1 (MQ=255)
                                                                                                                                                  cccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGc   >  2:2973288/1‑150 (MQ=255)
                                                                                                                                                   ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGaa                                                                             <  2:5038363/75‑1 (MQ=255)
                                                                                                                                                   ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggggg                                             <  2:4710457/107‑1 (MQ=255)
                                                                                                                                                   ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggggg                                             >  1:4710457/1‑107 (MQ=255)
                                                                                                                                                   ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGcc  <  2:1613108/150‑1 (MQ=255)
                                                                                                                                                   ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGcc  >  2:1085646/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCC  >  NC_002947/4552726‑4553022

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: