Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,986 A→T 6.6% T164S (ACC→TCC)  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,9860AT 6.6% 688.1 / 28.1 227T164S (ACC→TCC) paaFputative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (110/102);  new base T (7/8);  total (117/110)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-02

Only 100 of 231 total aligned reads displayed.

GCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCAG  >  NC_000913/1457837‑1458135
                                                                                                                                                        |                                                                                                                                                     
gCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCa                                                                                                                                                       >  2:1876043/1‑150 (MQ=255)
 cGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCTCTCTCGGCTTCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAc                                                                                                                                                      <  2:1072304/150‑1 (MQ=255)
 cGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAc                                                                                                                                                      <  2:966750/150‑1 (MQ=255)
 cGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAc                                                                                                                                                      <  1:860239/150‑1 (MQ=255)
  gATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCATAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCAcc                                                                                                                                                     <  2:1922999/150‑1 (MQ=255)
     gtggtgGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCt                                                                                                                                                  >  1:1391696/1‑150 (MQ=255)
      tggtggTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGTGATACTTTCTCCGCTc                                                                                                                                                 <  1:1839050/150‑1 (MQ=255)
      tggtggTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTc                                                                                                                                                 <  1:938133/150‑1 (MQ=255)
       ggtggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCa                                                                                                                                                <  2:558033/150‑1 (MQ=255)
           tttGCCGGAGATAACGCGCGTTTTGGGTTGCCGTAAATCACTCTCGGCATCATGCCTGGCGCAGTCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGTGATACTTTCTCCGCTCagcac                                                                                                                                            <  1:852301/149‑2 (MQ=255)
           gtTGCCGTATAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCTGAACGCAACTTTTAATCCGTAGTGTCGGTAAATCGTTAGACAGCAAAATGG‑‑TTCTGAGCGGTTATACTTTCTCCGCTCagcac                                                                                                                                            <  2:43551/150‑2 (MQ=255)
              gCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCagcaagca                                                                                                                                         >  2:586880/1‑150 (MQ=255)
                  aagagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagc                                                                                                                                     >  1:1681932/2‑150 (MQ=255)
                      gaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGc                                                                                                                                 >  2:1138807/1‑150 (MQ=255)
                      gaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGc                                                                                                                                 >  2:465991/1‑150 (MQ=255)
                         cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCagcaa                                                                                                                                            >  1:255339/1‑136 (MQ=255)
                         cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCagcaa                                                                                                                                            <  2:255339/136‑1 (MQ=255)
                               ttttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGt                                                                                                                        <  1:144862/150‑1 (MQ=255)
                               ttttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGt                                                                                                                        <  1:1715529/150‑1 (MQ=255)
                                tttGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGtt                                                                                                                       <  1:684236/150‑1 (MQ=255)
                                  tGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAg                                                                                                                     >  1:908182/1‑150 (MQ=255)
                                   gggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                                    >  1:106879/1‑150 (MQ=255)
                                        gCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGt                                                                                                               >  1:740856/1‑150 (MQ=255)
                                         ccGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTc                                                                                                              >  2:1673064/1‑150 (MQ=255)
                                         ccGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTc                                                                                                              >  1:1046214/1‑150 (MQ=255)
                                            gAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTc                                                                                                           <  2:1200362/150‑1 (MQ=255)
                                                tCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTccccc                                                                                                       >  2:483383/1‑149 (MQ=255)
                                                   ctctCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCg                                                                                                    >  2:1904726/1‑150 (MQ=255)
                                                     ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAt                                                                                                  >  2:1111565/1‑150 (MQ=255)
                                                     ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAt                                                                                                  >  2:817562/1‑150 (MQ=255)
                                                      tcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAtt                                                                                                 <  1:1528614/150‑1 (MQ=255)
                                                       cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                                <  1:1983747/150‑1 (MQ=255)
                                                       cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                                <  1:1449888/150‑1 (MQ=255)
                                                        ggCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTa                                                                                               >  2:627914/1‑150 (MQ=255)
                                                         gCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTaa                                                                                              <  2:548282/150‑1 (MQ=255)
                                                          catcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAc                                                                                             >  1:1748243/1‑150 (MQ=255)
                                                           ttcatGGCTGGCGGTGGCGGAACGCTAGGTTTATTGCGTAGTGTGGGTATTTCGTTAGCGAGCTATATGG‑‑TGCTGAGCGGAGAAAGTATGTCGGCTCTGGTAGGTCAGGAGGGCGGGCTGGTTAGCGACGTCTTCCGGAGCGATTTaagc                                                                                            <  1:804774/149‑3 (MQ=255)
                                                           atcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAcc                                                                                            <  1:1325002/150‑1 (MQ=255)
                                                           atcatGCCTGGCGCAGGCGGAACGCAACGTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                           >  1:1032946/1‑150 (MQ=255)
                                                            tcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                           <  1:1093592/150‑1 (MQ=255)
                                                             catGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCt                                                                                          >  2:702342/1‑150 (MQ=255)
                                                               tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                        <  1:702342/150‑1 (MQ=255)
                                                                gCCTGGCGCAGGCGGAACGCAAGGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGa                                                                                       <  2:1685571/150‑1 (MQ=255)
                                                                gCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGa                                                                                       >  1:1644667/1‑150 (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‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                                <  1:465988/150‑1 (MQ=255)
                                                                          ggCGGAACGCAAAGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCCTGCTGTGCTTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                                >  2:1915899/1‑150 (MQ=255)
                                                                               aaCGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCt                                                                        >  2:1828367/1‑150 (MQ=255)
                                                                                  gCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagcag                                                                                                                                   <  2:217352/88‑1 (MQ=255)
                                                                                     aCGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATc                                                                  <  1:1154811/150‑1 (MQ=255)
                                                                                      cGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCAGCg                                                                 <  1:302773/150‑1 (MQ=255)
                                                                                         ttAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGaaa                                                              >  2:791907/1‑150 (MQ=255)
                                                                                          tAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGaaaa                                                             >  2:248348/1‑150 (MQ=255)
                                                                                          tAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGaaaa                                                             >  1:1214629/1‑150 (MQ=255)
                                                                                             tCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATgg                                                          >  2:639551/1‑150 (MQ=255)
                                                                                             tCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATgg                                                          <  1:1502065/150‑1 (MQ=255)
                                                                                             tCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATgg                                                          >  2:1278567/1‑150 (MQ=255)
                                                                                               cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                                        <  1:1594406/150‑1 (MQ=255)
                                                                                               cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                                        >  1:1072215/1‑150 (MQ=255)
                                                                                                gTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCAc                                                       <  1:1917695/150‑1 (MQ=255)
                                                                                                 tAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGTGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACg                                                      <  1:138215/150‑1 (MQ=255)
                                                                                                 tAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACg                                                      <  2:863038/150‑1 (MQ=255)
                                                                                                    tgtCGGTAAATCGTTAGCCAGAAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCa                                                   >  2:944302/1‑150 (MQ=255)
                                                                                                        ggTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgc                                              >  2:1334502/1‑150 (MQ=255)
                                                                                                         gTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcgc                                              <  2:1451554/150‑1 (MQ=255)
                                                                                                             aTCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCt                                          >  1:533038/1‑150 (MQ=255)
                                                                                                               cGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTgg                                        <  1:777290/150‑1 (MQ=255)
                                                                                                                  tAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCttgtaa                                     >  2:995537/1‑145 (MQ=255)
                                                                                                                  tAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCttgtaa                                     >  2:1002263/1‑145 (MQ=255)
                                                                                                                  tAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCt                                     <  2:109531/150‑1 (MQ=255)
                                                                                                                     ccAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTAc                                  <  1:1876043/150‑1 (MQ=255)
                                                                                                                     ccAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTAc                                  <  1:48561/150‑1 (MQ=255)
                                                                                                                        gCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAg                               >  1:121804/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‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGG                            <  2:1446565‑M1/150‑1 (MQ=255)
                                                                                                                           aaaTGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGG                            >  2:1994306‑M1/1‑150 (MQ=255)
                                                                                                                                 g‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAG                       <  1:1828367‑M1/150‑1 (MQ=255)
                                                                                                                                    tGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGC                      <  1:1278567‑M1/150‑1 (MQ=255)
                                                                                                                                    tGCTGAGCGGAGAAAGTATCACCCCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGC                      <  1:1334654‑M1/150‑1 (MQ=255)
                                                                                                                                     gCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTgg                                        <  2:1863116/131‑1 (MQ=255)
                                                                                                                                      cTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa                    >  2:309551/1‑150 (MQ=255)
                                                                                                                                      cTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa                    >  2:1267973/1‑150 (MQ=255)
                                                                                                                                      cTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaa                    <  1:58257/150‑1 (MQ=255)
                                                                                                                                       tGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaag                   >  1:1130335/1‑150 (MQ=255)
                                                                                                                                        gAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaagc                  >  1:199586/1‑150 (MQ=255)
                                                                                                                                        gAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagcaagc                  >  1:329504/1‑150 (MQ=255)
                                                                                                                                          gCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAgcgc                >  1:1361417/1‑150 (MQ=255)
                                                                                                                                          gCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAgcgc                >  1:786450/1‑150 (MQ=255)
                                                                                                                                            ggAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTg              >  1:435118/1‑150 (MQ=255)
                                                                                                                                                aaaGTATCACCGCTCAGCAAGCACTCCAGGCCGGGCTGGTTAGCGACGTCTGCCCCAGCGAATTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGcc          <  1:1959083/150‑1 (MQ=255)
                                                                                                                                                 aaGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCa         >  1:91577/1‑150 (MQ=255)
                                                                                                                                                  aGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCACGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAg        <  1:992090/150‑1 (MQ=255)
                                                                                                                                                   gTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTACCCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCACCAGt       >  2:1603108/1‑150 (MQ=255)
                                                                                                                                                       cACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCa   >  2:1648957/1‑150 (MQ=255)
                                                                                                                                                        aCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCAg  >  1:49699/1‑150 (MQ=255)
                                                                                                                                                        aCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCCGCGCCAGTCGCAg  >  2:89274/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGG‑‑TGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCAG  >  NC_000913/1457837‑1458135

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: