Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,670,992 T→C 5.7% F29L (TTC→CTC)  hcaC → 3‑phenylpropionate dioxygenase, ferredoxin subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,670,9920TC 5.7% 349.0 / 8.4 107F29L (TTC→CTC) hcaC3‑phenylpropionate dioxygenase, ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base T (50/50);  new base C (3/3);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 110 total aligned reads displayed.

GAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTATGC  >  NC_000913/2670846‑2671141
                                                                                                                                                  |                                                                                                                                                     
gAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCa                                                                                                                                                    <  2:338568/150‑1 (MQ=255)
gAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCa                                                                                                                                                    <  2:338650/150‑1 (MQ=255)
     ggATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACgat                                                                                                                                               >  2:1686692/1‑148 (MQ=255)
     ggATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGtt                                                                                                                                               >  2:333128/1‑150 (MQ=255)
        tatCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTggc                                                                                                                                            <  2:1493195/150‑1 (MQ=255)
            gTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGgcggcg                                                                                                                                        >  2:1536828/1‑150 (MQ=255)
                   aTCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTa                                                                                                                                 <  2:78045/150‑1 (MQ=255)
                       aGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGcc                                                                                                                             >  1:1106806/1‑150 (MQ=255)
                         cGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCAt                                                                                                                           <  1:1012999/150‑1 (MQ=255)
                         cGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCAt                                                                                                                           <  2:1581105/150‑1 (MQ=255)
                           gTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTa                                                                                                                         >  2:1304763/1‑150 (MQ=255)
                           gTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTa                                                                                                                         >  2:143098/1‑150 (MQ=255)
                           gTAATCACTTCACATAACCTGAGTGAACTGTTCTGATGAATCGAATTTATGCGTTTCCAGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTa                                                                                                                         >  2:1246867/1‑150 (MQ=255)
                            tAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTaa                                                                                                                        <  1:1260852/150‑1 (MQ=255)
                                cACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGAt                                                                                                                    <  1:356359/150‑1 (MQ=255)
                                cACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGAt                                                                                                                    <  1:356393/150‑1 (MQ=255)
                                  cTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCg                                                                                                                  >  1:257264/1‑150 (MQ=255)
                                    tCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGtt                                                                                                                <  2:562193/150‑1 (MQ=255)
                                       cATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGTCATTAACGATCGTTGCa                                                                                                             <  1:255160/150‑1 (MQ=255)
                                       cATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCa                                                                                                             >  2:574965/1‑150 (MQ=255)
                                       cATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCa                                                                                                             >  1:1606558/1‑150 (MQ=255)
                                       cATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCa                                                                                                             >  1:558597/1‑150 (MQ=255)
                                         tAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGc                                                                                                           <  2:787292/150‑1 (MQ=255)
                                             cTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATg                                                                                                       >  2:629526/1‑150 (MQ=255)
                                             cTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGCTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATg                                                                                                       >  1:652028/1‑150 (MQ=255)
                                             cTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGCTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATg                                                                                                       >  2:1454474/1‑150 (MQ=255)
                                              tGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATgg                                                                                                      >  2:808132/1‑150 (MQ=255)
                                              tGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATgg                                                                                                      >  2:1661734/1‑150 (MQ=255)
                                              tGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATgg                                                                                                      >  2:1661896/1‑150 (MQ=255)
                                                   gtACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATg                                                                                                 >  1:200841/1‑150 (MQ=255)
                                                    tACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGc                                                                                                >  2:56110/1‑150 (MQ=255)
                                                    tACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGc                                                                                                >  1:577959/1‑150 (MQ=255)
                                                      cTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGt                                                                                              <  2:114681/150‑1 (MQ=255)
                                                       tGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTc                                                                                             <  2:1608319/150‑1 (MQ=255)
                                                       tGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGAGCGTTGCAGCCATGGTAATGCCTc                                                                                             <  2:926005/150‑1 (MQ=255)
                                                       tGTTCTGAGGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTc                                                                                             <  2:1266915/150‑1 (MQ=255)
                                                       tGTTCTGAGGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTc                                                                                             <  2:1266907/150‑1 (MQ=255)
                                                         ttCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa                                                                                           <  2:302335/150‑1 (MQ=255)
                                                            tgatgaATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGt                                                                                        >  1:466558/1‑150 (MQ=255)
                                                               tgaATCGAATTTATGCGTGACCCGTAGCGGATGTGCACGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGATCAACGTTGGCGGCGAGTTGTATTCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTAAg                                                                                     >  2:102024/1‑150 (MQ=255)
                                                                   tCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAgg                                                                                 >  2:349789/1‑150 (MQ=255)
                                                                   tCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAgg                                                                                 >  2:464995/1‑150 (MQ=255)
                                                                     gAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGt                                                                               >  2:322226/1‑150 (MQ=255)
                                                                      aaTTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTa                                                                              >  2:1122620/1‑150 (MQ=255)
                                                                      aaTTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTa                                                                              >  2:1122618/1‑150 (MQ=255)
                                                                           ctGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTgg                                                                         <  2:728855/149‑1 (MQ=255)
                                                                             gCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGaa                                                                       >  1:341862/1‑150 (MQ=255)
                                                                             gCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGaa                                                                       <  2:443583/150‑1 (MQ=255)
                                                                             gCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGaa                                                                       >  1:341980/1‑150 (MQ=255)
                                                                             gCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGaa                                                                       >  1:880499/1‑150 (MQ=255)
                                                                                         gCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACg                                                           >  1:363239/1‑150 (MQ=255)
                                                                                          cGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACgg                                                          <  2:1553000/150‑1 (MQ=255)
                                                                                            gATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTg                                                        >  2:1644562/1‑150 (MQ=255)
                                                                                              tgtgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGa                                                      <  1:759660/150‑1 (MQ=255)
                                                                                              tgtgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGa                                                      <  2:784328/150‑1 (MQ=255)
                                                                                              tgtgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGa                                                      <  2:1168869/150‑1 (MQ=255)
                                                                                              tgtgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTGGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGa                                                      <  2:363239/150‑1 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGATTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGTGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    <  1:26069/150‑1 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCATTGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    >  2:628891/1‑150 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCATTGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    >  2:9432/1‑150 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    >  2:722899/1‑150 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    >  2:1696720/1‑150 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    >  2:86152/1‑150 (MQ=255)
                                                                                                tgCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGt                                                    >  1:385569/1‑150 (MQ=255)
                                                                                                  ccGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGc                                                  >  2:1521269/1‑150 (MQ=255)
                                                                                                  ccGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGc                                                  >  1:1259561/1‑150 (MQ=255)
                                                                                                   cGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGcc                                                 <  2:1134518/150‑1 (MQ=255)
                                                                                                      aGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGc                                              <  1:1522574/150‑1 (MQ=255)
                                                                                                       gggTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCATCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCt                                             <  2:150424/150‑1 (MQ=255)
                                                                                                       gcgTGAGGCTCTCCGGATCGATACCTCGCCCCTCCTCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGCATCTGGAAGATGACGCCACGGTGGAGTGCCCGcc                                             <  1:1203300/148‑2 (MQ=255)
                                                                                                             ggCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACAcgc                                       >  1:1264537/1‑150 (MQ=255)
                                                                                                                tctcCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACAcgccgc                                    >  2:435650/1‑150 (MQ=255)
                                                                                                                  tcCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCa                                  >  1:202445/1‑150 (MQ=255)
                                                                                                                          gATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGc                          <  2:1350574/150‑1 (MQ=255)
                                                                                                                          gATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGc                          >  1:1500918/1‑150 (MQ=255)
                                                                                                                          gATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGc                          <  2:1637442/150‑1 (MQ=255)
                                                                                                                          gATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGc                          <  1:718768/150‑1 (MQ=255)
                                                                                                                            tACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCt                        >  2:285221/1‑150 (MQ=255)
                                                                                                                                  gCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAAc                  >  2:1053499/1‑150 (MQ=255)
                                                                                                                                   cccGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACg                 >  1:1367539/1‑150 (MQ=255)
                                                                                                                                    ccGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACgg                <  2:1264537/150‑1 (MQ=255)
                                                                                                                                    ccGTCATCGCCCTGCTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACgg                <  2:652028/150‑1 (MQ=255)
                                                                                                                                    ccGTCATCGCCCTGCTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACgg                <  1:1087778/150‑1 (MQ=255)
                                                                                                                                    ccGTCATCGCCCTGCTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACgg                <  1:1454474/150‑1 (MQ=255)
                                                                                                                                          tCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTGAGTAGGGTTTCAGCAAGATGTGGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAg          <  2:1505672/150‑1 (MQ=255)
                                                                                                                                          tCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAg          <  1:333128/150‑1 (MQ=255)
                                                                                                                                          tCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAg          <  1:1686692/150‑1 (MQ=255)
                                                                                                                                             cccTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGt       <  1:464995/150‑1 (MQ=255)
                                                                                                                                             cccTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGt       <  2:812124/150‑1 (MQ=255)
                                                                                                                                               cTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTa     <  1:919639/150‑1 (MQ=255)
                                                                                                                                               cTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTa     <  2:303701/150‑1 (MQ=255)
                                                                                                                                                tGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAt    <  1:808132/150‑1 (MQ=255)
                                                                                                                                                tGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAt    <  1:1661896/150‑1 (MQ=255)
                                                                                                                                                tGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAt    <  1:1661734/150‑1 (MQ=255)
                                                                                                                                                tGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAt    <  2:424656/150‑1 (MQ=255)
                                                                                                                                                 gTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAtg   >  2:1631141/1‑150 (MQ=255)
                                                                                                                                                 gTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAtg   >  2:1631136/1‑150 (MQ=255)
                                                                                                                                                  ttCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAtgc  <  2:1494643/150‑1 (MQ=255)
                                                                                                                                                  ttCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAtgc  <  1:1176247/150‑1 (MQ=255)
                                                                                                                                                  ttCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTAtgc  <  2:1583166/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTCAGAAGGGTATCTGGAAGATGACGCCACGGTGGAGTGCCCGCTACACGCCGCCAGTTTTTGCCTGAAAACGGGGAAAGCGTTATGC  >  NC_000913/2670846‑2671141

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