Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,823,835 A→C 7.6% W53G (TGG→GGG)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,823,8350AC 7.6% 1121.4 / 8.0 396W53G (TGG→GGG) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (243/107);  new base C (3/40);  total (246/150)
Fisher's exact test for biased strand distribution p-value = 7.92e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 412 total aligned reads displayed.

GATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATA  >  CP009273/1823686‑1823984
                                                                                                                                                     |                                                                                                                                                     
gATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCa                                                                                                                                                       >  1:4957558/1‑150 (MQ=255)
gATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCa                                                                                                                                                       <  2:5143733/150‑1 (MQ=255)
gATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCa                                                                                                                                                       >  2:3064108/1‑150 (MQ=255)
gATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCa                                                                                                                                                       <  2:1034437/150‑1 (MQ=255)
  tcntcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTGATTTCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgc                                                                                                                                                     >  1:3319643/1‑150 (MQ=255)
  tcatcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgc                                                                                                                                                     >  2:3863636/1‑150 (MQ=255)
   catcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcg                                                                                                                                                    >  2:2876378/1‑150 (MQ=255)
    atcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcgc                                                                                                                                                   >  2:159994/1‑150 (MQ=255)
    atcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcgc                                                                                                                                                   >  2:729218/1‑150 (MQ=255)
    atcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcgc                                                                                                                                                   >  2:4861384/1‑150 (MQ=255)
    atcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcgc                                                                                                                                                   >  1:33991/1‑150 (MQ=255)
      cGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCgg                                                                                                                                                 <  2:2707514/150‑1 (MQ=255)
          gTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTGCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAc                                                                                                                                             >  2:4771478/1‑150 (MQ=255)
           tCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcg                                                                                                                                            >  1:2784834/1‑150 (MQ=255)
            cACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgc                                                                                                                                           >  1:606468/1‑150 (MQ=255)
            cACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgc                                                                                                                                           >  2:3731797/1‑150 (MQ=255)
             aCTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGAATGACGTTCAGCAATTGAAAGCCGCGCCCAGCGCGGAAAcgcc                                                                                                                                          <  1:2059957/150‑1 (MQ=255)
                tCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCa                                                                                                                                       >  1:2393198/1‑150 (MQ=255)
                              ccAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACa                                                                                                                         >  2:3153909/1‑150 (MQ=255)
                                     ggCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGc                                                                                                                  >  2:253732/1‑150 (MQ=255)
                                     ggCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGc                                                                                                                  >  1:2839369/1‑150 (MQ=255)
                                     ggCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGc                                                                                                                  >  1:1365170/1‑150 (MQ=255)
                                        ttACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCg                                                                                                               >  1:628794/1‑150 (MQ=255)
                                          aCCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                                             >  1:5247398/1‑150 (MQ=255)
                                          aCCCGTTTCTCTGGAAATAATCGCGTTTATTTCGGCTTTTTTGCTTTCCAGCAGTGCGGAAACGCGTTCGACAACGGAATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                                             <  2:146616/150‑1 (MQ=255)
                                           cccGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAcc                                                                                                            >  1:1590178/1‑150 (MQ=255)
                                            ccGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCt                                                                                                           >  2:2002923/1‑150 (MQ=255)
                                             cGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTg                                                                                                          >  1:2991630/1‑150 (MQ=255)
                                             cGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTg                                                                                                          >  1:4869157/1‑150 (MQ=255)
                                                ttCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGc                                                                                                       >  1:1812162/1‑150 (MQ=255)
                                                 tctcTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcg                                                                                                      >  2:1740618/1‑150 (MQ=255)
                                                 tctcTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcg                                                                                                      >  1:145324/1‑150 (MQ=255)
                                                 tctcTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcg                                                                                                      >  1:3674468/1‑150 (MQ=255)
                                                      ggCAATAAGCGCGGTTAATTCGGCTTTATTGCTTGCCAGCTGTGCGGCAAAGCGTTCGACACCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatc                                                                                                 <  2:3257324/150‑1 (MQ=255)
                                                       gCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcg                                                                                                >  2:2181929/1‑150 (MQ=255)
                                                       caaataatCGCGTTTTATCCGGCTTTTTGGCTTCCCAGCGGTGCGCCACTGGGTTCGACACCGGCATGACTTTCACAAAAGGAACCCCGCCCCCCGCGCGCACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcg                                                                                                <  1:3739074/148‑1 (MQ=255)
                                                        cAATAATCGGGTTTAATTCGGTTTTATTGCTTTCCAGCTGTGCGGAAACGCGTTCGACAACGGCATGACGTTCAGCAAATGAACGCCGCGCCCAGCGCGTAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcgg                                                                                               <  1:4733554/150‑1 (MQ=255)
                                                        cAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcgg                                                                                               >  2:1171233/1‑150 (MQ=255)
                                                        cAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcgg                                                                                               >  1:2226936/1‑150 (MQ=255)
                                                         aataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggc                                                                                              >  2:3343190/1‑150 (MQ=255)
                                                         aataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggc                                                                                              >  2:4065472/1‑150 (MQ=255)
                                                         aataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggc                                                                                              >  2:1232150/1‑150 (MQ=255)
                                                         aataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggc                                                                                              >  1:16030/1‑150 (MQ=255)
                                                          ataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggca                                                                                             >  2:1245813/1‑150 (MQ=255)
                                                            aatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCggcatcggcatc                                                                                           >  1:4506285/1‑150 (MQ=255)
                                                               cgcgGTTAAGGCGGCTTTTTTGGTTTCCCGCAGGGCGGCAAAGCGTTCGACACCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                                                                        <  2:2782160/150‑1 (MQ=255)
                                                                   gTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcct                                                                                    <  2:4258682/150‑1 (MQ=255)
                                                                            gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAAGGGCATGACGTTCAGCAAAGGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCCTTGCCTTGCCATAac                                                                           <  2:2476283/150‑1 (MQ=255)
                                                                             cTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAaca                                                                          <  2:2614155/150‑1 (MQ=255)
                                                                             cTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAaca                                                                          >  2:2106576/1‑150 (MQ=255)
                                                                              tttATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAacac                                                                         >  2:1871733/1‑150 (MQ=255)
                                                                                 tttGCTTTCCAGCAGTGCGGCAAAGGGTTCGACAACGGAATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTc                                                                      <  1:2565205/149‑1 (MQ=255)
                                                                                  ttGCTTTCCCGCAGTGCGGAAATGCGTTCGAAAACGGCATGACGTTCAGAAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                     <  2:632304/150‑1 (MQ=255)
                                                                                  ttGCTTTCCAGCTGTGCGGCAACGCGTGCGACACCGGATTGACGTTCACCAACTGAACGCCGCGCCCCGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                     <  2:3948401/150‑1 (MQ=255)
                                                                                  ttGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                     >  2:1782214/1‑150 (MQ=255)
                                                                                  ttGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCg                                                                     >  2:4606105/1‑150 (MQ=255)
                                                                                    gCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGcc                                                                   >  1:3251696/1‑150 (MQ=255)
                                                                                     cTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGccc                                                                  >  1:425090/1‑150 (MQ=255)
                                                                                         cccgctgTGCGGCAAAGCTTTCGACAACGGCTTGACGTTCAGCAAATGAACGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                                                              <  2:606661/144‑1 (MQ=255)
                                                                                         ccAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                                                              >  1:2689906/1‑150 (MQ=255)
                                                                                          cagcagTGCGGCAACGGGTTCGACAAGGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc                                                             <  1:4794610/150‑1 (MQ=255)
                                                                                          cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc                                                             >  1:3173291/1‑150 (MQ=255)
                                                                                           agccgggCGGCAAAGCGTTCGACAACGGAATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAcc                                                            <  2:2252307/144‑1 (MQ=255)
                                                                                              agTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGa                                                         >  1:66011/1‑150 (MQ=255)
                                                                                                 ggggaaaaGGGTTCGAAATCGGATTTACGTTCAGAAAATGAATGCCGCGCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTa                                                      <  1:1167128/145‑1 (MQ=255)
                                                                                                 gCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATTAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTa                                                      <  1:1262865/150‑1 (MQ=255)
                                                                                                   ggCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACg                                                    >  1:1539251/1‑150 (MQ=255)
                                                                                                     cAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCt                                                  >  1:3974583/1‑150 (MQ=255)
                                                                                                     cAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCt                                                  >  2:3526030/1‑150 (MQ=255)
                                                                                                       aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc                                                >  1:2596490/1‑150 (MQ=255)
                                                                                                       aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc                                                >  2:4598614/1‑150 (MQ=255)
                                                                                                         gCGTTCGACACCGGCATGAGGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcac                                              <  2:2318557/150‑1 (MQ=255)
                                                                                                          ggtgcGACACGGGCTTGACGTGCATAAAATGAAAGCCGCGCCCCGCGGGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcaca                                             <  1:2911031/146‑1 (MQ=255)
                                                                                                            ttCGACACCGGCTTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACAcg                                           <  1:769551/150‑1 (MQ=255)
                                                                                                               gACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGa                                        >  2:5085224/1‑150 (MQ=255)
                                                                                                                 cAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATg                                      >  2:4443368/1‑150 (MQ=255)
                                                                                                                  aaCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATgc                                     >  2:4166886/1‑150 (MQ=255)
                                                                                                                        aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTg                               >  2:4980061/1‑150 (MQ=255)
                                                                                                                           aCGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGcc                            >  1:113861/1‑150 (MQ=255)
                                                                                                                            cGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGccc                           >  1:2122056/1‑150 (MQ=255)
                                                                                                                                cAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTa                       >  2:2021337/1‑150 (MQ=255)
                                                                                                                                cAGAAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTa                       <  1:4648077/150‑1 (MQ=255)
                                                                                                                                  gCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACTAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATc                     <  2:1816152/150‑1 (MQ=255)
                                                                                                                                    aaaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCa                   >  2:4501834/1‑150 (MQ=255)
                                                                                                                                    aaaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCa                   >  1:3040131/1‑150 (MQ=255)
                                                                                                                                     aaTGAAAGCCGCGCCCAGCGATGAAACTCCGCACGGGCTGCCCGACAAGCATGCTCGACCTGAGCGGCATCGGAATCATAGCCATGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGAGCCCAGGCCCGTTATCCAg                  >  2:2493866/1‑150 (MQ=255)
                                                                                                                                          aaGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGGTATCCAGTCAcc             >  1:823805/1‑150 (MQ=255)
                                                                                                                                           aGCCGCGCCCCGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCg            <  2:4401843/150‑1 (MQ=255)
                                                                                                                                           aGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCg            >  2:3441202/1‑150 (MQ=255)
                                                                                                                                              cgcgCCCCGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTa         <  1:2464151/150‑1 (MQ=255)
                                                                                                                                               gcgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa        <  2:842568/150‑1 (MQ=255)
                                                                                                                                               gcgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa        <  2:2404760/150‑1 (MQ=255)
                                                                                                                                               gcgcCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTaa        <  2:2788409/150‑1 (MQ=255)
                                                                                                                                                cgcCCAGCGCGGACACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAAt       <  1:5085224/150‑1 (MQ=255)
                                                                                                                                                  cccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATcc     >  1:357987/1‑150 (MQ=255)
                                                                                                                                                  cccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATcc     >  1:4400551/1‑150 (MQ=255)
                                                                                                                                                  cccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATcc     <  1:5161719/150‑1 (MQ=255)
                                                                                                                                                   ccAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCa    >  2:2805647/1‑150 (MQ=255)
                                                                                                                                                    cAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGCTGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCAt   <  2:2967311/150‑1 (MQ=255)
                                                                                                                                                     aGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATa  >  2:503925/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATA  >  CP009273/1823686‑1823984

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