Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,646 G→T 15.4% R385S (CGC→AGC)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,6460GT15.4% 821.9 / 96.2 263R385S (CGC→AGC) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (108/114);  new base T (23/18);  total (131/132)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

ATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGT  >  CP009273/2808497‑2808795
                                                                                                                                                     |                                                                                                                                                     
atcaCACACTTTCTGATACGCTTCGCCGCCGTTAATACTAAACAGCGTTTGCGCGACGCGTTTCAGAACGCGGAACAGATCTTTAACCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCACCCTCAGACCCGGTTTGCg                                                                                                                                                       >  1:1935576/1‑150 (MQ=255)
  cacacaCTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGAc                                                                                                                                                     >  1:3445509/1‑150 (MQ=255)
  cacacaCTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGAc                                                                                                                                                     >  1:3610825/1‑150 (MQ=255)
  cacacaCTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGAc                                                                                                                                                     >  1:3980695/1‑150 (MQ=255)
     acacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCtt                                                                                                                                                  <  1:1233466/150‑1 (MQ=255)
      cacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTc                                                                                                                                                 <  2:4331363/150‑1 (MQ=255)
      cacTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTc                                                                                                                                                 <  2:1099741/150‑1 (MQ=255)
       acTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCCGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCg                                                                                                                                                >  1:164529/1‑150 (MQ=255)
       acTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCCGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCg                                                                                                                                                >  1:470717/1‑150 (MQ=255)
         tttCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAg                                                                                                                                              >  2:1763166/1‑150 (MQ=255)
           tCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGa                                                                                                                                            <  1:533238/150‑1 (MQ=255)
              gATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCa                                                                                                                                         >  1:146586/1‑150 (MQ=255)
              gATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCa                                                                                                                                         <  2:297501/150‑1 (MQ=255)
               aTACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCAc                                                                                                                                        >  2:2489982/1‑150 (MQ=255)
               aTACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCAc                                                                                                                                        >  1:2950762/1‑150 (MQ=255)
                 aCGCGTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCAccc                                                                                                                                      <  2:3662417/150‑1 (MQ=255)
                  cGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCg                                                                                                                                     >  1:3444456/1‑150 (MQ=255)
                  cGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCg                                                                                                                                     >  2:1861979/1‑150 (MQ=255)
                    cTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGt                                                                                                                                   >  2:3315659/1‑150 (MQ=255)
                    cTTCGCAGCCGTGAATACTATCCAGCGGTTGCGCGACGCGTTTCAGATCGCGGTACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCCGCCGATTCCCAGCGTCGGACCCGGGGTGCGTCCGTCGGGGAGCGCCCggg                                                                                                                                   >  2:2491675/1‑149 (MQ=255)
                         ccgccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAg                                                                                                                              <  1:2983257/150‑1 (MQ=255)
                          cgccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGt                                                                                                                             <  2:2741357/150‑1 (MQ=255)
                            ccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTa                                                                                                                           <  2:1935576/150‑1 (MQ=255)
                            ccgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTa                                                                                                                           <  2:3570309/150‑1 (MQ=255)
                             cgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTaa                                                                                                                          <  1:2183853/150‑1 (MQ=255)
                             cgTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTaa                                                                                                                          >  1:3030654/1‑150 (MQ=255)
                                 aaTACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAAcgcg                                                                                                                      >  2:1079324/1‑150 (MQ=255)
                                 aaTACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAAcgcg                                                                                                                      >  2:4439533/1‑150 (MQ=255)
                                 aaTACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAAcgcg                                                                                                                      <  1:3925580/150‑1 (MQ=255)
                                    aCTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCgtgt                                                                                                                   >  2:4151061/1‑150 (MQ=255)
                                     cTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTa                                                                                                                  <  2:1414088/150‑1 (MQ=255)
                                       aTCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACa                                                                                                                <  1:2422739/150‑1 (MQ=255)
                                         ccAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAgg                                                                                                              >  2:596481/1‑150 (MQ=255)
                                         ccAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAgg                                                                                                              >  2:816533/1‑150 (MQ=255)
                                          cAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGc                                                                                                             <  2:2077936/150‑1 (MQ=255)
                                          cAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGc                                                                                                             <  2:1209811/150‑1 (MQ=255)
                                             cGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACCGGCAAg                                                                                                          >  2:3948501/1‑150 (MQ=255)
                                                ttGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTc                                                                                                       <  2:493962/150‑1 (MQ=255)
                                                 tGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCg                                                                                                      <  2:3179006/150‑1 (MQ=255)
                                                  gcgcGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGc                                                                                                     <  2:225931/150‑1 (MQ=255)
                                                      gACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACt                                                                                                 >  1:1122832/1‑150 (MQ=255)
                                                      gACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACt                                                                                                 >  1:3274998/1‑150 (MQ=255)
                                                      gACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACt                                                                                                 <  2:2088505/150‑1 (MQ=255)
                                                             ttCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAtt                                                                                          <  2:2672231/150‑1 (MQ=255)
                                                             ttCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAtt                                                                                          >  1:480971/1‑150 (MQ=255)
                                                              tCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCAttt                                                                                         <  2:2809525/150‑1 (MQ=255)
                                                               cAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTc                                                                                        >  1:826355/1‑150 (MQ=255)
                                                               cAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTc                                                                                        >  1:3465631/1‑150 (MQ=255)
                                                               cAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTc                                                                                        >  1:3166180/1‑150 (MQ=255)
                                                                        gAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATc                                                                               >  2:2435318/1‑150 (MQ=255)
                                                                            aGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagcc                                                                           >  2:2524966/1‑150 (MQ=255)
                                                                             gATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagcca                                                                          <  1:122900/150‑1 (MQ=255)
                                                                             gATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagcca                                                                          <  2:690060/150‑1 (MQ=255)
                                                                              aTCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATcagccag                                                                         <  2:1202143/150‑1 (MQ=255)
                                                                                  ttACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGca                                                                     >  2:1234996/1‑150 (MQ=255)
                                                                                      ccACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACAcc                                                                 >  2:25541/1‑150 (MQ=255)
                                                                                          cTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGAc                                                             <  1:4151061/150‑1 (MQ=255)
                                                                                          cTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGAc                                                             <  2:1148089/150‑1 (MQ=255)
                                                                                               gTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATaa                                                        >  2:1406123/1‑150 (MQ=255)
                                                                                                    gggAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAgg                                                   >  1:20761/1‑150 (MQ=255)
                                                                                                      gAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTc                                                 >  1:2325784/1‑150 (MQ=255)
                                                                                                       aaCTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCg                                                <  1:662339/150‑1 (MQ=255)
                                                                                                        aCTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGa                                               >  2:2475988/1‑150 (MQ=255)
                                                                                                        aCTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGa                                               <  1:3995403/150‑1 (MQ=255)
                                                                                                         cTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAg                                              >  1:943784/1‑150 (MQ=255)
                                                                                                           gtgCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGa                                            <  1:3315659/150‑1 (MQ=255)
                                                                                                           gtgCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGa                                            >  2:3837239/1‑150 (MQ=255)
                                                                                                            tgCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGaa                                           <  1:2855656/150‑1 (MQ=255)
                                                                                                               ggTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCg                                        >  2:583349/1‑150 (MQ=255)
                                                                                                               ggTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCg                                        >  2:583312/1‑150 (MQ=255)
                                                                                                                    cGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACct                                   <  2:3038076/150‑1 (MQ=255)
                                                                                                                    cGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACct                                   <  1:1778059/150‑1 (MQ=255)
                                                                                                                     gCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACctg                                  >  2:1452903/1‑150 (MQ=255)
                                                                                                                       aGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACctgct                                <  2:312790/150‑1 (MQ=255)
                                                                                                                          cGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGtt                             <  1:2728135/150‑1 (MQ=255)
                                                                                                                           gATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCNTGTACANGCAAGTTNGCTACTGCTCATTTCCGGTGCATCNGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTc                            <  1:253654/150‑1 (MQ=255)
                                                                                                                             tACCCAGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCAt                          >  1:2466171/1‑150 (MQ=255)
                                                                                                                             tACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCAt                          >  2:3817916/1‑150 (MQ=255)
                                                                                                                             tACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCAt                          >  2:1500207/1‑150 (MQ=255)
                                                                                                                               cccAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAACCAGGTCGAGGAATCGCACCTGCTGTTCATCt                        <  2:237597/150‑1 (MQ=255)
                                                                                                                               cccAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCt                        <  1:582529/150‑1 (MQ=255)
                                                                                                                               cccAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCt                        >  1:3102635/1‑150 (MQ=255)
                                                                                                                                 cAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTac                      >  1:1089723/1‑150 (MQ=255)
                                                                                                                                  aGCGTCAGACCCGGTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTaca                     >  1:2251797/1‑150 (MQ=255)
                                                                                                                                     gTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCa                  >  2:470589/1‑150 (MQ=255)
                                                                                                                                     gTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCa                  >  2:216033/1‑150 (MQ=255)
                                                                                                                                       cAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGGGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAAt                <  2:1930596/150‑1 (MQ=255)
                                                                                                                                          aCCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCgg             >  2:1686403/1‑150 (MQ=255)
                                                                                                                                          aCCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCgg             <  1:1268811/150‑1 (MQ=255)
                                                                                                                                          aCCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCgg             >  2:1687206/1‑150 (MQ=255)
                                                                                                                                           cccGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGc            >  1:694818/1‑150 (MQ=255)
                                                                                                                                           cccGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGc            <  1:1826231/150‑1 (MQ=255)
                                                                                                                                             cGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCga          <  2:4507068/150‑1 (MQ=255)
                                                                                                                                             cGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCga          >  1:4323156/1‑150 (MQ=255)
                                                                                                                                              ggTTTGCTACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCgag         <  1:2457411/150‑1 (MQ=255)
                                                                                                                                              ggTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCgag         <  2:1663278/150‑1 (MQ=255)
                                                                                                                                              ggTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCgag         >  2:905815/1‑150 (MQ=255)
                                                                                                                                                 ttGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAAc      <  2:2596237/150‑1 (MQ=255)
                                                                                                                                                    cGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACggg   <  2:3427238/150‑1 (MQ=255)
                                                                                                                                                    cGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACggg   >  2:1108918/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCACACACTTTCTGATACGCTTCGCCGCCGTTAATACTATCCAGCGTTTGCGCGACGCGTTTCAGATCGCGGAACAGATCTTTACCCACCTGCGGTAACGGGAACTGTGCGGTTTCGCAGCCGATACCCAGCGTCAGACCCGGTTTGCGACCTTCGAGGATCACCCGGTTCCAGTTAACGCGTGTACAGGCAAGTTCGCTACTGCTCATTTCCGGTGCATCAGCCAGCGCACACCAGACCATAAACAGGTCGAGGAATCGCACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGT  >  CP009273/2808497‑2808795

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