Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,976 A→T 5.9% G160G (GGA→GGT 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,9760AT 5.9% 740.1 / 22.8 228G160G (GGA→GGTpaaFputative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (106/108);  new base T (7/7);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

GCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCG  >  NC_000913/1457827‑1458125
                                                                                                                                                        |                                                                                                                                                     
 cATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGgag                                                                                                                                                      <  2:904293/150‑1 (MQ=255)
    tgttgtGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGGTGCCGGAAATCACTCTAGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGaaa                                                                                                                                                   <  2:173006/150‑1 (MQ=255)
     gttgtGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGTAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATTG‑‑TGCTGAGCGGAGAAAg                                                                                                                                                  <  2:1018237/150‑1 (MQ=255)
      ttgtGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGt                                                                                                                                                 <  2:81174/150‑1 (MQ=255)
        gtgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTAt                                                                                                                                               >  2:1528614/1‑150 (MQ=255)
        gtgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTAt                                                                                                                                               <  1:1414583/150‑1 (MQ=255)
        gtgCGATGTGGTGGGTGCCGGAGAGCACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTAt                                                                                                                                               <  2:1677005/150‑1 (MQ=255)
         tgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATc                                                                                                                                              <  1:571122/150‑1 (MQ=255)
         tgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATc                                                                                                                                              >  1:1200362/1‑150 (MQ=255)
           cGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAc                                                                                                                                            <  2:966750/150‑1 (MQ=255)
           cGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAc                                                                                                                                            <  1:860239/150‑1 (MQ=255)
           cGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAc                                                                                                                                            >  2:1983747/1‑150 (MQ=255)
            gATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCATAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAcc                                                                                                                                           <  2:1922999/150‑1 (MQ=255)
            gATGTGGTGGTTGCCGGAGAGAACGCGCGGTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAcc                                                                                                                                           <  1:328845/150‑1 (MQ=255)
              tgtgGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGc                                                                                                                                         >  1:995930/1‑150 (MQ=255)
                tggtggTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGTGATACTTTCTCCGCTc                                                                                                                                       <  1:1839050/150‑1 (MQ=255)
                tggtggTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTc                                                                                                                                       <  1:1215895/150‑1 (MQ=255)
                tggtggTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTc                                                                                                                                       >  1:1072151/1‑150 (MQ=255)
                      ttGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCagcaag                                                                                                                                 <  1:33490/150‑1 (MQ=255)
                        gCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCagcaagca                                                                                                                               >  2:586880/1‑150 (MQ=255)
                         ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGATGAGCGGAGAAAGTATCACCGCTCAGCAAGcac                                                                                                                              <  2:266133/150‑1 (MQ=255)
                          cGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcaca                                                                                                                             <  2:251725/150‑1 (MQ=255)
                            aagagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagc                                                                                                                           >  1:1681932/2‑150 (MQ=255)
                                gaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGc                                                                                                                       <  2:242473/150‑1 (MQ=255)
                                 aaCGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGcc                                                                                                                      >  2:667431/1‑150 (MQ=255)
                                   cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCagcaa                                                                                                                                  >  1:255339/1‑136 (MQ=255)
                                     cgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGc                                                                                                                  >  1:1326967/1‑150 (MQ=255)
                                        gTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTgg                                                                                                               >  2:929081/1‑150 (MQ=255)
                                        gTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGAAGGCGGAACGCAACGTTTAATACGTAGTGTAGGTATATCGTTAGCCAGCAAAATCG‑‑AGCTGAGCGGATCAATTATCACCGCTCAGCAAGCACAGCAGGCCGGTCTgg                                                                                                               >  2:839520/1‑150 (MQ=255)
                                         ttttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGt                                                                                                              >  1:437705/1‑150 (MQ=255)
                                         ttttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGt                                                                                                              >  1:1451554/1‑150 (MQ=255)
                                         ttttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGt                                                                                                              <  1:1715529/150‑1 (MQ=255)
                                          tttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGtt                                                                                                             <  1:684236/150‑1 (MQ=255)
                                           ttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTa                                                                                                            <  1:520201/150‑1 (MQ=255)
                                             gggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                          <  2:1391696/150‑1 (MQ=255)
                                             gggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                          <  1:1782299/150‑1 (MQ=255)
                                              ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCg                                                                                                         <  1:1972598/150‑1 (MQ=255)
                                                 tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACg                                                                                                      <  1:1048066/150‑1 (MQ=255)
                                                 tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACg                                                                                                      <  2:331917/150‑1 (MQ=255)
                                                 tGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGGTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACg                                                                                                      <  1:793391/150‑1 (MQ=255)
                                                  gCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGt                                                                                                     >  1:740856/1‑150 (MQ=255)
                                                      gAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTc                                                                                                 <  2:1200362/150‑1 (MQ=255)
                                                        aaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTccc                                                                                               >  2:1136508/1‑150 (MQ=255)
                                                             ctctCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCg                                                                                          >  2:1904726/1‑150 (MQ=255)
                                                               ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGTGATACTTTCTCCGCTCagca                                                                                                                                   >  1:1448674/1‑107 (MQ=38)
                                                               ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAt                                                                                        >  2:817562/1‑150 (MQ=255)
                                                                tcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAtt                                                                                       <  1:1528614/150‑1 (MQ=255)
                                                                   gCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTaa                                                                                    <  2:548282/150‑1 (MQ=255)
                                                                     atcatGCCTGGCGCAGGCGGAACGCAACGTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                 >  1:1032946/1‑150 (MQ=255)
                                                                      tcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                 <  2:1110448/150‑1 (MQ=255)
                                                                       catGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCt                                                                                >  2:702342/1‑150 (MQ=255)
                                                                        atGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTc                                                                               >  1:1368239/1‑150 (MQ=255)
                                                                         tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                              <  1:1138807/150‑1 (MQ=255)
                                                                         tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGGAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCATAGCGATTTTACCCTCg                                                                              <  2:1629634/150‑1 (MQ=255)
                                                                          gCCTGGCGCAGGCGGAACGCAAGGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGa                                                                             <  2:1685571/150‑1 (MQ=255)
                                                                          gCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGa                                                                             >  1:1644667/1‑150 (MQ=255)
                                                                            cTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAAt                                                                           <  1:1357765/150‑1 (MQ=255)
                                                                                cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATagga                                                                       >  1:43139/1‑147 (MQ=255)
                                                                                cgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGc                                                                       >  2:1795591/1‑150 (MQ=255)
                                                                                 gcAGGCTGAACGCAACGTTTAATCCGTAGTTTCGGTAACTCGTTCGCCAGCAAACTGG‑‑TGCTGAGCGGAGATAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGAGGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                      <  2:914036/150‑1 (MQ=255)
                                                                                 gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGaa                                                                            <  1:1272125/144‑1 (MQ=255)
                                                                                 gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGaa                                                                            >  2:1272125/1‑144 (MQ=255)
                                                                                 gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                      <  1:465988/150‑1 (MQ=255)
                                                                                    ggCGGAACGCAAAGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCCTGCTGTGCTTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                      >  2:1915899/1‑150 (MQ=255)
                                                                                         aaCGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCt                                                              <  1:1102936/150‑1 (MQ=255)
                                                                                         aaCGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCt                                                              >  2:1828367/1‑150 (MQ=255)
                                                                                            gCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagcag                                                                                                                         >  1:217352/1‑88 (MQ=255)
                                                                                            gCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagcag                                                                                                                         <  2:217352/88‑1 (MQ=255)
                                                                                             cAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                              <  1:829041/130‑1 (MQ=255)
                                                                                             cAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                          >  2:47278/1‑150 (MQ=255)
                                                                                                      aTCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATg                                                 <  1:884142/150‑1 (MQ=255)
                                                                                                      aTCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAACGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATg                                                 >  2:635431/1‑150 (MQ=255)
                                                                                                       tCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATgg                                                <  1:1502065/150‑1 (MQ=255)
                                                                                                         cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                              <  1:1594406/150‑1 (MQ=255)
                                                                                                                  ggTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGCCATCGAAAATGGAACGTCACTcg                                     >  2:1255428/1‑150 (MQ=255)
                                                                                                                  ggTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcg                                     <  2:442926/150‑1 (MQ=255)
                                                                                                                  ggTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcg                                     >  2:48561/1‑150 (MQ=255)
                                                                                                                  ggTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcg                                     >  1:704909/1‑150 (MQ=255)
                                                                                                                  ggTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgc                                    >  2:1334502/1‑150 (MQ=255)
                                                                                                                   gTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgc                                    <  2:1451554/150‑1 (MQ=255)
                                                                                                                     aaaTCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgccg                                  >  1:1933776/1‑150 (MQ=255)
                                                                                                                     aaaTCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTACCCCAGCGATTTAACCCACGAAACCGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgccg                                  >  1:1083078/1‑150 (MQ=255)
                                                                                                                         cGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTgg                              <  1:777290/150‑1 (MQ=255)
                                                                                                                            tAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCttgtaa                           >  2:995537/1‑145 (MQ=255)
                                                                                                                            tAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCt                           <  2:109531/150‑1 (MQ=255)
                                                                                                                              gCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                         >  1:1305377/1‑150 (MQ=255)
                                                                                                                               ccAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTAc                        <  1:48561/150‑1 (MQ=255)
                                                                                                                                cAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACa                       >  1:985753/1‑150 (MQ=255)
                                                                                                                                   cAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGc                    <  1:1670393/150‑1 (MQ=255)
                                                                                                                                   cAAAACGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGc                    <  1:635431/150‑1 (MQ=255)
                                                                                                                                    aaaaTGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCg                   <  1:1734444/150‑1 (MQ=255)
                                                                                                                                     aaaTGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCCCGTCACTCGCCGCTGGCCTTACAAGCGG                  >  2:1883960‑M1/1‑150 (MQ=255)
                                                                                                                                     aaaTGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGG                  <  2:576088‑M1/150‑1 (MQ=255)
                                                                                                                                           g‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAG             <  1:915084‑M1/150‑1 (MQ=255)
                                                                                                                                           g‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAG             <  1:1828367‑M1/150‑1 (MQ=255)
                                                                                                                                              tGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGC            <  1:1904726‑M1/150‑1 (MQ=255)
                                                                                                                                               gCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTgg                              <  2:1863116/131‑1 (MQ=255)
                                                                                                                                                cTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa          >  1:1391337/1‑150 (MQ=255)
                                                                                                                                                cTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa          <  1:58257/150‑1 (MQ=255)
                                                                                                                                                  gAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaagc        >  1:199586/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GCATTGTTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCG  >  NC_000913/1457827‑1458125

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: