Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,981 G→C 5.2% S162T (AGT→ACT)  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,9810GC 5.2% 306.8 / 6.5 97S162T (AGT→ACT) paaF2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base G (43/49);  new base C (2/3);  total (45/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 = 9.31e-01

TTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAG  >  NC_000913/1457833‑1458129
                                                                                                                                                    |                                                                                                                                                    
ttgtGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGt                                                                                                                                                     <  1:679099/150‑1 (MQ=255)
  gtgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAATTCGTTAGCCAGCAAAATGGTGCTGAGCGGAGCTAGTAt                                                                                                                                                   >  1:261990/1‑150 (MQ=255)
   tgCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATc                                                                                                                                                  >  1:1728409/1‑150 (MQ=255)
    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:469719/1‑150 (MQ=255)
                             cgcgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCgg                                                                                                                        >  2:748097/1‑150 (MQ=255)
                               cgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGc                                                                                                                      >  2:558109/1‑150 (MQ=255)
                               cgcgTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGc                                                                                                                      >  2:1547013/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:1578212/150‑1 (MQ=255)
                                                 aaaTCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTcc                                                                                                    <  1:1556793/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:1798014/1‑150 (MQ=255)
                                                           cGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGAttt                                                                                          <  1:584458/150‑1 (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:1414004/150‑1 (MQ=255)
                                                                           gcAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGcc                                                                          <  2:952471/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)
                                                                                                                                                    |                                                                                                                                                    
TTGTGCGATGTGGTGGTTGCCGGAGAGAACGCGCGTTTTGGGTTGCCGGAAATCACTCTCGGCATCATGCCTGGCGCAGGCGGAACGCAACGTTTAATCCGTAGTGTCGGTAAATCGTTAGCCAGCAAAATGGTGCTGAGCGGAGAAAGTATCACCGCTCAGCAAGCACAGCAGGCCGGGCTGGTTAGCGACGTCTTCCCCAGCGATTTAACCCTCGAATACGCCTTACAGCTGGCATCGAAAATGGCACGTCACTCGCCGCTGGCCTTACAAGCGGCAAAGCAAGCGCTGCGCCAG  >  NC_000913/1457833‑1458129

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: ‑