Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,532 A→C 8.5% V12V (GTA→GTC) ‡ 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,5320AC 8.5% 195.1 / 8.7 69V12V (GTA→GTC) ‡paaFputative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base C (2/4);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTGG  >  NC_000913/1457389‑1457681
                                                                                                                                               |                                                                                                                                                     
cAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAttgttg                                                                                                                                                 >  2:224433/1‑150 (MQ=255)
   ttACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTg                                                                                                                                              <  1:224433/150‑1 (MQ=255)
     aCAGTCTGGAACCGGATGAACTGGCCGCAGTTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTTACTTTGACGCGATGGGGATGGCATGAGCGAACTTATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                            <  2:152667/150‑1 (MQ=255)
     aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                            <  1:1303482/150‑1 (MQ=255)
     aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                            <  2:342119/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                                                                                                       >  1:1399642/1‑150 (MQ=255)
                                          gtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaataat                                                                                                       <  2:1406128/150‑1 (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:1188355/1‑150 (MQ=255)
                                                                                                               aGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgtgtg                                  >  1:218180/1‑150 (MQ=255)
                                                                                                                        aTCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGt                         <  1:1443428/150‑1 (MQ=255)
                                                                                                                           gtcagccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGGTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAAt                      >  1:162356/1‑150 (MQ=255)
                                                                                                                             taagggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTc                                            >  2:697152/6‑126 (MQ=255)
                                                                                                                             taagggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTc                                            <  1:697152/121‑1 (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           <  2:495339/150‑1 (MQ=255)
                                                                                                                                      gcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCttttt           >  1:1463188/1‑150 (MQ=255)
                                                                                                                                        aaCGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTg         <  1:1111117/150‑1 (MQ=255)
                                                                                                                                           gAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCg      <  2:26164/150‑1 (MQ=255)
                                                                                                                                               aTTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTgg  >  1:397320/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGCTGG  >  NC_000913/1457389‑1457681

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: