Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,824,100 T→A 5.2% D308V (GAT→GTT)  astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,824,1000TA 5.2% 799.5 / 22.7 228D308V (GAT→GTT) astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (110/106);  new base A (6/6);  total (116/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 229 total aligned reads displayed.

CCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAG  >  CP009273/1823951‑1824248
                                                                                                                                                     |                                                                                                                                                    
cccTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGAt                                                                                                                                                      >  1:2955646/1‑150 (MQ=255)
  cTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCg                                                                                                                                                    >  1:3840529/1‑150 (MQ=255)
   tGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAATTCAAACGCTCGGTTGCCGGATCgg                                                                                                                                                   >  1:799512/1‑150 (MQ=255)
   tGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCgg                                                                                                                                                   >  1:2795850/1‑150 (MQ=255)
     gCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGACCGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTa                                                                                                                                                 <  1:1075082/150‑1 (MQ=255)
     gCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGTAGTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTa                                                                                                                                                 >  2:2838733/1‑150 (MQ=255)
     gCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTa                                                                                                                                                 <  2:3617440/150‑1 (MQ=255)
     gCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTa                                                                                                                                                 >  1:2523892/1‑150 (MQ=255)
       ccGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACg                                                                                                                                               >  1:2224074/1‑150 (MQ=255)
        cGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACCGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACgc                                                                                                                                              <  1:3974243/150‑1 (MQ=255)
            aTCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTNCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACgcacca                                                                                                                                          >  1:4786633/1‑150 (MQ=255)
            aTCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACgcacca                                                                                                                                          <  2:1174895/150‑1 (MQ=255)
                aGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACgcaccagcac                                                                                                                                      <  2:849029/150‑1 (MQ=255)
                 gTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAGATCAAACGCTCGGTTGCCGGATCGGTACgcaccagcacc                                                                                                                                     >  1:3323084/1‑150 (MQ=255)
                  tCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACgcaccagcacca                                                                                                                                    <  1:636763/150‑1 (MQ=255)
                  tCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACgcaccagcacca                                                                                                                                    <  2:3353945/150‑1 (MQ=255)
                    aCCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCAcg                                                                                                                                  <  2:1331280/150‑1 (MQ=255)
                       gTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAAAGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCgg                                                                                                                               >  2:153334/1‑150 (MQ=255)
                         tAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGaa                                                                                                                             >  2:5018287/1‑150 (MQ=255)
                         tAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGaa                                                                                                                             >  2:3394129/1‑150 (MQ=255)
                          aaTCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGaaa                                                                                                                            >  2:481527/1‑150 (MQ=255)
                              cATACAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGt                                                                                                                        <  1:2234738/150‑1 (MQ=255)
                                tAACGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGa                                                                                                                      <  1:2389148/150‑1 (MQ=255)
                                tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGa                                                                                                                      >  2:100746/1‑150 (MQ=255)
                                tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGa                                                                                                                      <  1:4977351/150‑1 (MQ=255)
                                tAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGa                                                                                                                      >  2:479580/1‑150 (MQ=255)
                                 aaaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGAt                                                                                                                     >  2:2344121/1‑150 (MQ=255)
                                   aGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATaa                                                                                                                   <  2:4530846/150‑1 (MQ=255)
                                     tCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAAtt                                                                                                                 >  1:58985/1‑150 (MQ=255)
                                          cTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCAt                                                                                                            <  1:3155381/150‑1 (MQ=255)
                                           tGTTTTCTCCTCTGCGCACCGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCAtt                                                                                                           <  1:832711/150‑1 (MQ=255)
                                           tGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCCCGCGGAAATGGTGATAATTTTCAtt                                                                                                           >  2:2644866/1‑150 (MQ=255)
                                             ttttCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTgg                                                                                                         >  1:4347739/1‑150 (MQ=255)
                                              tttCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGc                                                                                                        <  2:1585354/150‑1 (MQ=255)
                                               ttCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATACTTTTCATTGGCg                                                                                                       >  2:2728594/1‑150 (MQ=255)
                                                tctcCTCTGCGCACAGGCGCACCATACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGa                                                                                                      >  2:2495384/1‑150 (MQ=255)
                                                    ctctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACcag                                                                                                  >  1:1770066/1‑150 (MQ=255)
                                                     tctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACcagg                                                                                                 <  2:1058757/150‑1 (MQ=255)
                                                       tGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACcaggca                                                                                               >  2:1933485/1‑150 (MQ=255)
                                                       tGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACcaggca                                                                                               >  1:748113/1‑150 (MQ=255)
                                                        gcgcACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACcaggcag                                                                                              <  1:5056670/150‑1 (MQ=255)
                                                          gcACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCaggcaggc                                                                                            >  2:2387869/1‑150 (MQ=255)
                                                           cacaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCt                                                                                           <  1:2584927/150‑1 (MQ=255)
                                                           cacaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCt                                                                                           <  1:4285643/150‑1 (MQ=255)
                                                           cacaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCt                                                                                           <  1:2688437/150‑1 (MQ=255)
                                                            acaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTg                                                                                          >  2:752091/1‑150 (MQ=255)
                                                            acaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTg                                                                                          >  2:3022908/1‑150 (MQ=255)
                                                            acaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGAACGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTg                                                                                          >  1:4241147/1‑150 (MQ=255)
                                                             caGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTgg                                                                                         >  2:3208120/1‑150 (MQ=255)
                                                                  gcACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAAt                                                                                    >  1:4202558/1‑150 (MQ=255)
                                                                  gcACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAAt                                                                                    >  1:4835590/1‑150 (MQ=255)
                                                                       aGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGccc                                                                               >  1:1295717/1‑150 (MQ=255)
                                                                        gACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCt                                                                              <  2:4179182/150‑1 (MQ=255)
                                                                        gACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCt                                                                              <  1:4235923/150‑1 (MQ=255)
                                                                        gACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCt                                                                              <  1:3994298/150‑1 (MQ=255)
                                                                         aCGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTg                                                                             <  1:4693319/150‑1 (MQ=255)
                                                                           gAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCAGGATCGGAACGAACCAGCACCACTCGTAAATGGTGATAATTTTCATTGGCTAACAGGCAGGCTGGGAAATCTCCCTgcg                                                                           >  1:5175259/1‑150 (MQ=255)
                                                                                cgATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGc                                                                      <  1:3651252/150‑1 (MQ=255)
                                                                                 gATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCt                                                                     >  1:325259/1‑150 (MQ=255)
                                                                                 gATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCt                                                                     >  1:1513067/1‑150 (MQ=255)
                                                                                 gATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCt                                                                     >  1:3908680/1‑150 (MQ=255)
                                                                                     cccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGccc                                                                 <  2:4899791/150‑1 (MQ=255)
                                                                                      ccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCt                                                                <  2:2536951/150‑1 (MQ=255)
                                                                                      ccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCt                                                                <  1:4932823/150‑1 (MQ=255)
                                                                                          cGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTg                                                            <  2:2438836/150‑1 (MQ=255)
                                                                                           gTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGc                                                           <  1:1411137/150‑1 (MQ=255)
                                                                                            tGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGcc                                                          <  2:3914519/150‑1 (MQ=255)
                                                                                              gCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCAc                                                        <  1:4061933/150‑1 (MQ=255)
                                                                                                atgtGAGGGCATCCAGTTGTGCGGCGGTTAAAATCACACGCTCGGTGGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCAGTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACtt                                                      <  2:2887410/147‑1 (MQ=255)
                                                                                                aTTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACtt                                                      >  2:856046/1‑150 (MQ=255)
                                                                                                  ttGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTcc                                                    >  2:3065474/1‑150 (MQ=255)
                                                                                                     aGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTccacc                                                 >  2:941788/1‑150 (MQ=255)
                                                                                                     aGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTccacc                                                 >  2:3540059/1‑150 (MQ=255)
                                                                                                          aTCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCg                                            >  2:4667277/1‑150 (MQ=255)
                                                                                                          aTCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGGGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCg                                            >  1:94339/1‑150 (MQ=255)
                                                                                                              aGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACtt                                        <  2:3283609/150‑1 (MQ=255)
                                                                                                                ttGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACtttt                                      <  2:1817700/150‑1 (MQ=255)
                                                                                                                 tgtgCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTa                                     >  2:4719931/1‑150 (MQ=255)
                                                                                                                    gcggcgGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACgg                                  <  1:4747062/150‑1 (MQ=255)
                                                                                                                        cggTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATgg                              <  2:4475736/150‑1 (MQ=255)
                                                                                                                          gTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGgcg                            >  1:4554649/1‑150 (MQ=255)
                                                                                                                            tAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGgcgcg                          <  1:1630534/150‑1 (MQ=255)
                                                                                                                            tAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGgcgcg                          <  1:1442832/150‑1 (MQ=255)
                                                                                                                             aaaaTCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGgcgcgc                         >  2:4370797/1‑150 (MQ=255)
                                                                                                                               aaTCAAACGCTCGGTTGCCGGAACGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAc                       <  1:374018/150‑1 (MQ=255)
                                                                                                                               aaTCAAACGCTCGGTTGCCGGAACGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAc                       <  1:2985461/150‑1 (MQ=255)
                                                                                                                                aTCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAcg                      >  1:4722996/1‑150 (MQ=255)
                                                                                                                                  cAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAcgcg                    <  2:4815706/150‑1 (MQ=255)
                                                                                                                                   aaaCGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCgg                   >  2:536230/1‑150 (MQ=255)
                                                                                                                                   aaaCGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCgg                   >  2:2066476/1‑150 (MQ=255)
                                                                                                                                    aaCGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGt                  <  2:1275589/150‑1 (MQ=255)
                                                                                                                                      cGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCg                <  2:3739517/150‑1 (MQ=255)
                                                                                                                                      cGCTCGGGTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCg                <  2:4621191/150‑1 (MQ=255)
                                                                                                                                            gTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTca          >  1:1228495/1‑150 (MQ=255)
                                                                                                                                              tGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTcaca        >  1:1948903/1‑150 (MQ=255)
                                                                                                                                              tGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTcaca        >  2:3623912/1‑150 (MQ=255)
                                                                                                                                               gCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCacac       >  2:700132/1‑150 (MQ=255)
                                                                                                                                                ccGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACt      <  2:4504084/150‑1 (MQ=255)
                                                                                                                                                 cGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTc     >  2:2195267/1‑150 (MQ=255)
                                                                                                                                                 cGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTc     >  2:2614618/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAG  >  CP009273/1823951‑1824248

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: