Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,986 A→T 5.1% T164S (ACC→TCC)  paaF → 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 5.1% 250.8 / 3.1 94T164S (ACC→TCC) paaF2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base A (38/48);  new base T (2/3);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

GCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCAG  >  NC_000913/1457837‑1458135
                                                                                                                                                     |                                                                                                                                                     
gCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCa                                                                                                                                                       >  2:1749316/1‑150 (MQ=255)
    tgtgGTGGGTGCCGGAGAGAACGCGCGTTTTGGGGTGCCGGAAATCACTCTAGGCAGCATGCCTGGCGCAGGCGGAACGCAACGTATAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGc                                                                                                                                                   <  1:1141974/150‑1 (MQ=255)
         tggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagc                                                                                                                                              <  2:476402/150‑1 (MQ=255)
         tggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACACTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAAGCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagc                                                                                                                                              <  1:103324/150‑1 (MQ=255)
          ggtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagca                                                                                                                                             >  1:683829/1‑150 (MQ=255)
          ggatGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGAAACGTTTATTCCGTATTGTCTGTAAATCGTTAGCCAGCAAAATCGTTCTGAGCGGTGATACTTTCTCCGCTCagca                                                                                                                                             <  1:390852/147‑1 (MQ=255)
           gtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTCGCCAGCAAAATGGTGCAGAGCGGAGAAAGTATCCCCGCTCagcaa                                                                                                                                            >  2:1181178/1‑150 (MQ=255)
           gtTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGTGATACTTTCTCCGCTCagcac                                                                                                                                            <  2:1088252/150‑2 (MQ=255)
             tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                          >  2:979009/1‑150 (MQ=255)
             tGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCagcaagc                                                                                                                                          <  2:972935/150‑1 (MQ=255)
               ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcac                                                                                                                                        <  2:758703/150‑1 (MQ=255)
               ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGcac                                                                                                                                        <  2:1087816/150‑1 (MQ=255)
               ccGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCACTCAGCAAGcac                                                                                                                                        <  2:645660/150‑1 (MQ=255)
               ccGGAGAGAACGCGCGGTTTGGGTTGCCGGAAATCACTCTCGTCATCCTGCCTGGCGCAGGCGGACCGCAACGTTTAAACCGTAGTTACGGTAACTCGGTAGCCAGCAACATGGTGCTGAGCGGCGAAAGAATCACCGCTAATCCAGcac                                                                                                                                        >  1:106386/1‑150 (MQ=255)
                 ggAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcag                                                                                                                                      <  2:894766/150‑1 (MQ=255)
                  gagagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagc                                                                                                                                     <  1:1117597/150‑1 (MQ=255)
                   agagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagca                                                                                                                                    >  2:1440557/1‑150 (MQ=255)
                   agagaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCAcagca                                                                                                                                    <  1:367164/150‑1 (MQ=255)
                     agaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAgg                                                                                                                                  <  1:1061399/150‑1 (MQ=255)
                      gaACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTATATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGc                                                                                                                                 >  2:1856686/1‑150 (MQ=255)
                         cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                              >  2:748097/1‑150 (MQ=255)
                         cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                              >  2:469719/1‑150 (MQ=255)
                           cgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGc                                                                                                                            >  2:1547013/1‑150 (MQ=255)
                           cgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGc                                                                                                                            >  2:558109/1‑150 (MQ=255)
                           cgcgCTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGc                                                                                                                            <  2:654460/150‑1 (MQ=255)
                                   gggTTGCCGGACAGCACTCTCGGCATCATGCCTGGCGCACGCGGAACGCAACGTGTAATCCGTAGTGGCGGCAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGGATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGc                                                                                                                    <  2:1347703/150‑1 (MQ=255)
                                    ggTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCg                                                                                                                   <  1:520570/150‑1 (MQ=255)
                                     gTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGa                                                                                                                  >  2:927979/1‑150 (MQ=255)
                                     gTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGa                                                                                                                  >  2:1074057/1‑150 (MQ=255)
                                     gTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAAAGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTTAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGa                                                                                                                  >  1:1154954/1‑150 (MQ=255)
                                       tGCCTGAAATCACTCTCGGCTTCATGCCTGGCGCAGGCGGAACGAAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGGTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACg                                                                                                                <  2:1235415/150‑1 (MQ=255)
                                           ggAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCtt                                                                                                            >  1:585646/1‑150 (MQ=255)
                                             aaaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcc                                                                                                          <  1:1556793/150‑1 (MQ=255)
                                             aaaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcc                                                                                                          <  1:1578212/150‑1 (MQ=255)
                                             aaaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcc                                                                                                          <  1:748097/150‑1 (MQ=255)
                                              aaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTccc                                                                                                         >  1:1638043/1‑150 (MQ=255)
                                                 cACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAg                                                                                                      <  1:284734/150‑1 (MQ=255)
                                                  aCTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACATCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGc                                                                                                     <  2:372714/150‑1 (MQ=255)
                                                    tctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGa                                                                                                   <  2:1359851/150‑1 (MQ=255)
                                                     ctcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGATAGTATCACCGCTCAGCAAGCAAAGCAGGCTGGGCTGGTTAGCGACGTCTTCCCCAGCGAt                                                                                                  >  2:340936/1‑150 (MQ=255)
                                                     ctcGGCATAATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAAAGGTGCTTAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGGTAGCGACGTCTTCCCCAGCGAt                                                                                                  >  2:654390/1‑150 (MQ=255)
                                                      tcGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGTTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAtt                                                                                                 <  1:1114334/150‑1 (MQ=255)
                                                       cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                                <  1:584458/150‑1 (MQ=255)
                                                       cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                                >  1:1798014/1‑150 (MQ=255)
                                                         gCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATTGTGCTGCGCGGAGAAAGTATCACCCCTCAGCAAGCACAGCAGGCCGGGCTGTATAGCGACGTCTTCCCCAGCGATTTaa                                                                                              >  1:971713/1‑150 (MQ=255)
                                                          catcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAc                                                                                             <  2:1728409/150‑1 (MQ=255)
                                                            tcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                           >  2:1577200/1‑150 (MQ=255)
                                                            tcatGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAAccc                                                                                           <  2:1297864/150‑1 (MQ=255)
                                                              atGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTc                                                                                         <  2:461278/150‑1 (MQ=255)
                                                               tGCCTGGCGCAGGCGGAACGCAAGGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCTGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                        <  1:15937/150‑1 (MQ=255)
                                                               tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                        <  2:1091604/150‑1 (MQ=255)
                                                               tGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGCGAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCg                                                                                        <  2:95967/150‑1 (MQ=255)
                                                                    ggCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTACCCCTCGAATAc                                                                                   >  1:1525846/1‑150 (MQ=255)
                                                                     gcgcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCCCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACg                                                                                  >  2:1068184/1‑150 (MQ=255)
                                                                     gagcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCTCTAACCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCAAGCGATTTAACCCTCGAATACg                                                                                  >  1:1342840/3‑150 (MQ=255)
                                                                      ggcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGc                                                                                 <  2:1841798/149‑1 (MQ=255)
                                                                       gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                                <  1:196305/150‑1 (MQ=255)
                                                                       gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                                <  2:952471/150‑1 (MQ=255)
                                                                       gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                                <  2:1414004/150‑1 (MQ=255)
                                                                         aGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCtt                                                                              >  2:633455/1‑150 (MQ=255)
                                                                          ggCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTa                                                                             >  2:1696983/1‑150 (MQ=255)
                                                                                 cGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTgg                                                                      <  2:1311301/150‑1 (MQ=255)
                                                                                   cAACGTTTAATCCGTAGTGTCGGTACATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                                    <  2:688587/150‑1 (MQ=255)
                                                                                   cAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                                    >  1:1184537/1‑150 (MQ=255)
                                                                                   cAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCa                                                                    >  2:1823296/1‑150 (MQ=255)
                                                                                      cGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCg                                                                 <  2:1663489/150‑1 (MQ=255)
                                                                                          tAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGaaaa                                                             >  1:748705/1‑150 (MQ=255)
                                                                                             tCCGTAGCGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAAAgg                                                          <  2:688950/150‑1 (MQ=255)
                                                                                               cGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCa                                                        <  1:927979/150‑1 (MQ=255)
                                                                                                gTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCAc                                                       >  1:1391370/1‑150 (MQ=255)
                                                                                                     gtCGGTAAATCGTTTGCCAGCAAATTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGACTACGCCTTACAGCTGGCATCGAAAATGGCACGTCAc                                                  <  2:264791/150‑1 (MQ=255)
                                                                                                     gtCGGTAAATCGATAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCAc                                                  <  1:472116/150‑1 (MQ=255)
                                                                                                        ggTAAATCGTTAGTCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCCGCAAGCACAGCAGGCCGGGCTGTTTAGCGACGTCTTCCCCAGCGAGTTAACCCGCGAATACGCCTTACAGCGGGCATCGAAAATGGCACGTCACTcc                                               >  2:310820/1‑149 (MQ=255)
                                                                                                        ggTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTcg                                               >  1:1715704/1‑150 (MQ=255)
                                                                                                            aaTCGTTAGCCAGCAAAATGGTGCTGAGCCGAGAAAGTATCACCGCTCAGCAAGCACAGCAGTCCGGGCTGGATAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCATTAAAGCTGGCATCGAAAATGCAACGTCACTcgccgc                                           >  1:169554/1‑150 (MQ=255)
                                                                                                               cGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTgg                                        <  1:1429853/150‑1 (MQ=255)
                                                                                                                gTTAGCCAGCAACATGGTGCTGAGCGGTGAAAGTATCACCGCTCAGCAAGCAGAGCTGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGc                                       <  2:1680486/150‑1 (MQ=255)
                                                                                                                 ttAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCATTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGcc                                      <  2:1265137/150‑1 (MQ=255)
                                                                                                                 ttAGCCAGCAAAATGGTGCTGAGAGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACTCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGcc                                      <  1:54011/150‑1 (MQ=255)
                                                                                                                  tAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCt                                     >  2:1061541/1‑150 (MQ=255)
                                                                                                                   aGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTACCCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCCCTCGCCGCTGGCCtt                                    >  1:337867/1‑150 (MQ=255)
                                                                                                                    gCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTa                                   >  1:1622418/1‑150 (MQ=255)
                                                                                                                      cAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACa                                 >  1:403566/1‑150 (MQ=255)
                                                                                                                       aGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCGAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACaa                                <  2:21496/150‑1 (MQ=255)
                                                                                                                         cAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGCCGTCTTCCCCAGCGATTTACCCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGc                              >  1:74638/1‑150 (MQ=255)
                                                                                                                           aaaTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCCTACAAGCgg                            >  2:137753/1‑150 (MQ=255)
                                                                                                                             aTGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCa                          >  1:1292045/1‑150 (MQ=255)
                                                                                                                                  gCTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCAGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagca                     >  1:1397553/1‑150 (MQ=255)
                                                                                                                                  gCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAaagca                     <  1:1074057/150‑1 (MQ=255)
                                                                                                                                   cTGAGCGGTGATACTTTCTCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCTaagcaa                    >  1:1016489/1‑150 (MQ=255)
                                                                                                                                            gaAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCCgcgc           >  2:460482/1‑150 (MQ=255)
                                                                                                                                               aGTATCCCCGCTCAGCAAGCACAGCCGGCCGGGGTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAg        <  2:131029/150‑1 (MQ=255)
                                                                                                                                               aGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAg        <  1:618148/150‑1 (MQ=255)
                                                                                                                                               aGTATCACCGCGCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGGCACTCGCCGCGGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAg        <  1:1160929/150‑1 (MQ=255)
                                                                                                                                                   tCACCGCTAAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGc    <  2:1342840/150‑1 (MQ=255)
                                                                                                                                                     aCCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCAg  <  2:1153925/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAGTCGCAG  >  NC_000913/1457837‑1458135

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑