Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,195,635 G→T T82T (ACC→ACA U712_11000 ← Chorismate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,195,6350GT98.7% 246.1 / NA 79T82T (ACC→ACAU712_11000Chorismate mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (35/43);  minor base A (1/0);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTCAGTATTTTTTGTCAATGATAAATCGGGTCTTAATACAACGGCTTTTTCTAAATATACATGCCTGATCCGTTCCTGAGGAACATCTGTCTGGACGGTCATCATGACTCTGATGCACTTTTTCAGACCGCCGGTTACGTCCATTTCCTGCATACAGGTTACCGGCACATACTGCCATCCTGAAAGCTCACGGACAGCTTTTGCCGGGAAAACGGAGTGCAAATCAGGTGTCGCCGAAAGAAGCATTTGAACAACATCTTCCGGTT  >  CP006881/2195501‑2195768
                                                                                                                                      |                                                                                                                                     
aaTTCAGTATTTTTTGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTa                                                                                                                                   <  4:177380/139‑1 (MQ=255)
aaTTCAGTATTTTTTGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTa                                                                                                                                   <  4:196212/139‑1 (MQ=255)
aaTTCAGTATTTTTTGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTa                                                                                                                                   <  3:186780/139‑1 (MQ=255)
  ttCAGTATTTTTTGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACg                                                                                                                                 >  8:425989/1‑139 (MQ=255)
  ttCAGTATTTTTTGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACg                                                                                                                                 <  5:84797/139‑1 (MQ=255)
       tATTTTTTGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCAt                                                                                                                            >  3:184988/1‑139 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  4:12029/138‑1 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  8:3864/138‑1 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  8:175973/138‑1 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  5:346066/138‑1 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  6:165317/138‑1 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  6:429594/138‑1 (MQ=255)
           ttttGTCAATGATAAATCGGGCCTCAATACGACAACATTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  2:293798/138‑1 (MQ=255)
           ttttGGCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  1:145387/138‑1 (MQ=255)
                tCAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCa                                                                                                                    <  7:7245/138‑1 (MQ=255)
                 cAATGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATa                                                                                                                  >  3:131125/1‑139 (MQ=255)
                  aaTGATAAATCGTGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATa                                                                                                                  <  6:64691/138‑1 (MQ=255)
                  aaTGATAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATAc                                                                                                                 <  8:375637/139‑1 (MQ=255)
                    tGATATATCGGGCCTCAATACGACATCTTTATCTAAATATACCTGTCTGATCTGATCCAGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACat                                                                                                               <  8:175356/139‑2 (MQ=255)
                       tAAATCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGt                                                                                                             <  2:358976/138‑1 (MQ=255)
                        aaaTCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGGCCGCCTGTTACGTCCATTTCCTGCATACATGTTa                                                                                                           <  2:361199/139‑1 (MQ=255)
                         aaTCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTa                                                                                                           <  3:325250/138‑1 (MQ=255)
                         aaTCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAc                                                                                                          <  1:304816/139‑1 (MQ=255)
                           tCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        <  4:359204/139‑1 (MQ=255)
                           tCGGGCCTCAATACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        <  3:428782/139‑1 (MQ=255)
                               gcctcaaTACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        <  7:425989/130‑1 (MQ=255)
                                cctcaaTACGACAACTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGtaca                                                                                                   >  7:57/5‑139 (MQ=255)
                                           caactatttCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATcc                                                                                        >  8:250709/7‑139 (MQ=255)
                                           caacTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATcc                                                                                        >  4:56559/4‑139 (MQ=255)
                                           caacTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATcc                                                                                        >  5:50131/4‑139 (MQ=255)
                                              cTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        >  6:231032/1‑120 (MQ=255)
                                              cTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        <  5:231032/120‑1 (MQ=255)
                                              cTTTTTCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTg                                                                                      >  4:26857/1‑138 (MQ=255)
                                                 tttCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATa                                                                                                 >  3:394660/1‑124 (MQ=255)
                                                 tttCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATa                                                                                                 <  4:394660/124‑1 (MQ=255)
                                                 tttCTAAATATACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGaaa                                                                                   <  7:250709/138‑1 (MQ=255)
                                                        atatACATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTcgcga                                                                           <  3:26857/139‑2 (MQ=255)
                                                         tataCATGTCTGATCTGATCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTcgcgaa                                                                          <  4:184988/139‑3 (MQ=255)
                                                                   tgatctgatcCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        >  7:32960/1‑99 (MQ=255)
                                                                   tgatctgatcCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCg                                                                                                        <  8:32960/99‑1 (MQ=255)
                                                                   tgatctgatcCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                >  5:176726/1‑139 (MQ=255)
                                                                   tgatctgatcCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                >  3:251046/1‑139 (MQ=255)
                                                                   tgatctgatcCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                >  3:25913/1‑139 (MQ=255)
                                                                   tgatctgatcCTGAGGGACATCTGTCTGGACCATCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                >  3:136726/1‑139 (MQ=255)
                                                                     atctgatCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGcc                                                               <  3:56559/132‑1 (MQ=255)
                                                                     atctgatCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                <  5:6208/131‑1 (MQ=255)
                                                                     atctgatCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCg                                                              >  6:364124/7‑139 (MQ=255)
                                                                          atCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCCACCTGAAAGCTCTCGAACAGCTTTTGCCGGGaa                                                          >  7:208839/2‑138 (MQ=255)
                                                                          atCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                <  1:12903/131‑1 (MQ=255)
                                                                          atCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGaa                                                          >  3:374874/2‑138 (MQ=255)
                                                                          atCCTGAGGGACATCTGTCTCGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGc                                                                >  2:12903/2‑132 (MQ=255)
                                                                           tCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         <  4:240189/74‑1 (MQ=255)
                                                                           tCCTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTc                                                                                                                         >  3:240189/1‑74 (MQ=255)
                                                                            ccTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAcc                                                                                                         <  4:18380/89‑1 (MQ=255)
                                                                            ccTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAcc                                                                                                         >  3:18380/1‑89 (MQ=255)
                                                                             cTGAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAAca                                                      <  6:50131/139‑2 (MQ=255)
                                                                               gAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAc                                                                                                          >  1:5105/1‑85 (MQ=255)
                                                                               gAGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAc                                                                                                          <  2:5105/85‑1 (MQ=255)
                                                                                aGGGACATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAAcagcg                                                   <  4:251046/139‑5 (MQ=255)
                                                                                  ggaCATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTg                                                 <  4:374874/137‑1 (MQ=255)
                                                                                    aCATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGc                                                >  4:378098/1‑138 (MQ=255)
                                                                                     cATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAAca                                                      >  8:178537/1‑130 (MQ=255)
                                                                                     cATCTGTCTGGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAAca                                                      <  7:178537/131‑2 (MQ=255)
                                                                                          gtctgGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCGTGTTACGTCCATTTCCTGCATACATGTTACCGTTACATACTGCCATCCTGAAAGCACGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAg                                         >  2:136420/1‑139 (MQ=255)
                                                                                             tgGACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtg                                      >  6:133972/1‑139 (MQ=255)
                                                                                               gACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtgta                                    <  4:131125/139‑2 (MQ=255)
                                                                                               gACCGTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtgta                                    <  4:25913/139‑2 (MQ=255)
                                                                                                accgTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGgtgtag                                   >  3:68638/4‑137 (MQ=255)
                                                                                                   gTCATCATGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGcc                                 >  8:6153/1‑138 (MQ=255)
                                                                                                    tcatcaTGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAcc                                                                                                         <  4:237050/65‑1 (MQ=255)
                                                                                                    tcatcaTGACTCTTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTAcc                                                                                                         >  3:237050/1‑65 (MQ=255)
                                                                                                              tctTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTa                                                                                                           <  6:85924/53‑1 (MQ=255)
                                                                                                              tctTATGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTa                                                                                                           >  5:85924/1‑53 (MQ=255)
                                                                                                                ttaTGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCAttt                    >  1:160555/3‑138 (MQ=255)
                                                                                                                  aTGCACTTTTTCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGaa                 >  6:250637/1‑139 (MQ=255)
                                                                                                                        tttttCAGACCGCCTGTTACGTCCATTTCCTGCATACATTTTACAGGTACATACTGCCATCCTGAAAGCTCTCGACCAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTCGCCGACAGCAGCATTTGAACAACAt           >  5:264582/1‑139 (MQ=255)
                                                                                                                        tttttCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTATTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAacaaca            >  6:185776/1‑138 (MQ=255)
                                                                                                                          tttCAGACCGCCTGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCt         >  3:362523/1‑139 (MQ=255)
                                                                                                                          tttCAGACCGCCAGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGCACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCt         >  4:340457/1‑139 (MQ=255)
                                                                                                                                 ccgcctGTTACGTCGATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgtt  <  3:110466/139‑1 (MQ=255)
                                                                                                                                 ccgcctGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgtt  <  5:133972/139‑1 (MQ=255)
                                                                                                                                 ccgcctGTTACGTCCATTTCCTGCATACATGTTACCGGTACATACTGCCATCCTGAAAGCTCGCGAACAGCTTTTGCCGGGAAAACAGCGTGCAAATCAGGTGTAGCCGACAGAAGCATTTGAACAACATCTTCCGgtt  <  5:250637/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
AATTCAGTATTTTTTGTCAATGATAAATCGGGTCTTAATACAACGGCTTTTTCTAAATATACATGCCTGATCCGTTCCTGAGGAACATCTGTCTGGACGGTCATCATGACTCTGATGCACTTTTTCAGACCGCCGGTTACGTCCATTTCCTGCATACAGGTTACCGGCACATACTGCCATCCTGAAAGCTCACGGACAGCTTTTGCCGGGAAAACGGAGTGCAAATCAGGTGTCGCCGAAAGAAGCATTTGAACAACATCTTCCGGTT  >  CP006881/2195501‑2195768

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