Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,983 A→T 8.2% I163F (ATC→TTC)  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,9830AT 8.2% 1018.6 / 54.8 333I163F (ATC→TTC) paaFputative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (161/144);  new base T (12/16);  total (173/160)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 336 total aligned reads displayed.

TGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCG  >  NC_000913/1457834‑1458132
                                                                                                                                                     |                                                                                                                                                     
tgtgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTa                                                                                                                                                       <  1:550185/150‑1 (MQ=255)
 gtgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTAt                                                                                                                                                      >  1:1764522/1‑150 (MQ=255)
 gtgCGAGGTGGTGGGTGCCGGAGAGAACGCGCGGTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTAt                                                                                                                                                      <  1:1372534/150‑1 (MQ=255)
       tgtgGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGc                                                                                                                                                <  1:288175/150‑1 (MQ=255)
        gtggtgGTTTCCTGAAAGAACACACGTTAATAGTTTCAGGAAATCACTCTACGCATAATCAATGACGCAACCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCt                                                                                                                                               >  2:395410/1‑150 (MQ=255)
        gtggtgGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCAc                                                                                                                                                   <  2:109994/146‑1 (MQ=255)
         gggtggtTGCCGGAGTGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCTGGGGGAACGCAACGTTTAATCCGTAGGGTCGGTAAATGGTTAGCCAGCAAAATGGTGATGAGCGTAGAAAGTATCACCGCTc                                                                                                                                              <  2:1744911/149‑1 (MQ=255)
          ggtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                             >  1:1974775/1‑150 (MQ=255)
          ggtggtTGCCGGAGAGAACGAGCGTATTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                             >  2:2121114/1‑150 (MQ=255)
           gtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCag                                                                                                                                            <  1:1090209/150‑1 (MQ=255)
            gggtTGCCGGAGAGTACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGTTAAATCGTTAGCCTGCAAAATGGTGCTGAGCGGAGAAAGTTTCACCGCTCagc                                                                                                                                           <  2:2375847/149‑1 (MQ=255)
            gggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagc                                                                                                                                           <  1:1223805/149‑1 (MQ=255)
              gtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagcac                                                                                                                                         <  1:2257028/150‑2 (MQ=255)
               ttGCCGGAGAGAACGCGCGTTTTGGTTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGACCGCAACGTTTAATCCTTAGTGTCGGGAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagcacc                                                                                                                                        <  2:1593494/150‑3 (MQ=255)
                tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       >  1:1000551/1‑150 (MQ=255)
                tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       <  1:951836/150‑1 (MQ=255)
                tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                       <  2:255089/150‑1 (MQ=255)
                  ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcac                                                                                                                                     >  1:2410383/1‑150 (MQ=255)
                  ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcac                                                                                                                                     >  2:2242743/1‑150 (MQ=255)
                   cGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcaca                                                                                                                                    <  2:2184777/150‑1 (MQ=255)
                    ggAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcag                                                                                                                                   >  2:1462904/1‑150 (MQ=255)
                    ggAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGAGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcag                                                                                                                                   >  1:1843104/1‑150 (MQ=255)
                          aaCGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGcc                                                                                                                             <  1:1264409/150‑1 (MQ=255)
                           aCGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCg                                                                                                                            >  2:149317/1‑150 (MQ=255)
                            cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                           >  2:852427/1‑150 (MQ=255)
                                 gTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTgg                                                                                                                      >  2:587630/1‑150 (MQ=255)
                                 gTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTgg                                                                                                                      <  1:1347820/150‑1 (MQ=255)
                                 gTTTTGGGTTGCCGGAAATCACTCTCGCCATAATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTgg                                                                                                                      >  2:1942393/1‑150 (MQ=255)
                                  ttttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGt                                                                                                                     <  1:1450950/150‑1 (MQ=255)
                                     tGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAg                                                                                                                  <  2:311336/150‑1 (MQ=255)
                                       ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCg                                                                                                                <  1:2232686/150‑1 (MQ=255)
                                         ttGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGAc                                                                                                              >  2:1515895/1‑150 (MQ=255)
                                          tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACg                                                                                                             <  2:69515/150‑1 (MQ=255)
                                             cGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCt                                                                                                          >  2:1704612/1‑150 (MQ=255)
                                             cGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCt                                                                                                          >  1:2478137/1‑150 (MQ=255)
                                                  aTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcccc                                                                                                     >  2:1312399/1‑150 (MQ=255)
                                                    cACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAg                                                                                                   >  2:1903749/1‑150 (MQ=255)
                                                     aCTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGc                                                                                                  <  2:366979/150‑1 (MQ=255)
                                                          cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                             <  1:2092499/150‑1 (MQ=255)
                                                          cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                             <  1:1837784/150‑1 (MQ=255)
                                                          cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                             <  2:640294/150‑1 (MQ=255)
                                                           ggCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTa                                                                                            >  1:1701117/1‑150 (MQ=255)
                                                             catcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAc                                                                                          <  1:574163/150‑1 (MQ=255)
                                                               tcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                        >  1:1030912/1‑150 (MQ=255)
                                                                  tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                     >  1:1359658/1‑150 (MQ=255)
                                                                   gCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGa                                                                                    <  2:761817/150‑1 (MQ=255)
                                                                     cTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAAt                                                                                  <  2:104516/150‑1 (MQ=255)
                                                                      tGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATa                                                                                 >  1:1308015/1‑150 (MQ=255)
                                                                         cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGc                                                                              >  2:920098/1‑150 (MQ=255)
                                                                         cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGc                                                                              >  1:1421028/1‑150 (MQ=255)
                                                                          gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                             <  1:2016281/150‑1 (MQ=255)
                                                                            aGGCAGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCtt                                                                           >  2:164114/1‑150 (MQ=255)
                                                                              gCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTAc                                                                         >  1:810662/1‑150 (MQ=255)
                                                                               cGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACa                                                                        <  1:2340402/150‑1 (MQ=255)
                                                                               cGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACa                                                                        <  2:1284535/150‑1 (MQ=255)
                                                                                      cAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                                 <  2:1846845/150‑1 (MQ=255)
                                                                                        aCGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATc                                                               <  2:1026937/150‑1 (MQ=255)
                                                                                         cGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCg                                                              >  1:1802479/1‑150 (MQ=255)
                                                                                             tAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGaaaa                                                          >  1:213585/1‑150 (MQ=255)
                                                                                              aaTCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAAt                                                         <  2:454458/150‑1 (MQ=255)
                                                                                              aaTCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAAt                                                         >  1:2389674/1‑150 (MQ=255)
                                                                                                 ccGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGc                                                      <  1:916792/150‑1 (MQ=255)
                                                                                                  cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                                     <  2:507231/150‑1 (MQ=255)
                                                                                                  cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                                     >  2:2149467/1‑150 (MQ=255)
                                                                                                   gTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCAc                                                    <  1:2043985/150‑1 (MQ=255)
                                                                                                      gtgtCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTc                                                 <  1:1572066/150‑1 (MQ=255)
                                                                                                       tgtCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCa                                                <  1:2312645/150‑1 (MQ=255)
                                                                                                       tgtCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAACGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCa                                                <  1:366658/150‑1 (MQ=255)
                                                                                                        ttCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGAGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCCCGTCAc                                               >  2:1177416/2‑150 (MQ=255)
                                                                                                        gtCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCAc                                               >  2:2266944/1‑150 (MQ=255)
                                                                                                         tCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACt                                              <  1:351307/150‑1 (MQ=255)
                                                                                                            gTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgc                                           >  2:2407560/1‑150 (MQ=255)
                                                                                                                  cGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTgg                                     <  2:9589/150‑1 (MQ=255)
                                                                                                                   gTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGc                                    <  1:88925/150‑1 (MQ=255)
                                                                                                                    ttAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGcc                                   <  2:1555892/150‑1 (MQ=255)
                                                                                                                       gCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                                >  2:1010123/1‑150 (MQ=255)
                                                                                                                       gCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                                >  1:708841/1‑150 (MQ=255)
                                                                                                                       gCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                                >  2:1353732/1‑150 (MQ=255)
                                                                                                                       gCCACCAAAAACGTGCTGATCGGAGAAAGTATAACCGCTCCGCAAGCACAGCAGGACGGGCTGGATAGCCAGGTCTTCCACAGCCATTTAACCCTCGAATAAGCCTTACAGCTGGCATCGAAAATGGCACGTCACACGCAGCTGGCCTTa                                >  2:1409958/1‑150 (MQ=255)
                                                                                                                             aaaaTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCg                          <  2:1342229/150‑1 (MQ=255)
                                                                                                                             aaaaTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCg                          <  1:1506693/150‑1 (MQ=255)
                                                                                                                                  cttgcTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAA                     >  1:2119805‑M1/3‑150 (MQ=255)
                                                                                                                                    tGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGC                   <  1:39258‑M1/150‑1 (MQ=255)
                                                                                                                                      cTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa                 >  2:2041386/1‑150 (MQ=255)
                                                                                                                                       tGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTAc                               >  2:122854/1‑135 (MQ=255)
                                                                                                                                       tGAGCGGGGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaag                >  1:782321/1‑150 (MQ=255)
                                                                                                                                       tGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaag                >  2:2437879/1‑150 (MQ=255)
                                                                                                                                        gAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaagc               >  2:2279827/1‑150 (MQ=255)
                                                                                                                                         aGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAgcg              <  2:66389/150‑1 (MQ=255)
                                                                                                                                         aGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAgcg              >  2:1984853/1‑150 (MQ=255)
                                                                                                                                          gCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAgcgc             >  1:489771/1‑150 (MQ=255)
                                                                                                                                          gCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCACCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAgcgc             >  2:328912/1‑150 (MQ=255)
                                                                                                                                             gagaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgc          <  1:654807/150‑1 (MQ=255)
                                                                                                                                             gagaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgc          <  1:1496612/150‑1 (MQ=255)
                                                                                                                                               gaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgcgc        <  2:557785/150‑1 (MQ=255)
                                                                                                                                               gaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgcgc        <  1:1923500/150‑1 (MQ=255)
                                                                                                                                               gaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTgcgc        >  2:1773509/1‑150 (MQ=255)
                                                                                                                                                aaaGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGcc       <  1:1903749/150‑1 (MQ=255)
                                                                                                                                                aaaGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGcc       >  2:2074313/1‑150 (MQ=255)
                                                                                                                                                    tATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTc   <  1:1695007/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCG  >  NC_000913/1457834‑1458132

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: