Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,050,591 C→A 17.5% A393A (GCG→GCT serA ← D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,050,5910CA17.5% 621.6 / 92.1 213A393A (GCG→GCTserAD‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (83/92);  new base A (21/17);  total (104/109)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

CGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTG  >  CP009273/3050442‑3050740
                                                                                                                                                     |                                                                                                                                                     
 gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg                                                                                                                                                      <  1:743968/150‑1 (MQ=255)
  cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc                                                                                                                                                     <  2:4525199/150‑1 (MQ=255)
   aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCt                                                                                                                                                    >  2:837038/1‑150 (MQ=255)
    aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt                                                                                                                                                   <  2:1135154/150‑1 (MQ=255)
       tAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                                                                                                                                >  1:1444706/1‑150 (MQ=255)
        aGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTc                                                                                                                                               <  2:1368016/150‑1 (MQ=255)
         gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg                                                                                                                                              <  2:2698915/150‑1 (MQ=255)
          aaTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCgg                                                                                                                                             <  1:4078621/150‑1 (MQ=255)
              cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAAc                                                                                                                                         <  1:2920915/150‑1 (MQ=255)
                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt                                                                                                                                       <  2:1463617/150‑1 (MQ=255)
                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt                                                                                                                                       <  1:2829605/150‑1 (MQ=255)
                 aaCGCGCCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTc                                                                                                                                      <  1:2642548/150‑1 (MQ=255)
                  aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCt                                                                                                                                     <  1:1804781/150‑1 (MQ=255)
                    gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                   <  2:1347846/150‑1 (MQ=255)
                    gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                   <  2:2198411/150‑1 (MQ=255)
                    gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTtc                                                                                                                                   <  2:5009223/150‑1 (MQ=255)
                     cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcg                                                                                                                                  <  2:2319264/150‑1 (MQ=255)
                      ggCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgt                                                                                                                                 <  1:641098/150‑1 (MQ=255)
                      ggCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTtcgt                                                                                                                                 >  2:1572061/1‑150 (MQ=255)
                       gCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTtcgtc                                                                                                                                <  1:2529184/150‑1 (MQ=255)
                        cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtcg                                                                                                                               >  1:1738151/1‑150 (MQ=255)
                             gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtc                                                                                                                          >  1:3219148/1‑150 (MQ=255)
                             gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtc                                                                                                                          <  1:4208794/150‑1 (MQ=255)
                              tgtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTtca                                                                                                                         >  2:1448583/1‑150 (MQ=255)
                               gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaa                                                                                                                        >  2:484751/1‑150 (MQ=255)
                               gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaa                                                                                                                        <  1:2133861/150‑1 (MQ=255)
                                 ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaata                                                                                                                      <  1:1971091/150‑1 (MQ=255)
                                  gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatat                                                                                                                     >  2:3129264/1‑150 (MQ=255)
                                   agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatc                                                                                                                    >  1:2871899/1‑150 (MQ=255)
                                   agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatc                                                                                                                    >  1:4158593/1‑150 (MQ=255)
                                    gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGACTAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatca                                                                                                                   >  1:3617673/1‑150 (MQ=255)
                                    gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatca                                                                                                                   <  2:864579/150‑1 (MQ=255)
                                      aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaat                                                                                                                 <  1:4124231/150‑1 (MQ=255)
                                      aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaat                                                                                                                 <  1:2399179/150‑1 (MQ=255)
                                      aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaat                                                                                                                 >  2:1704571/1‑150 (MQ=255)
                                         gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGTAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAc                                                                                                              <  2:4545189/150‑1 (MQ=255)
                                           tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACca                                                                                                            <  2:3615738/150‑1 (MQ=255)
                                                aaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGTCGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                                                                                                       <  2:4204566/150‑1 (MQ=255)
                                                aaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATaa                                                                                                       <  2:2026061/150‑1 (MQ=255)
                                                    gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCa                                                                                                   <  1:4017010/150‑1 (MQ=255)
                                                     ggCAAGTCAGTGACCTGCCCGTTGATTTTCGGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCAt                                                                                                  <  1:5093238/150‑1 (MQ=255)
                                                       cAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCt                                                                                                <  1:1466308/150‑1 (MQ=255)
                                                           tCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGc                                                                                            <  2:2631830/150‑1 (MQ=255)
                                                            cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCg                                                                                           >  1:1830172/1‑150 (MQ=255)
                                                              gTGACCTGCCCGTTGGTTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGa                                                                                         >  2:1947132/1‑150 (MQ=255)
                                                               tGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGaa                                                                                        >  1:1606846/1‑150 (MQ=255)
                                                               tGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCGTCAATATCAATAACCACATAACCCATCTGGGCGGaa                                                                                        <  2:4029686/150‑1 (MQ=255)
                                                                gACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAg                                                                                       >  1:2469747/1‑150 (MQ=255)
                                                                   cTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGttt                                                                                    >  2:4770386/1‑150 (MQ=255)
                                                                      cccGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCa                                                                                 <  2:4729122/150‑1 (MQ=255)
                                                                       ccGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAg                                                                                <  1:3364269/150‑1 (MQ=255)
                                                                       ccGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAg                                                                                <  2:2898992/150‑1 (MQ=255)
                                                                       ccGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAg                                                                                <  2:5114085/150‑1 (MQ=255)
                                                                       ccGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAg                                                                                <  2:3370204/150‑1 (MQ=255)
                                                                       acGTTGATTTTCAGAGAAGGGGAATGAGTACAGCAGACGGGCGAGAATGGTACCCTGAAGAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATACCCCATCGGGGCGGAAGTTTGCAg                                                                                <  2:5114610/149‑1 (MQ=255)
                                                                        cGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGa                                                                               >  1:1809530/1‑150 (MQ=255)
                                                                            gATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATAtt                                                                           <  2:1252272/150‑1 (MQ=255)
                                                                               tttCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTgcg                                                                        <  1:647016/150‑1 (MQ=255)
                                                                                 tCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGcgcg                                                                      >  2:589565/1‑150 (MQ=255)
                                                                                  cAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcgg                                                                     <  2:3804788/150‑1 (MQ=255)
                                                                                  cAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcggc                                                                    >  2:3589363/1‑150 (MQ=255)
                                                                                  cAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCgcgg                                                                     >  1:941008/1‑150 (MQ=255)
                                                                                      gaAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGAt                                                                 >  1:3042976/1‑150 (MQ=255)
                                                                                      gaAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGAt                                                                 <  1:3670858/150‑1 (MQ=255)
                                                                                      gaAGGGGAATGAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGAt                                                                 <  2:1787937/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                           <  1:1718346/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                           <  1:964554/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                           <  1:2769236/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                           >  2:2553640/1‑150 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                           <  1:3144315/150‑1 (MQ=255)
                                                                                            gAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGAc                                                           >  2:3432131/1‑150 (MQ=255)
                                                                                               ttAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGcc                                                        >  2:2421372/1‑150 (MQ=255)
                                                                                                tAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGccc                                                       >  1:2889114/1‑150 (MQ=255)
                                                                                                 aGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATCTTGCGCGGCGATGTTGACGCCCt                                                      >  2:3183924/1‑150 (MQ=255)
                                                                                                       gcagACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCgg                                                <  1:3865603/150‑1 (MQ=255)
                                                                                                           aCGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGaa                                            >  2:134086/1‑150 (MQ=255)
                                                                                                                cgcgAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGAttt                                       <  1:4877877/150‑1 (MQ=255)
                                                                                                                   gAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGt                                    >  2:4392556/1‑150 (MQ=255)
                                                                                                                   gAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAAAAACCACATAAACCATCTGGGCGGAAGTTAGCAGATATACAGCGGCGATGTTGACTCCCTGCTCGGCGAAGAtgttgt                                    >  2:4932587/1‑150 (MQ=255)
                                                                                                                    aaTGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGtt                                   >  2:4645219/1‑150 (MQ=255)
                                                                                                                     aTGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTc                                  >  1:4404030/1‑150 (MQ=255)
                                                                                                                      tGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCa                                 >  1:889312/1‑150 (MQ=255)
                                                                                                                             cGGAATAGCTTTCATTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGt                          >  1:2188395/1‑150 (MQ=255)
                                                                                                                                aaTAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                       >  1:4667728/1‑150 (MQ=255)
                                                                                                                                aaTAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                       >  1:866246/1‑150 (MQ=255)
                                                                                                                                aaTAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                       >  1:3471090/1‑150 (MQ=255)
                                                                                                                                aaTAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAg                       >  2:2613554/1‑150 (MQ=255)
                                                                                                                                  tAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCa                     >  1:4821083/1‑150 (MQ=255)
                                                                                                                                   aGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCAc                    <  1:1306447/150‑1 (MQ=255)
                                                                                                                                       ttCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGccc                <  2:225679/150‑1 (MQ=255)
                                                                                                                                        tCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCg               >  2:4264391/1‑150 (MQ=255)
                                                                                                                                          aTTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGa             >  2:3715055/1‑150 (MQ=255)
                                                                                                                                          aTTGCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGa             >  1:4736325/1‑150 (MQ=255)
                                                                                                                                            tGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACg           >  1:4977581/1‑150 (MQ=255)
                                                                                                                                             gCCTGCAGAGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACgg          <  2:3789567/150‑1 (MQ=255)
                                                                                                                                              ccTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGt         >  2:2138056/1‑150 (MQ=255)
                                                                                                                                                tGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGttt       <  1:3196134/150‑1 (MQ=255)
                                                                                                                                                 gCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGtttt      >  1:131615/1‑150 (MQ=255)
                                                                                                                                                   aGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCg    <  1:1533363/150‑1 (MQ=255)
                                                                                                                                                     cgcTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTg  >  1:3034973/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTG  >  CP009273/3050442‑3050740

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: