Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,793,940 A→T 9.4% I209N (ATT→AAT)  PP_3355 ← beta‑ketothiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,793,9400AT 9.4% 246.7 / 14.0 86I209N (ATT→AAT) PP_3355beta‑ketothiolase
Reads supporting (aligned to +/- strand):  ref base A (42/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 = 9.16e-01

TGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAG  >  NC_002947/3793792‑3794089
                                                                                                                                                     |                                                                                                                                                      
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                                                                                                                            >  2:1410811/1‑150 (MQ=255)
                          aaCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCt                                                                                                                            >  1:565762/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:146814/1‑150 (MQ=255)
                                 gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                                     <  2:2402597/150‑1 (MQ=255)
                                 gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                                     <  2:1000208/150‑1 (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                                                                                                                                                >  2:2345657/1‑81 (MQ=255)
                                                                           ggCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCATTGCATTGCcgc‑g                                                                                                                                                <  1:2345657/81‑1 (MQ=255)
                                                                              cAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑TTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCAc                                                                        <  2:565762/150‑1 (MQ=255)
                                                                                gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGa                                                                                                    >  2:193306/1‑120 (MQ=255)
                                                                                gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGa                                                                                                    <  1:193306/120‑1 (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                                                                                                                                       <  2:1581677/53‑1 (MQ=255)
                                                                                                                tCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCa                                                                                                                                       >  1:1581677/1‑53 (MQ=255)
                                                                                                                 ccTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAggg                                     <  1:1410811/150‑1 (MQ=255)
                                                                                                                  cTGGCTG‑ACCTCACTCACACCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTATACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGc                                    >  1:425981/1‑150 (MQ=255)
                                                                                                                   tGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGc                                                              >  1:1579135/1‑123 (MQ=255)
                                                                                                                   tGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGc                                                              <  2:1579135/123‑1 (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                           <  2:876477/150‑1 (MQ=255)
                                                                                                                            ctcactcaCGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTc                           <  1:1939404/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                    >  1:777457/1‑150 (MQ=255)
                                                                                                                                   acgccacgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt                    <  2:535689/150‑1 (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:639968/1‑150 (MQ=255)
                                                                                                                                        acgcGGCATTGCATTGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg               >  1:921481/1‑150 (MQ=255)
                                                                                                                                        acgcGGCATTGCATTGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg               >  2:1006998/1‑150 (MQ=255)
                                                                                                                                        acgcGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg               >  2:1886960/1‑150 (MQ=255)
                                                                                                                                         cgcgGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcgc              <  1:1870084/150‑1 (MQ=255)
                                                                                                                                              caatgcaatgCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAg         <  2:947447/150‑1 (MQ=255)
                                                                                                                                              caatgcaatgCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAg         >  1:1890420/1‑150 (MQ=255)
                                                                                                                                                atgcaatgcCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCATGCGCCAGGc       <  2:841248/150‑1 (MQ=255)
                                                                                                                                                    aatncCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATa   >  2:1505490/1‑150 (MQ=255)
                                                                                                                                                    aatgcCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATa   >  2:547061/1‑150 (MQ=255)
                                                                                                                                                     atgcCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAg  >  1:1910713/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTG‑ACCTCACTCACGCCACGCGGCAATGCAATGCCGC‑GTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAG  >  NC_002947/3793792‑3794089

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