Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,976 A→T 8.0% G160G (GGA→GGT)  paaF → 2,3‑dehydroadipyl‑CoA hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,457,9760AT 8.0% 261.3 / 10.6 87G160G (GGA→GGT) paaF2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base T (5/2);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCT  >  NC_000913/1457828‑1458122
                                                                                                                                                    |                                                                                                                                                   
caatgttgtGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCACGCGGAACGCAACGTTTAATCCCTAGTGTCGGAAAATCGTTAGCCAGAAAAATGGTGCTGAGCGga                                                                                                                                                     >  2:314083/4‑149 (MQ=255)
cATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGga                                                                                                                                                     >  1:1582318/1‑149 (MQ=255)
cATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGga                                                                                                                                                     >  2:782739/1‑149 (MQ=255)
  ttgttgTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGgaga                                                                                                                                                   >  2:1579373/1‑149 (MQ=255)
    gttgtGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGaaa                                                                                                                                                 <  2:378210/149‑1 (MQ=255)
        tgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTAt                                                                                                                                             >  1:943144/1‑149 (MQ=255)
                 gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCa                                                                                                                                    <  2:1258932/149‑1 (MQ=255)
                 gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                    <  1:515902/149‑1 (MQ=255)
                 gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                    <  1:557617/149‑1 (MQ=255)
                 gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                    <  1:1349414/149‑1 (MQ=255)
                 gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                    >  2:1085268/1‑149 (MQ=255)
                  tggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCag                                                                                                                                   <  2:1395149/149‑1 (MQ=255)
                     gtGCCGGAGAGAACGCGCGTGTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTAGCACCGCTCagcaa                                                                                                                                <  1:1501636/148‑1 (MQ=255)
                            agagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagc                                                                                                                         <  2:900520/149‑1 (MQ=255)
                                   gcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                  >  1:1438804/1‑149 (MQ=255)
                                   gcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                  >  1:264632/1‑149 (MQ=255)
                                    cgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaa                                                                                                                                >  1:972929/1‑134 (MQ=255)
                                    agcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaa                                                                                                                                <  2:972929/133‑1 (MQ=255)
                                       gTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTg                                                                                                              <  2:699938/149‑1 (MQ=255)
                                             ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                        >  2:999255/1‑149 (MQ=255)
                                             ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                        >  2:478438/1‑149 (MQ=255)
                                             ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                        >  2:878435/1‑149 (MQ=255)
                                              gTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCg                                                                                                       >  2:1162130/1‑149 (MQ=255)
                                                tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGAc                                                                                                     <  1:1388959/149‑1 (MQ=255)
                                                tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGAc                                                                                                     <  1:782739/149‑1 (MQ=255)
                                                ggCCGGAAATCACTCTCGTCATCACGCCCGGCGCAGGCGGAACGCATCGTTTAATGCGTCGTGTCGGTAAATCGTTAGCCAGCAAACTGGTGCTGAGCGGAGAATGTTTCAACGCTCAGCAAGCTGAGCAGGCCGGGCTGTTTAGCGAc                                                                                                     <  1:314083/148‑1 (MQ=255)
                                                  ccGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGt                                                                                                   >  2:590691/1‑149 (MQ=255)
                                                      aactCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGGTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAACTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑c                                                                                              <  1:129748/146‑1 (MQ=255)
                                                      aaaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑c                                                                                              <  2:1566490/149‑1 (MQ=255)
                                                      aaaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑c                                                                                              <  1:1397353/149‑1 (MQ=255)
                                                        atctgtctcTGCATCATGCGAGGCGCAGGCGGAACGCAACGCTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAATTGGTGCTGAGCGGATATAGTATCACCGCTGAGCAAGCACAGCAGTCATGGCTGGTTATCGACGTCTT‑ccc                                                                                            <  2:1280392/144‑1 (MQ=255)
                                                         tCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑cccc                                                                                           <  1:1523710/149‑1 (MQ=255)
                                                           aCTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAg                                                                                         >  2:230111/1‑149 (MQ=255)
                                                            ctctCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGc                                                                                        >  2:1150194/1‑149 (MQ=255)
                                                               tcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGAt                                                                                     <  1:675922/149‑1 (MQ=255)
                                                                   catcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTaa                                                                                 >  1:1568659/1‑149 (MQ=255)
                                                                       atGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCt                                                                             >  1:669724/1‑149 (MQ=255)
                                                                        tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTc                                                                            <  1:544562/149‑1 (MQ=255)
                                                                        tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTc                                                                            <  1:1115770/149‑1 (MQ=255)
                                                                        tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGGTAGCGACGTCTT‑CCCCAGCGATTTAACCCTc                                                                            <  2:261057/149‑1 (MQ=255)
                                                                        tGCCTGGCGCAGGCGGAACGAAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTc                                                                            <  1:1388498/149‑1 (MQ=255)
                                                                        tGCCTGGCGCAAGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTc                                                                            <  1:1430702/149‑1 (MQ=255)
                                                                          ccTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCAT‑CCCCAGCGATTTAACCCTCGa                                                                          >  2:1098773/1‑149 (MQ=255)
                                                                          ccTGGCGCAGGCGGAACGCAAAGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGGACAGCAGGACGGGCTGGTGAGCGACGTCTG‑CCCCAGCGCTTTAACCCTcta                                                                          >  2:852501/1‑147 (MQ=255)
                                                                           cTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGaa                                                                         <  1:1508900/149‑1 (MQ=255)
                                                                            tGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAAt                                                                        <  1:1579373/149‑1 (MQ=255)
                                                                             ggCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCGGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATa                                                                       >  2:543381/1‑149 (MQ=255)
                                                                             ggCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCt                                                                             <  2:80120/143‑1 (MQ=255)
                                                                             ggCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCt                                                                             >  1:80120/1‑143 (MQ=255)
                                                                               cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACg                                                                     <  1:878435/149‑1 (MQ=255)
                                                                               cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACg                                                                     <  1:999255/149‑1 (MQ=255)
                                                                               cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACg                                                                     <  1:478438/149‑1 (MQ=255)
                                                                                gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGc                                                                    <  1:145860/149‑1 (MQ=255)
                                                                                    gCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCTGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTa                                                                <  1:543381/149‑1 (MQ=255)
                                                                                    gCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGAAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTa                                                                >  1:90545/1‑149 (MQ=255)
                                                                                     cGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCCAGCGATTTAACCCTCGAATACGCCTTa                                                                <  2:943144/149‑1 (MQ=255)
                                                                                      ggAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACa                                                              >  1:1571709/1‑149 (MQ=255)
                                                                                      ggAACGCAACGTTTAATACGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTT‑CACCAGCGCTTTGACCCTCGCATACGCCTTACa                                                              >  1:791710/1‑149 (MQ=255)
                                                                                           gCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTgg                                                         <  2:1334824/149‑1 (MQ=255)
                                                                                             aaCGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                       >  2:1205981/1‑149 (MQ=255)
                                                                                             aaCGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                       >  1:978232/1‑149 (MQ=255)
                                                                                              aCGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCAt                                                      >  2:1086253/1‑149 (MQ=255)
                                                                                               cGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATc                                                     <  2:1438804/149‑1 (MQ=255)
                                                                                               cGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATc                                                     <  2:264632/149‑1 (MQ=255)
                                                                                                     aTCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAAt                                               <  2:592471/149‑1 (MQ=255)
                                                                                                      tCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATg                                              >  1:164107/1‑149 (MQ=255)
                                                                                                       ccGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATgg                                             <  2:350301/149‑1 (MQ=255)
                                                                                                        cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGc                                            <  1:230111/149‑1 (MQ=255)
                                                                                                         gTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                           >  1:1240962/1‑149 (MQ=255)
                                                                                                         gTAGTGTAGGTAACTCGTTAGCCAGCAAAATGGTGCGGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                           <  1:1436771/149‑1 (MQ=255)
                                                                                                               tCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTACCCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCAc                                     >  2:415495/1‑149 (MQ=255)
                                                                                                               tCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCAc                                     >  1:1492112/1‑149 (MQ=255)
                                                                                                                 ggTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTc                                   >  1:1509117/1‑149 (MQ=255)
                                                                                                                  gTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcg                                  >  2:1511619/1‑149 (MQ=255)
                                                                                                                       tCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTGACAGCTGGCATCGAAAATGGCACGTCACTcgccgcc                             <  1:1063773/149‑2 (MQ=255)
                                                                                                                          ttAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGc                          <  2:835158/149‑1 (MQ=255)
                                                                                                                             gCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCtt                       >  2:535648/1‑149 (MQ=255)
                                                                                                                              ccAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTT‑CCCCAGCGATTTACCCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                      >  1:63820/1‑149 (MQ=255)
                                                                                                                              ccAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                      >  1:818632/1‑149 (MQ=255)
                                                                                                                               cAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTAc                     >  2:193164/1‑149 (MQ=255)
                                                                                                                                     aaTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCgg               <  1:1162130/149‑1 (MQ=255)
                                                                                                                                         ttGCTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCaaa           >  1:1145768/2‑149 (MQ=255)
                                                                                                                                         gttcTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCaaa           <  2:1456470/146‑1 (MQ=255)
                                                                                                                                            cTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagca        >  1:921569/1‑149 (MQ=255)
                                                                                                                                            cTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagca        >  2:1308680/1‑149 (MQ=255)
                                                                                                                                             tGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa       >  1:641789/1‑149 (MQ=255)
                                                                                                                                             tGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa       >  2:948306/1‑149 (MQ=255)
                                                                                                                                             tGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa       >  2:1408675/1‑149 (MQ=255)
                                                                                                                                             tGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa       <  1:535648/149‑1 (MQ=255)
                                                                                                                                               aGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaagc     >  2:1531263/1‑149 (MQ=255)
                                                                                                                                                  ggAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCt  >  2:96533/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTT‑CCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCT  >  NC_000913/1457828‑1458122

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