Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,461,488 G→C 8.5% G199A (GGC→GCC)  paaJ → beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,461,4880GC 8.5% 325.9 / 18.4 106G199A (GGC→GCC) paaJbeta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (46/51);  new base C (3/6);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 109 total aligned reads displayed.

TGTGAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATT  >  NC_000913/1461339‑1461636
                                                                                                                                                     |                                                                                                                                                    
tgtgAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAgg                                                                                                                                                      >  1:1569457/1‑150 (MQ=255)
tgtgAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAgg                                                                                                                                                      <  2:1613036/150‑1 (MQ=255)
     aCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTc                                                                                                                                                 >  1:1050720/1‑150 (MQ=255)
        cGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTgg                                                                                                                                              >  1:414888/1‑150 (MQ=255)
         gCTCATGGCTCAGCAATTTGGAACTGACAGCCTGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGc                                                                                                                                             <  2:1239015/150‑1 (MQ=255)
         gCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAAGGGCATAGAATGTAGCTGAACTGTTATAAATCTCACGAGAAGATCAAGATAGTTTTTCGCTACGCAGTCAGGAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGc                                                                                                                                             <  2:844064/150‑1 (MQ=255)
         gCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGc                                                                                                                                             <  1:1264765/150‑1 (MQ=255)
          cTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCt                                                                                                                                            <  1:1192403/150‑1 (MQ=255)
           tCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTg                                                                                                                                           >  2:1207572/1‑150 (MQ=255)
            cATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTga                                                                                                                                          <  2:878210/150‑1 (MQ=255)
             aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                         <  1:638505/150‑1 (MQ=255)
              tGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgagg                                                                                                                                        <  1:1191544/150‑1 (MQ=255)
              tGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgagg                                                                                                                                        <  1:1257102/150‑1 (MQ=255)
               ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCGAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                       >  1:518154/1‑150 (MQ=255)
               ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgagga                                                                                                                                       <  1:1234386/150‑1 (MQ=255)
               ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgagga                                                                                                                                       <  1:31922/150‑1 (MQ=255)
                gCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggat                                                                                                                                      <  2:114704/150‑2 (MQ=255)
                     gCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGt                                                                                                                                 >  1:1501316/1‑150 (MQ=255)
                        aTTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTcc                                                                                                                              <  2:1050720/150‑1 (MQ=255)
                          ttGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTCCGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCgg                                                                                                                            >  1:1480428/1‑150 (MQ=255)
                             gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                         <  1:267994/150‑1 (MQ=255)
                             gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                         <  2:852126/150‑1 (MQ=255)
                               aCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtgtg                                                                                                                       <  2:1362643/150‑1 (MQ=255)
                                cTGACAGCATGCCTGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTgtgt                                                                                                                      <  1:711238/150‑1 (MQ=255)
                                  gACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTg                                                                                                                    >  2:958928/1‑150 (MQ=255)
                                      gCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGaaaa                                                                                                                >  1:496602/1‑150 (MQ=255)
                                         tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                             >  2:1626425/1‑150 (MQ=255)
                                         tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                             <  1:1207572/150‑1 (MQ=255)
                                          gCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACaa                                                                                                            >  2:451757/1‑150 (MQ=255)
                                           ccGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAg                                                                                                           <  1:643044/150‑1 (MQ=255)
                                                  cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                    >  1:1486578/1‑150 (MQ=255)
                                                  cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                    >  1:483877/1‑150 (MQ=255)
                                                     cAGAGATTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCTAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                 <  1:216830/150‑1 (MQ=255)
                                                     cAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                 <  1:258636/150‑1 (MQ=255)
                                                     cAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                 <  2:546934/150‑1 (MQ=255)
                                                         gaATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAAc                                                                                             >  2:740338/1‑150 (MQ=255)
                                                         gaATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTCTTGTAAc                                                                                             >  1:554154/1‑150 (MQ=255)
                                                          aaTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACa                                                                                            <  1:788046/150‑1 (MQ=255)
                                                           aTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAg                                                                                           <  2:1178231/150‑1 (MQ=255)
                                                            tGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGa                                                                                          >  1:773711/1‑150 (MQ=255)
                                                               aGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAAt                                                                                       <  1:277398/150‑1 (MQ=255)
                                                                  tGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATaca                                                                                    >  2:638424/1‑150 (MQ=255)
                                                                  tGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATaca                                                                                    >  1:1028714/1‑150 (MQ=255)
                                                                   gAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                             >  2:703567/1‑124 (MQ=255)
                                                                   gAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                             <  1:703567/124‑1 (MQ=255)
                                                                        gTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatga                                                                              >  1:195493/1‑150 (MQ=255)
                                                                        cttAAAAATATAACTAGAAGATCAACATAGTTTTAAGCAACGCACTCACCAACTTACCTCAAAAGAGCAATCCTCAGGAATTCTGACTGAGGAGATTGTTACGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatga                                                                              >  2:1015231/2‑150 (MQ=255)
                                                                         ttAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgat                                                                             >  1:18231/1‑150 (MQ=255)
                                                                          tAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgatg                                                                            <  1:1225403/150‑1 (MQ=255)
                                                                          tAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgatg                                                                            >  1:1616227/1‑150 (MQ=255)
                                                                          tAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGGGTTGTAACAGAAATACAACatgatg                                                                            <  2:1375959/150‑1 (MQ=255)
                                                                          tAAAAATCTCACGAGAACATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgatg                                                                            <  2:1058419/150‑1 (MQ=255)
                                                                           aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                           <  1:535034/150‑1 (MQ=255)
                                                                           aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                           >  2:441958/1‑150 (MQ=255)
                                                                            aaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAg                                                                          <  1:1034650/150‑1 (MQ=255)
                                                                            aaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAg                                                                          >  1:1172210/1‑150 (MQ=255)
                                                                                  tcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTg                                                                    >  1:1318084/1‑150 (MQ=255)
                                                                                       agaagaTCAAGATAGTTTCGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGccc                                                               <  2:1213205/150‑1 (MQ=255)
                                                                                        gaatttCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGAGCAATCCTCAGGCATTCTGGCTGATGAGATTGTTTCGGTTTTTTTGAATAACAAGAAAGGTTTTGTTACAGAAATAGATCATGATTATCATCTGCGcctg                                                              <  2:1081104/145‑3 (MQ=255)
                                                                                         aataTCAAGATAGTTTTGCGTTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAATACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCgg                                                             <  2:875326/147‑1 (MQ=255)
                                                                                         aagaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCgg                                                             <  2:1172210/150‑1 (MQ=255)
                                                                                           gaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGaa                                                           >  2:477812/1‑150 (MQ=255)
                                                                                           gaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGaa                                                           >  2:1152451/1‑150 (MQ=255)
                                                                                            aTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGaaa                                                          <  1:1221728/150‑1 (MQ=255)
                                                                                              cAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacg                                                        <  1:961169/150‑1 (MQ=255)
                                                                                              cAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacg                                                        >  2:1257215/1‑150 (MQ=255)
                                                                                                  aTAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAGACGACGc                                                    <  2:165537/150‑1 (MQ=255)
                                                                                                  aTAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGc                                                    <  1:1057195/150‑1 (MQ=255)
                                                                                                     gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                 <  2:539436/150‑1 (MQ=255)
                                                                                                     gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                 >  1:411748/1‑150 (MQ=255)
                                                                                                     gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                 <  2:18231/150‑1 (MQ=255)
                                                                                                      ttttGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGa                                                <  2:414888/150‑1 (MQ=255)
                                                                                                          gcgcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAg                                            <  2:1028714/150‑1 (MQ=255)
                                                                                                           cgcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGt                                           >  1:1137935/1‑150 (MQ=255)
                                                                                                            gcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGtt                                          >  1:1176088/1‑150 (MQ=255)
                                                                                                            gcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGtt                                          >  1:664770/1‑150 (MQ=255)
                                                                                                             cTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTa                                         <  2:1589812/150‑1 (MQ=255)
                                                                                                                     tCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGtt                                 <  2:528148/150‑1 (MQ=255)
                                                                                                                        gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                              >  1:739951/1‑150 (MQ=255)
                                                                                                                          aaCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAg                            <  1:1447538/150‑1 (MQ=255)
                                                                                                                            cGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGGAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCa                          >  2:898360/1‑150 (MQ=255)
                                                                                                                               aCGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCa                       >  1:772859/1‑150 (MQ=255)
                                                                                                                                cGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAt                      >  1:1620422/1‑150 (MQ=255)
                                                                                                                                        gcgcAATCCTCAGGCATTCTGTCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGcc              <  2:763988/150‑1 (MQ=255)
                                                                                                                                        gcgcAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGcc              <  2:411748/150‑1 (MQ=255)
                                                                                                                                        gcgcAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGcc              <  2:1537942/150‑1 (MQ=255)
                                                                                                                                           cAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt           >  2:1221107/1‑150 (MQ=255)
                                                                                                                                           atctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt           >  1:1628891/4‑150 (MQ=255)
                                                                                                                                            aaTCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATg          <  1:451757/150‑1 (MQ=255)
                                                                                                                                            aaTCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATg          <  2:739951/150‑1 (MQ=255)
                                                                                                                                             ctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                              >  2:1551994/2‑129 (MQ=255)
                                                                                                                                             ctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                              <  1:1551994/128‑1 (MQ=255)
                                                                                                                                              tCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATggg        <  2:772859/150‑1 (MQ=255)
                                                                                                                                              tCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATggg        >  1:997449/1‑150 (MQ=255)
                                                                                                                                                cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGt      <  2:354965/150‑1 (MQ=255)
                                                                                                                                                cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGt      <  1:1274370/150‑1 (MQ=255)
                                                                                                                                                 tCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTg     >  2:317038/1‑150 (MQ=255)
                                                                                                                                                  cAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGa    <  1:1152451/150‑1 (MQ=255)
                                                                                                                                                  cAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGa    >  1:1295355/1‑150 (MQ=255)
                                                                                                                                                    ggCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAtt  >  1:846363/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGTGAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATT  >  NC_000913/1461339‑1461636

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: