Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,793,935 A→T 9.3% L211M (TTG→ATG)  PP_3355 ← beta‑ketothiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,793,9350AT 9.3% 248.4 / 14.6 86L211M (TTG→ATG) PP_3355beta‑ketothiolase
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base T (6/2);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

AAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGC  >  NC_002947/3793787‑3794084
                                                                                                                                                     |                                                                                                                                                      
aaGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAATCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑GCCTCACTCACGCCACGCGgcaat                                                                                                                                                       <  2:1186042/150‑1 (MQ=255)
  ggCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGgcattgc                                                                                                                                                     <  1:676416/150‑1 (MQ=255)
   gCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGgcattgca                                                                                                                                                    <  1:2064444/150‑1 (MQ=255)
     tGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCCGTCGGCCCAGTGCGGCCGCGTCGGGAGGCCTTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGgcaatgcaat                                                                                                                                                  <  1:886775/150‑1 (MQ=255)
     tGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGgcaatgcaat                                                                                                                                                  >  2:1885831/1‑150 (MQ=255)
        cgccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCATTGCATTGcc                                                                                                                                               >  2:2233773/1‑150 (MQ=255)
         gccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCcg                                                                                                                                              >  1:188189/1‑150 (MQ=255)
           cgcgCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCAGCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATCCcgc‑g                                                                                                                                           >  2:460736/1‑150 (MQ=255)
           cgcgCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCcgc‑g                                                                                                                                           <  1:1815017/150‑1 (MQ=255)
             cgcCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCATTGCATTGCCGC‑Gtg                                                                                                                                         <  1:2233773/150‑1 (MQ=255)
             cgcCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtg                                                                                                                                         <  1:1885831/150‑1 (MQ=255)
              gcCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtgg                                                                                                                                        >  1:1122522/1‑150 (MQ=255)
               ccATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtggc                                                                                                                                       <  2:188189/150‑1 (MQ=255)
                 aTCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtggctg                                                                                                                                     >  1:953806/1‑150 (MQ=255)
                 aTCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtgg                                                                                                                                        <  2:2243202/147‑1 (MQ=255)
                 aTCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtgg                                                                                                                                        >  1:2243202/1‑147 (MQ=255)
                  tCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑Gtggctgg                                                                                                                                    >  1:51736/1‑150 (MQ=255)
                    gACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCa                                                                                                                                  <  1:442740/150‑1 (MQ=255)
                        accGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCt                                                                                                                              >  1:287010/1‑150 (MQ=255)
                          cgcgCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTc                                                                                                                            >  2:2386747/1‑150 (MQ=255)
                          cgcgCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTc                                                                                                                            <  1:1299868/150‑1 (MQ=255)
                               aaCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCt                                                                                                                       >  1:565762/1‑150 (MQ=255)
                               aaCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCt                                                                                                                       >  2:1410811/1‑150 (MQ=255)
                               aaCTGTTGCCAGCCGTCTGAACGCCACCGGTGCAGGCGAGGGGCATTCGGGCCGGTGCGGACGCGTCGGTAGGAGGAAAAAAGGTCAGC‑GCCTGTACCGCACTCACGCCACGCGGCAATGCAATGCCGA‑GTGGCTGGCGGCCTGCTACCt                                                                                                                       >  1:2389230/1‑150 (MQ=255)
                                aCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTc                                                                                                                      <  2:2389230/150‑1 (MQ=255)
                                  tGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGACAGTACGTC‑CTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCACCCTACTACCTccc                                                                                                                    >  2:564921/1‑149 (MQ=255)
                                    ttGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGt                                                                                                                  <  1:1845158/150‑1 (MQ=255)
                                      gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                                <  2:2402597/150‑1 (MQ=255)
                                      gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                                <  2:1000208/150‑1 (MQ=255)
                                      gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                                >  2:146814/1‑150 (MQ=255)
                                            cGTCTGCACGCGATCGGGGTGCGCCGGGCGGAGTCGGGCGAGTGCGGCCGAGTCGGTAGGCCGTGGGTGGCTGTCCTCGCTG‑CCCTTCCTAACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGCCAGTCTTCTACCACctgcgcctgcgc                                                                                                          <  2:473708/150‑1 (MQ=255)
                                            cGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTgcgc                                                                                                          <  1:48828/150‑1 (MQ=255)
                                                  cACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGGCGGTAGGCCGTGGGGGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGAc                                                                                                    >  2:26046/1‑150 (MQ=255)
                                                     gCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCt                                                                                                 >  2:1527764/1‑150 (MQ=255)
                                                                 cACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtc                                                                                     <  2:122101/150‑1 (MQ=255)
                                                                 cACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGAGGACATG‑ACGCCATCAGGGGCAATGATCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtc                                                                                     >  1:876477/1‑150 (MQ=255)
                                                                 cACCGGGCGCAGTCGGGCCAGTGAGGCCGGGTCGGTGGGAGGTGCGTGATTATTGCTGCGC‑GCCTCACTCGCGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtc                                                                                     <  2:522090/150‑1 (MQ=255)
                                                                  aCCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcg                                                                                    >  1:522090/1‑150 (MQ=255)
                                                                    cGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtngtc                                                                                  <  2:1507794/150‑1 (MQ=255)
                                                                    cGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtc                                                                                  <  2:287010/150‑1 (MQ=255)
                                                                                ggCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCATTGCATTGCcgc‑g                                                                                                                                           <  1:2345657/81‑1 (MQ=255)
                                                                                ggCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCATTGCATTGCcgc‑g                                                                                                                                           >  2:2345657/1‑81 (MQ=255)
                                                                                   cAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑TTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCAc                                                                   <  2:565762/150‑1 (MQ=255)
                                                                                     gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGa                                                                                               <  1:193306/120‑1 (MQ=255)
                                                                                     gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGa                                                                                               >  2:193306/1‑120 (MQ=255)
                                                                                     gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGGCGATTTCCTCGTGGAACCAGCCCTGCCCCt                                                                 >  1:418988/1‑150 (MQ=255)
                                                                                       gCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTgg                                                               <  1:1132187/150‑1 (MQ=255)
                                                                                        cGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTggg                                                              <  2:1168006/150‑1 (MQ=255)
                                                                                               gTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAg                                                       >  2:362374/1‑150 (MQ=255)
                                                                                                tCGGTAGGCCGTGGGTGGCGGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAgg                                                      <  1:118867/150‑1 (MQ=255)
                                                                                                   tggacgcGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGcc                                                   >  2:1907030/7‑150 (MQ=255)
                                                                                                   gTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGcc                                                   >  2:1566230/1‑150 (MQ=255)
                                                                                                   gTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGAGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGcc                                                   >  1:719271/1‑150 (MQ=255)
                                                                                                      ggCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGc                                                <  1:1470037/150‑1 (MQ=255)
                                                                                                       gCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCt                                               >  1:817147/1‑150 (MQ=255)
                                                                                                        ccGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTg                                              >  1:2081922/1‑150 (MQ=255)
                                                                                                              gtggCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCt                                        >  1:535689/1‑150 (MQ=255)
                                                                                                               tggCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTTCAGGGCCCGCTGATTGCtt                                       >  2:137249/1‑149 (MQ=255)
                                                                                                                 gCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCt                                                                                                              <  1:1011665/77‑1 (MQ=255)
                                                                                                                 gCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCt                                                                                                              >  2:1011665/1‑77 (MQ=255)
                                                                                                                   tGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCGTGCGCGGTGACGCTGACGa                                                                                            >  2:1190470/1‑93 (MQ=255)
                                                                                                                   tGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGa                                                                                            <  1:1190470/93‑1 (MQ=255)
                                                                                                                   tGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCa                                   >  1:947447/1‑150 (MQ=255)
                                                                                                                     tCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCa                                                                                                                                  >  1:1581677/1‑53 (MQ=255)
                                                                                                                     tCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCa                                                                                                                                  <  2:1581677/53‑1 (MQ=255)
                                                                                                                      ccTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAggg                                <  1:1410811/150‑1 (MQ=255)
                                                                                                                       cTGGCTG‑ACCTCACTCACACCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTATACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGc                               >  1:425981/1‑150 (MQ=255)
                                                                                                                        tGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGc                                                         <  2:1579135/123‑1 (MQ=255)
                                                                                                                        tGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGc                                                         >  1:1579135/1‑123 (MQ=255)
                                                                                                                          gCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCa                            >  2:108473/1‑150 (MQ=255)
                                                                                                                          gCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCa                            >  2:781967/1‑150 (MQ=255)
                                                                                                                                ccTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCGGCGCAGGGCCCACGCGt                       >  2:1960579/1‑150 (MQ=255)
                                                                                                                                 ctcactcaCGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTc                      <  1:1939404/150‑1 (MQ=255)
                                                                                                                                 ctcactcaCGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTc                      <  2:876477/150‑1 (MQ=255)
                                                                                                                                  tcactcacGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTcc                     <  2:51736/150‑1 (MQ=255)
                                                                                                                                    actcacGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCAc                   >  1:2143951/1‑150 (MQ=255)
                                                                                                                                    actcacGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCAc                   >  1:1828707/1‑150 (MQ=255)
                                                                                                                                      tcacGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGACTGCGCGGTGACGCTGACGAGTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCg                 >  1:292012/1‑150 (MQ=255)
                                                                                                                                        acgccacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt               <  2:535689/150‑1 (MQ=255)
                                                                                                                                        acgccacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt               >  1:777457/1‑150 (MQ=255)
                                                                                                                                         cgccacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGttt              <  1:564921/150‑1 (MQ=255)
                                                                                                                                         caccacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGttt              <  2:425981/148‑1 (MQ=255)
                                                                                                                                          gccacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTc             >  1:1624258/1‑150 (MQ=255)
                                                                                                                                            cacgcGGCATTGCATTGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAc           >  1:176140/1‑150 (MQ=255)
                                                                                                                                            cacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAgg                                 >  1:1390848/1‑128 (MQ=255)
                                                                                                                                            cacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAgg                                 <  2:1390848/128‑1 (MQ=255)
                                                                                                                                             acgcGGCATTGCATTGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg          >  1:921481/1‑150 (MQ=255)
                                                                                                                                             acgcGGCATTGCATTGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg          >  2:1006998/1‑150 (MQ=255)
                                                                                                                                             acgcGGCATTGCATTGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg          >  1:639968/1‑150 (MQ=255)
                                                                                                                                             acgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg          >  2:1886960/1‑150 (MQ=255)
                                                                                                                                              cgcgGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcgc         <  1:1870084/150‑1 (MQ=255)
                                                                                                                                                   caatgcaatgCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAg    >  1:1890420/1‑150 (MQ=255)
                                                                                                                                                   caatgcaatgCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAg    <  2:947447/150‑1 (MQ=255)
                                                                                                                                                     atgcaatgcCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCATGCGCCAGGc  <  2:841248/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGC  >  NC_002947/3793787‑3794084

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