Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,461,493 C→A 5.8% L201M (CTG→ATG)  paaJ → beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,461,4930CA 5.8% 1219.4 / 34.5 357L201M (CTG→ATG) paaJbeta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (165/169);  new base A (16/6);  total (181/176)
Fisher's exact test for biased strand distribution p-value = 4.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 360 total aligned reads displayed.

ACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACCGCA  >  NC_000913/1461344‑1461642
                                                                                                                                                     |                                                                                                                                                     
   cGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTgg                                                                                                                                                    >  2:2108256/1‑150 (MQ=255)
    gCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGc                                                                                                                                                   <  1:595727/150‑1 (MQ=255)
     cTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCt                                                                                                                                                  >  2:396743/1‑150 (MQ=255)
      tCATGGCTCAGCAATTTGGAACTGACAGCATGCCTGAAAGGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTg                                                                                                                                                 <  1:2448330/150‑1 (MQ=255)
       cATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAACAGCGCAATCCTCAGGCATTCGGGcggg                                                                                                                                                >  2:1808776/1‑147 (MQ=255)
        aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                               <  2:1612976/150‑1 (MQ=255)
        aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                               <  2:1693051/150‑1 (MQ=255)
         tGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagg                                                                                                                                              <  2:1988443/150‑1 (MQ=255)
          ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                             >  2:1958292/1‑150 (MQ=255)
          ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                             >  2:178544/1‑150 (MQ=255)
           gCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgaggag                                                                                                                                            >  1:1792374/1‑150 (MQ=255)
            cTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGgaga                                                                                                                                           >  2:954427/1‑150 (MQ=255)
            cTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGgaga                                                                                                                                           >  2:902586/1‑150 (MQ=255)
             tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGAt                                                                                                                                          <  1:2092948/150‑1 (MQ=255)
                gCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGt                                                                                                                                       >  1:645678/1‑150 (MQ=255)
                   aTTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTcc                                                                                                                                    <  1:305034/150‑1 (MQ=255)
                      tGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                                 >  1:938291/1‑150 (MQ=255)
                      tGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                                 <  1:1170628/150‑1 (MQ=255)
                        gcaCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAACATCAAGATAGTTTTGCGCTACGCAGTCAGCACCGTACGGCAAACGCGCAATCCTCAGGCATTCTGGCTGCGGAGATTGTTCCGGTtg                                                                                                                               >  2:385711/3‑150 (MQ=255)
                        gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                               <  2:1068552/150‑1 (MQ=255)
                         aaCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGAGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtgt                                                                                                                              >  1:543294/1‑150 (MQ=255)
                          aCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCACGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtgtg                                                                                                                             <  2:1435820/150‑1 (MQ=255)
                          aCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtgtg                                                                                                                             >  2:1506093/1‑150 (MQ=255)
                          aCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtgtg                                                                                                                             <  1:183746/150‑1 (MQ=255)
                            tGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGtt                                                                                                                           >  2:1881852/1‑150 (MQ=255)
                                 gCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGaaaa                                                                                                                      <  2:33696/150‑1 (MQ=255)
                                  cATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGaaaaa                                                                                                                     >  1:1342019/1‑150 (MQ=255)
                                   aTGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAAc                                                                                                                    >  2:1657770/1‑150 (MQ=255)
                                   aTGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAAc                                                                                                                    >  2:1543519/1‑150 (MQ=255)
                                    tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                                   <  2:191481/150‑1 (MQ=255)
                                     gCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACaa                                                                                                                  <  2:1154551/150‑1 (MQ=255)
                                      ccGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAg                                                                                                                 <  1:1014187/150‑1 (MQ=255)
                                      ccGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGACATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAg                                                                                                                 >  2:449306/1‑150 (MQ=255)
                                            aCGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGTGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgt                                                                                                           <  2:2385989/150‑1 (MQ=255)
                                             cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                          >  2:572328/1‑150 (MQ=255)
                                             cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                          >  2:1280210/1‑150 (MQ=255)
                                                 agagAATGTAGCTGAACTGCTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttgt                                                                                                      <  2:2224380/150‑1 (MQ=255)
                                                      aTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAg                                                                                                 >  2:255876/1‑150 (MQ=255)
                                                         tAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGaaa                                                                                              >  2:1846008/1‑150 (MQ=255)
                                                             tGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATaca                                                                                          <  1:1524038/150‑1 (MQ=255)
                                                               aaattttAAAAATCTCACGAGAAGATCAAGATAGATTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGATTAGATTGTTCCGGTTGTGTTGAAAAACAACAAAGGTGTTTTAACAGAAATacaac                                                                                        <  1:1585165/145‑1 (MQ=255)
                                                               aaCTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATacaac                                                                                        <  1:1124001/150‑1 (MQ=255)
                                                                aCTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATacaaca                                                                                       <  1:35622/150‑1 (MQ=255)
                                                                   gTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatga                                                                                    >  2:671743/1‑150 (MQ=255)
                                                                   gTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatga                                                                                    <  2:1475380/150‑1 (MQ=255)
                                                                    ttAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgat                                                                                   >  1:2152916/1‑150 (MQ=255)
                                                                    ttAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgat                                                                                   >  1:2217350/1‑150 (MQ=255)
                                                                      aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                                 >  2:2418455/1‑150 (MQ=255)
                                                                      aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                                 >  2:1114517/1‑150 (MQ=255)
                                                                      aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                                 <  2:510904/150‑1 (MQ=255)
                                                                       aaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAg                                                                                <  1:2425370/150‑1 (MQ=255)
                                                                         aaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCa                                                                              <  2:2454349/150‑1 (MQ=255)
                                                                           tctcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATc                                                                            <  1:1562677/150‑1 (MQ=255)
                                                                           tctcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATc                                                                            <  1:2029980/150‑1 (MQ=255)
                                                                             tcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTg                                                                          >  2:272965/1‑150 (MQ=255)
                                                                                   gaagaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGATAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCCTCTGCGCCCg                                                                    >  2:596652/1‑150 (MQ=255)
                                                                                   gaagaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCg                                                                    >  1:814887/1‑150 (MQ=255)
                                                                                        tCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAac                                                               >  1:1694182/1‑150 (MQ=255)
                                                                                         cAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacg                                                              >  1:705554/1‑150 (MQ=255)
                                                                                            gATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacgacg                                                           >  1:7674/1‑150 (MQ=255)
                                                                                            gATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacgacg                                                           >  1:1902412/1‑150 (MQ=255)
                                                                                              tAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCt                                                         <  2:709584/150‑1 (MQ=255)
                                                                                                gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                       >  1:490022/1‑150 (MQ=255)
                                                                                                 ttttGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGa                                                      >  2:820226/1‑150 (MQ=255)
                                                                                                     gcgcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAg                                                  <  1:966979/150‑1 (MQ=255)
                                                                                                      cgcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGt                                                 <  2:894087/150‑1 (MQ=255)
                                                                                                        cTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTa                                               >  2:773325/1‑150 (MQ=255)
                                                                                                        cTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTa                                               >  1:2293486/1‑150 (MQ=255)
                                                                                                          aCGCAGTCAGCAACGTACGGCAAACGCGCAATCCTCAGGCAGTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGCCGCTGGCACAGTTACg                                             >  2:1206669/1‑150 (MQ=255)
                                                                                                          aCGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACg                                             <  1:999820/150‑1 (MQ=255)
                                                                                                          aCGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACg                                             <  1:1958292/150‑1 (MQ=255)
                                                                                                          aCGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACg                                             <  1:1431130/150‑1 (MQ=255)
                                                                                                            gCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTg                                           <  1:124328/150‑1 (MQ=255)
                                                                                                             cAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTgg                                          <  1:954427/150‑1 (MQ=255)
                                                                                                              aGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTggg                                         <  2:1643205/150‑1 (MQ=255)
                                                                                                              aGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTggg                                         <  2:1342019/150‑1 (MQ=255)
                                                                                                              aGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTggg                                         <  2:2133768/150‑1 (MQ=255)
                                                                                                              aGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTggg                                         >  2:865335/1‑150 (MQ=255)
                                                                                                                 cAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGa                                                                                                                         >  1:1078924/1‑67 (MQ=255)
                                                                                                                 cAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGa                                                                                                                         <  2:1078924/67‑1 (MQ=255)
                                                                                                                 cAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTa                                      >  1:1718983/1‑150 (MQ=255)
                                                                                                                   gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAGGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                                    >  1:835548/1‑150 (MQ=255)
                                                                                                                   gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                                    >  1:2063539/1‑150 (MQ=255)
                                                                                                                   gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                                    >  1:326266/1‑150 (MQ=255)
                                                                                                                      aCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAGGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGc                                 <  1:902586/150‑1 (MQ=255)
                                                                                                                      aCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGc                                 <  2:1757051/150‑1 (MQ=255)
                                                                                                                             gCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGc                                 <  2:1253490/143‑1 (MQ=255)
                                                                                                                              cAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTc                         >  1:1916608/1‑150 (MQ=255)
                                                                                                                               aaaaGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCg                        <  1:2395774/150‑1 (MQ=255)
                                                                                                                                aaaGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGt                       >  1:1418810/1‑150 (MQ=255)
                                                                                                                                 aaGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTg                      <  1:1655192/150‑1 (MQ=255)
                                                                                                                                   gcgcAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGcc                    >  2:428214/1‑150 (MQ=255)
                                                                                                                                     gcAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCaa                  >  2:2288761/1‑150 (MQ=255)
                                                                                                                                      atctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt                 >  2:2045711/4‑150 (MQ=255)
                                                                                                                                        ctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATgg               >  1:2183849/2‑150 (MQ=255)
                                                                                                                                         tCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATggg              >  1:1299654/1‑150 (MQ=255)
                                                                                                                                            tCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTg           <  1:2331005/150‑1 (MQ=255)
                                                                                                                                             cAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGa          <  1:2108256/150‑1 (MQ=255)
                                                                                                                                               ggCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGAGGTGAtt        >  2:1050684/1‑150 (MQ=255)
                                                                                                                                                     cTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACCGCa  <  1:396743/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACCGCA  >  NC_000913/1461344‑1461642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: