Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,207,558 G→T 31.5% A174A (GCG→GCT aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,207,5580GT31.5% 580.9 / 266.5 312A174A (GCG→GCTaceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (97/116);  new base T (47/52);  total (144/168)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

AACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCTGACGTGACGGG  >  CP009273/4207409‑4207707
                                                                                                                                                     |                                                                                                                                                     
  caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgat                                                                                                                                                     >  2:2780534/1‑150 (MQ=255)
  caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgat                                                                                                                                                     <  1:22400/150‑1 (MQ=255)
   aacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCgatg                                                                                                                                                    >  2:2504199/1‑150 (MQ=255)
    acacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatga                                                                                                                                                   >  2:1382873/1‑150 (MQ=255)
     cacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatgat                                                                                                                                                  >  2:2647308/1‑150 (MQ=255)
     cacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCtatgat                                                                                                                                                  <  2:3320067/150‑1 (MQ=255)
      acCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGAtt                                                                                                                                                 <  1:210243/150‑1 (MQ=255)
       ccTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTg                                                                                                                                                <  2:3550004/150‑1 (MQ=255)
       ccTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTg                                                                                                                                                <  2:269249/150‑1 (MQ=255)
               cgtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGt                                                                                                                                        <  1:3765740/150‑1 (MQ=255)
                  gCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCa                                                                                                                                     <  2:1456361/150‑1 (MQ=255)
                   ccGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAg                                                                                                                                    <  2:2305145/150‑1 (MQ=255)
                   ccGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAg                                                                                                                                    <  1:744330/150‑1 (MQ=255)
                     gATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCg                                                                                                                                  <  1:1239225/150‑1 (MQ=255)
                      aTCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCgg                                                                                                                                 >  2:1328944/1‑150 (MQ=255)
                      aTCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCgg                                                                                                                                 >  1:905246/1‑150 (MQ=255)
                        cAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCa                                                                                                                               <  2:377247/150‑1 (MQ=255)
                        cAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCa                                                                                                                               >  2:1165569/1‑150 (MQ=255)
                         aGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAg                                                                                                                              <  1:580836/150‑1 (MQ=255)
                         aGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAg                                                                                                                              >  1:2831140/1‑150 (MQ=255)
                             ccAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCa                                                                                                                          >  2:204869/1‑150 (MQ=255)
                             ccAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCa                                                                                                                          >  2:253974/1‑150 (MQ=255)
                              cAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCAc                                                                                                                         >  1:4760642/1‑150 (MQ=255)
                               aaTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACt                                                                                                                        <  2:2224502/150‑1 (MQ=255)
                                  ggTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCg                                                                                                                     >  2:658681/1‑150 (MQ=255)
                                    tCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGaa                                                                                                                   <  2:3698283/150‑1 (MQ=255)
                                       gcgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGAt                                                                                                                <  2:3461654/150‑1 (MQ=255)
                                         gggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCa                                                                                                              >  1:1462504/1‑150 (MQ=255)
                                             aTTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                                          >  2:1291507/1‑150 (MQ=255)
                                             aTTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                                          >  1:3819274/1‑150 (MQ=255)
                                              ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTgg                                                                                                         <  1:640308/150‑1 (MQ=255)
                                               tGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGc                                                                                                        <  2:2426377/150‑1 (MQ=255)
                                                     gggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAga                                                                                                  <  1:662344/150‑2 (MQ=255)
                                                      ggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTg                                                                                                 <  2:1882936/150‑1 (MQ=255)
                                                       gCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTga                                                                                                <  1:3428187/150‑1 (MQ=255)
                                                        cGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaa                                                                                               >  1:2073830/1‑150 (MQ=255)
                                                         gATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaag                                                                                              >  2:2945195/1‑150 (MQ=255)
                                                          aTCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaaga                                                                                             >  1:1622958/1‑150 (MQ=255)
                                                          aTCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaaga                                                                                             >  2:1724851/1‑150 (MQ=255)
                                                               cgcTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGc                                                                                        <  2:1696951/150‑1 (MQ=255)
                                                                gcTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCg                                                                                       <  1:2603115/150‑1 (MQ=255)
                                                                 cTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCgg                                                                                      <  2:4706702/150‑1 (MQ=255)
                                                                  tATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGt                                                                                     >  1:1520107/1‑150 (MQ=255)
                                                                   aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTc                                                                                    >  2:2012626/1‑150 (MQ=255)
                                                                   aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTc                                                                                    >  2:2627680/1‑150 (MQ=255)
                                                                       cGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACAt                                                                                >  2:1016003/1‑150 (MQ=255)
                                                                           tACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggc                                                                            <  2:1540962/150‑1 (MQ=255)
                                                                           tACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggc                                                                            <  2:2446773/150‑1 (MQ=255)
                                                                             cTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggcgg                                                                          >  2:3410231/1‑150 (MQ=255)
                                                                             cTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggcgg                                                                          <  2:1062642/150‑1 (MQ=255)
                                                                             cTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggcgg                                                                          <  2:3847482/150‑1 (MQ=255)
                                                                               tCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCa                                                                        <  1:3698857/150‑1 (MQ=255)
                                                                               tCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCa                                                                        >  1:2984353/1‑150 (MQ=255)
                                                                                 cTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCaaa                                                                      <  2:3579857/150‑1 (MQ=255)
                                                                                   gCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGt                                                                    <  2:4447771/150‑1 (MQ=255)
                                                                                     cGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGttt                                                                  >  1:871600/1‑150 (MQ=255)
                                                                                      gATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                                 <  2:4668400/150‑1 (MQ=255)
                                                                                       aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                                >  1:2559414/1‑150 (MQ=255)
                                                                                       aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTa                                                                <  2:888354/150‑1 (MQ=255)
                                                                                         cGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGt                                                              <  1:517213/150‑1 (MQ=255)
                                                                                          gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                             >  1:4511233/1‑150 (MQ=255)
                                                                                          gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTg                                                             >  1:4062099/1‑150 (MQ=255)
                                                                                             gCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCa                                                          <  1:347571/150‑1 (MQ=255)
                                                                                               cGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                                        >  1:2010310/1‑150 (MQ=255)
                                                                                                gATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACt                                                       >  1:3964878/1‑150 (MQ=255)
                                                                                                 aTGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCGAACTc                                                      <  2:2608545/150‑1 (MQ=255)
                                                                                                    cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAgg                                                   <  1:3615025/150‑1 (MQ=255)
                                                                                                    cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAgg                                                   <  1:1851916/150‑1 (MQ=255)
                                                                                                     ggAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGa                                                  >  2:575649/1‑150 (MQ=255)
                                                                                                      gAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGaa                                                 <  2:3919512/150‑1 (MQ=255)
                                                                                                       aaGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAg                                                >  1:2801875/1‑150 (MQ=255)
                                                                                                          ccGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                             <  2:3653888/150‑1 (MQ=255)
                                                                                                          ccGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                             >  1:316627/1‑150 (MQ=255)
                                                                                                            ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAtt                                           >  1:1571515/1‑150 (MQ=255)
                                                                                                             gTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTc                                          >  1:1942409/1‑150 (MQ=255)
                                                                                                                 tGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaa                                      >  1:1013442/1‑150 (MQ=255)
                                                                                                                    cGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACt                                   <  1:4024856/150‑1 (MQ=255)
                                                                                                                    cGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACt                                   <  1:2945195/150‑1 (MQ=255)
                                                                                                                    cGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACt                                   >  1:303617/1‑150 (MQ=255)
                                                                                                                     ggTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTg                                  <  1:833217/150‑1 (MQ=255)
                                                                                                                        gtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTc                               <  2:4483832/150‑1 (MQ=255)
                                                                                                                          ccTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTcgc                             >  1:4732916/1‑150 (MQ=255)
                                                                                                                             gAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcggc                          <  2:506702/150‑1 (MQ=255)
                                                                                                                              aaTGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcggcg                         <  2:2412653/150‑1 (MQ=255)
                                                                                                                               aTGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc                        >  2:894958/1‑150 (MQ=255)
                                                                                                                               aTGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc                        <  2:358764/150‑1 (MQ=255)
                                                                                                                               aTGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc                        >  2:3102586/1‑150 (MQ=255)
                                                                                                                               aTGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTATTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc                        >  1:3348079/1‑150 (MQ=255)
                                                                                                                                    tttGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTg                   <  2:1233092/150‑1 (MQ=255)
                                                                                                                                      tGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGc                 >  1:1893862/1‑150 (MQ=255)
                                                                                                                                      tGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGc                 >  1:2115444/1‑150 (MQ=255)
                                                                                                                                       gAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCa                >  1:4620748/1‑150 (MQ=255)
                                                                                                                                              ttaaaGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacg         <  2:3043685/148‑1 (MQ=255)
                                                                                                                                              tgaAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacg         >  1:4352015/1‑150 (MQ=255)
                                                                                                                                               gaAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacgt        >  2:1633017/1‑150 (MQ=255)
                                                                                                                                                aaaGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacgtg       <  1:3295838/150‑1 (MQ=255)
                                                                                                                                                aaaGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacgtg       <  1:4416952/150‑1 (MQ=255)
                                                                                                                                                 aaGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGTGTCTGGCAGCtgacgtga      <  1:1056094/150‑1 (MQ=255)
                                                                                                                                                  aGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacgtgac     <  1:3085304/150‑1 (MQ=255)
                                                                                                                                                  aGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtgacgtgac     <  1:2647308/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCTGACGTGACGGG  >  CP009273/4207409‑4207707

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: