Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,570,243 A→G 5.8% I352V (ATT→GTT)  pcaF‑I → beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,570,2430AG 5.8% 203.2 / 2.1 69I352V (ATT→GTT) pcaF‑Ibeta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base G (2/2);  total (35/34)
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.16e-01

GCAAACTGACCGAGCGTCTGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATC  >  NC_002947/1570097‑1570380
                                                                                                                                                  |                                                                                                                                          
gcAAACTGACCGAGCGTCTGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCAtt                                                                                                                                          >  2:337138/1‑149 (MQ=255)
        aCCGAGCGTCTGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTg                                                                                                                                    <  4:87739/147‑1 (MQ=255)
        aCCGAGCGTCTGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAAAGGCGGTGCCATTGCCCTgg                                                                                                                                   <  6:303325/148‑1 (MQ=255)
          cGAGCGTCTGGGGACAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCGCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGcc                                                                                                                                >  6:247483/1‑149 (MQ=255)
                tCTGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCAccccc                                                                                                                          >  7:53771/1‑149 (MQ=255)
                  tGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTg                                                                                                                        >  6:276347/1‑149 (MQ=255)
                    gggATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTggg                                                                                                                      >  1:299162/1‑149 (MQ=255)
                    gggATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTggg                                                                                                                      >  2:351729/1‑149 (MQ=255)
                          gCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGa                                                                                                                 <  2:299162/148‑1 (MQ=255)
                           cGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGc                                                                                                               >  7:330138/1‑149 (MQ=255)
                              tAAGTGATTTCGACGTGCTCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGt                                                                                                            <  3:272207/149‑1 (MQ=255)
                              tAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGt                                                                                                            >  8:302663/1‑149 (MQ=255)
                                  tGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAGTGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCAc                                                                                                        >  4:319628/1‑149 (MQ=255)
                                  tGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCaa                                                                                                        >  3:75176/1‑148 (MQ=255)
                                   gATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGACCTGGGCCACCCCCTGGGCATGAGCGGTGCACg                                                                                                       >  5:20348/1‑149 (MQ=255)
                                    aTTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGc                                                                                                      >  2:310208/1‑149 (MQ=255)
                                    aTTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGGGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACg                                                                                                       <  5:297360/148‑1 (MQ=255)
                                      ttCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCt                                                                                                    >  8:1924/1‑149 (MQ=255)
                                      ttCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACAATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCt                                                                                                    >  1:102153/1‑149 (MQ=255)
                                       tCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTg                                                                                                   >  2:327404/1‑149 (MQ=255)
                                        cGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTgg                                                                                                  >  3:134377/1‑149 (MQ=255)
                                         gACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACTCCCTGGGCATGAGCGGTGCACGCCTgg                                                                                                  <  8:53771/148‑1 (MQ=255)
                                         gACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTgg                                                                                                  <  8:330138/148‑1 (MQ=255)
                                              gATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTactga                                                                                            <  5:219829/149‑1 (MQ=255)
                                               aTCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTactga                                                                                            <  5:276347/148‑1 (MQ=255)
                                                 cGAGCTTAACGAAGCGTTTGCCAGCCAACGCCTGGCGGTGCTGCGTGAGCTGGGTCTGGCTGACGATGCGCCCCAGGTAAACCATAATGCCGGTGCCATTGCCCTGGGCCAcc                                                                                                                             >  5:157670/1‑113 (MQ=255)
                                                  gAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGc                                                                                        <  4:134377/149‑1 (MQ=255)
                                                      ttAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGtta                                                                                   >  4:64397/1‑148 (MQ=255)
                                                           gAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGa                                                                                                                 >  7:87668/1‑115 (MQ=255)
                                                           gAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGa                                                                                                                 <  8:87668/115‑1 (MQ=255)
                                                           gAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCAc                                                                                >  2:278357/1‑147 (MQ=255)
                                                           gAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCAc                                                                                >  2:88095/1‑147 (MQ=255)
                                                            aaGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAg                                                                             >  2:38150/1‑149 (MQ=255)
                                                            aaGCGTTTGCCAACCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAg                                                                             >  4:323691/1‑149 (MQ=255)
                                                             aGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGc                                                                            >  1:316113/1‑149 (MQ=255)
                                                                 tttGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAGTGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGga                                                                        <  3:319628/149‑1 (MQ=255)
                                                                 tttGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTg                                                                          >  4:262996/1‑147 (MQ=255)
                                                                 tttGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGga                                                                        <  1:351729/149‑1 (MQ=255)
                                                                     ccagccaAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGcc                                                                                                     <  8:311294/117‑1 (MQ=255)
                                                                     ccagccaAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGcc                                                                                                     >  7:311294/1‑117 (MQ=255)
                                                                          caAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTgg                                                               <  1:337138/149‑1 (MQ=255)
                                                                            aGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCg                                                             <  7:302663/149‑1 (MQ=255)
                                                                                cTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCg                                                         >  3:118531/1‑149 (MQ=255)
                                                                                 tGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCg                                                             <  3:12884/143‑1 (MQ=255)
                                                                                 tGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCg                                                             >  4:12884/1‑143 (MQ=255)
                                                                                   gCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCaa                                                      <  3:261481/149‑1 (MQ=255)
                                                                                   gCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCa                                                       <  7:312457/148‑1 (MQ=255)
                                                                                     ggTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTATGCACCAGCTGGAGAAGAGTGGCGGTCGCAAg                                                     <  1:38150/149‑1 (MQ=255)
                                                                                          tgcGTGAGCTGGGTGTGGCTGACAATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCGCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTg                                               <  2:102153/149‑1 (MQ=255)
                                                                                              tGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGagaaga                                                                   >  1:66506/1‑125 (MQ=255)
                                                                                              tGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGagaaga                                                                   <  2:66506/125‑1 (MQ=255)
                                                                                              tGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGa                                           <  6:20348/149‑1 (MQ=255)
                                                                                                aGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCCAGGCCCTTGCGAc                                          <  1:310208/148‑1 (MQ=255)
                                                                                                  cTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCAt                                       <  7:1924/149‑1 (MQ=255)
                                                                                                     ggTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCAtgt                                     <  3:262996/148‑1 (MQ=255)
                                                                                                      gtgtGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCAtgtgt                                   <  5:247483/149‑1 (MQ=255)
                                                                                                           gCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTgg                               <  1:88095/148‑1 (MQ=255)
                                                                                                             tGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCaa                                                      >  4:253548/1‑123 (MQ=255)
                                                                                                             tGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCaa                                                      <  3:253548/123‑1 (MQ=255)
                                                                                                                 gATGCGCCCCAGGTAAACCCTAATGGCGGTGCCGTTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTactgact                                                                                          <  6:319187/84‑1 (MQ=255)
                                                                                                                 gATGCGCCCCAGGTAAACCCTAATGGCGGTGCCGTTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTactgact                                                                                          <  6:313546/84‑1 (MQ=255)
                                                                                                                 gATGCGCCCCAGGTAAACCCTAATGGCGGTGCCGTTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTactgact                                                                                          >  5:319187/1‑84 (MQ=255)
                                                                                                                 gATGCGCCCCAGGTAAACCCTAATGGCGGTGCCGTTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTactgact                                                                                          >  5:313546/1‑84 (MQ=255)
                                                                                                                   tGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGcc                      <  1:278357/149‑1 (MQ=255)
                                                                                                                           aGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTg              <  2:316113/149‑1 (MQ=255)
                                                                                                                                  cccTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTgg       >  6:55338/1‑149 (MQ=255)
                                                                                                                                   ccTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGc      >  2:245452/1‑149 (MQ=255)
                                                                                                                                       aTGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATc  >  4:346748/1‑149 (MQ=255)
                                                                                                                                       aTGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAACGCCAGACCTTGGccgac  <  1:327404/149‑4 (MQ=39)
                                                                                                                                                  |                                                                                                                                          
GCAAACTGACCGAGCGTCTGGGGATAGCGGTAAGTGATTTCGACGTGATCGAGCTTAACGAAGCGTTTGCCAGCCAAGGCCTGGCGGTGCTGCGTGAGCTGGGTGTGGCTGACGATGCGCCCCAGGTAAACCCTAATGGCGGTGCCATTGCCCTGGGCCACCCCCTGGGCATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATC  >  NC_002947/1570097‑1570380

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