Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,534 T→C 8.7% L13S (TTG→TCG)  paaF → putative 2,3‑dehydroadipyl‑CoA hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,457,5340TC 8.7% 196.3 / 8.1 69L13S (TTG→TCG) paaFputative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/4);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTGGGG  >  NC_000913/1457389‑1457683
                                                                                                                                                 |                                                                                                                                                     
cAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAttgttg                                                                                                                                                   >  2:224433/1‑150 (MQ=255)
   ttACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTg                                                                                                                                                <  1:224433/150‑1 (MQ=255)
     aCAGTCTGGAACCGGATGAACTGGCCGCAGTTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTTACTTTGACGCGATGGGGATGGCATGAGCGAACTTATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                              <  2:152667/150‑1 (MQ=255)
     aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                              <  2:342119/150‑1 (MQ=255)
     aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                              <  1:1303482/150‑1 (MQ=255)
          cTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTa                                                                                                                                         >  2:1546778/1‑150 (MQ=255)
           tGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTaa                                                                                                                                        >  1:26164/1‑150 (MQ=255)
            ggAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAAc                                                                                                                                       >  2:312537/1‑150 (MQ=255)
            ggAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACAAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAAc                                                                                                                                       >  2:332456/1‑150 (MQ=255)
                 cGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCGCTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAAGGAGTATTGGTGGTGCCCGTTAACCGTcc                                                                                                                                  >  2:1458153/1‑150 (MQ=255)
                  ggATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTccc                                                                                                                                 <  2:670206/150‑1 (MQ=255)
                    aTGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCccgc                                                                                                                               >  1:1369705/1‑150 (MQ=255)
                     tGAAGTGGCCGCAGGTTATGTGTTGATTTGCGAGGAACTGCCGATGACCAGCGATGTGGTGGTTGACTTTTACGCGAAGGGGATGGGATGAGAGAACTGATAGTCAGCCGTCAGGTTCGAGTTTTGTTGCTGACCCTTAACCGTGgcgcc                                                                                                                              <  2:353473/150‑1 (MQ=255)
                      gAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCccgccg                                                                                                                             >  1:1354615/1‑150 (MQ=255)
                        aCTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCa                                                                                                                           <  1:1428220/150‑1 (MQ=255)
                         cTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCAc                                                                                                                          <  1:1159202/150‑1 (MQ=255)
                          tGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCAAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACg                                                                                                                         >  1:671194/1‑150 (MQ=255)
                               gCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATg                                                                                                                    >  1:931738/1‑150 (MQ=255)
                                cAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATgc                                                                                                                   <  1:530774/150‑1 (MQ=255)
                                  ggTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATgcgc                                                                                                                 <  2:827762/150‑1 (MQ=255)
                                   gTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCt                                                                                                                <  1:1090447/150‑1 (MQ=255)
                                   gTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCt                                                                                                                <  2:1486955/150‑1 (MQ=255)
                                     tATGTGTTGAGTTGCAAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTaa                                                                                                              >  2:1465347/1‑150 (MQ=255)
                                       tgtgTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaat                                                                                                            <  1:1089958/150‑1 (MQ=255)
                                         tgtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaataa                                                                                                          >  1:1385968/1‑150 (MQ=255)
                                          gtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaataat                                                                                                         <  2:1406128/150‑1 (MQ=255)
                                          gtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaataat                                                                                                         >  1:1399642/1‑150 (MQ=255)
                                           ttGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATg                                                                                                        <  1:1615443/150‑1 (MQ=255)
                                            tGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGc                                                                                                       <  1:1041584/150‑1 (MQ=255)
                                             gAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGcc                                                                                                      <  2:1434111/150‑1 (MQ=255)
                                                 tGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCctgc                                                                                                  <  1:1062217/150‑1 (MQ=255)
                                                    cAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGa                                                                                               >  2:481833/1‑150 (MQ=255)
                                                             ccGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGAGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGAGCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTgg                                                                                      <  2:1574259/150‑1 (MQ=255)
                                                               gCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTa                                                                                    >  2:1443428/1‑150 (MQ=255)
                                                               gCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTa                                                                                    >  2:1111117/1‑150 (MQ=255)
                                                                  gACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAAt                                                                                 >  2:275441/1‑150 (MQ=255)
                                                                    ccAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGa                                                                               >  2:889709/1‑150 (MQ=255)
                                                                      aGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAAc                                                                             <  1:889709/150‑1 (MQ=255)
                                                                        cGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCTGCTGATGCAACTGGTAAATGAACTgg                                                                          <  2:1369705/150‑1 (MQ=255)
                                                                          aTGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATg                                                                                <  1:1619069/143‑1 (MQ=255)
                                                                          aTGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATg                                                                                >  2:1619069/1‑143 (MQ=255)
                                                                           tgtgGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGaa                                                                        <  1:312537/150‑1 (MQ=255)
                                                                                 ggtTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCg                                                                  <  2:1555086/150‑1 (MQ=255)
                                                                                       cTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTAcc                                                            <  1:1546778/150‑1 (MQ=255)
                                                                                             cgcgAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATAcc                                                      <  2:931738/150‑1 (MQ=255)
                                                                                                 aaGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCa                                                  <  2:571837/150‑1 (MQ=255)
                                                                                                     ggATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTc                                              >  1:826149/1‑150 (MQ=255)
                                                                                                     ggATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGATGCGGCTACCGATACCAGCATTTc                                              >  2:1404888/1‑150 (MQ=255)
                                                                                                       aTGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCgg                                            >  1:495339/1‑150 (MQ=255)
                                                                                                        tGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGt                                           >  2:5463/1‑150 (MQ=255)
                                                                                                         ggCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTc                                          <  2:1399642/150‑1 (MQ=255)
                                                                                                            aTGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgt                                       <  1:332456/150‑1 (MQ=255)
                                                                                                               aGCGAACTGATCGTCAGCCGTCAGCACCGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgtgtg                                    <  2:1542748/150‑1 (MQ=255)
                                                                                                               aGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgtgtg                                    >  1:218180/1‑150 (MQ=255)
                                                                                                               aGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgtgtg                                    >  1:1188355/1‑150 (MQ=255)
                                                                                                                        aTCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGt                           <  1:1443428/150‑1 (MQ=255)
                                                                                                                           gtcagccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGGTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAAt                        >  1:162356/1‑150 (MQ=255)
                                                                                                                             taagggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTc                                              <  1:697152/121‑1 (MQ=255)
                                                                                                                             taagggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTc                                              >  2:697152/6‑126 (MQ=255)
                                                                                                                             cagccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGc                      <  2:518269/150‑1 (MQ=255)
                                                                                                                               gccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCAc                    >  2:323637/1‑150 (MQ=255)
                                                                                                                               gccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCAc                    <  2:671194/150‑1 (MQ=255)
                                                                                                                               agggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCAc                    >  1:778894/4‑150 (MQ=255)
                                                                                                                                    cagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCttt               >  1:986806/1‑150 (MQ=255)
                                                                                                                                      gcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGGAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCttttt             >  1:108421/1‑150 (MQ=255)
                                                                                                                                      gcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCttttt             >  1:1463188/1‑150 (MQ=255)
                                                                                                                                      gcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCttttt             <  2:495339/150‑1 (MQ=255)
                                                                                                                                        aaCGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTg           <  1:1111117/150‑1 (MQ=255)
                                                                                                                                           gAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCg        <  2:26164/150‑1 (MQ=255)
                                                                                                                                               aTTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTgg    >  1:397320/1‑150 (MQ=255)
                                                                                                                                                 tgttgCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTgggg  >  1:822360/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTGGGG  >  NC_000913/1457389‑1457683

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