Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 538,400 G→T 5.2% E353* (GAA→TAA)  glxK → glycerate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273538,4000GT 5.2% 830.7 / 16.7 211E353* (GAA→TAA) glxKglycerate kinase II
Reads supporting (aligned to +/- strand):  ref base G (82/118);  new base T (6/5);  total (88/123)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 218 total aligned reads displayed.

AAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGAT  >  CP009273/538251‑538549
                                                                                                                                                     |                                                                                                                                                     
   cgccgcTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAg                                                                                                                                                    >  2:1106617/1‑150 (MQ=255)
      cgcTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCTGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGc                                                                                                                                                 <  1:3105213/150‑1 (MQ=255)
            gtgtGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCa                                                                                                                                           <  2:2666277/150‑1 (MQ=255)
               tgGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCg                                                                                                                                        >  1:1299337/1‑150 (MQ=255)
                  cGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTg                                                                                                                                     <  2:3607894/150‑1 (MQ=255)
                  cGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTg                                                                                                                                     >  2:2510094/1‑150 (MQ=255)
                   gTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGa                                                                                                                                    <  1:2419064/150‑1 (MQ=255)
                    tCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGaa                                                                                                                                   <  2:1395898/150‑1 (MQ=255)
                        tGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCa                                                                                                                               >  2:2349623/1‑150 (MQ=255)
                           cGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAAtc                                                                                                                            >  1:755059/1‑150 (MQ=255)
                            gAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAAtct                                                                                                                           >  2:1248576/1‑150 (MQ=255)
                             aaGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAAtctc                                                                                                                          <  2:3749034/150‑1 (MQ=255)
                              aGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATctct                                                                                                                         >  1:4018471/1‑150 (MQ=255)
                               gCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCtt                                                                                                                        <  1:2839638/150‑1 (MQ=255)
                               gCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCtt                                                                                                                        <  1:186893/150‑1 (MQ=255)
                                      aaTGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAgc                                                                                                                 >  1:2647301/1‑150 (MQ=255)
                                            cctgTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGa                                                                                                           >  1:109464/4‑150 (MQ=255)
                                                tGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAata                                                                                                       >  2:2444954/1‑150 (MQ=255)
                                                 gATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAatat                                                                                                      >  2:3359194/1‑150 (MQ=255)
                                                  aTTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATAtt                                                                                                     <  1:1773944/150‑1 (MQ=255)
                                                   ttGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTg                                                                                                    <  2:1963364/150‑1 (MQ=255)
                                                     gggATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGcc                                                                                                  <  1:1075963/150‑1 (MQ=255)
                                                        aTTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTgc                                                                                               <  1:3286931/150‑1 (MQ=255)
                                                         ttGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTgcg                                                                                              >  2:636573/1‑150 (MQ=255)
                                                         ttGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTgcg                                                                                              >  2:636541/1‑150 (MQ=255)
                                                            cTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCa                                                                                           <  2:2399464/150‑1 (MQ=255)
                                                            cTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCa                                                                                           <  1:3671361/150‑1 (MQ=255)
                                                            cTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCa                                                                                           <  2:316215/150‑1 (MQ=255)
                                                            cTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCa                                                                                           <  2:2399741/150‑1 (MQ=255)
                                                              ggCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCAtt                                                                                         >  2:2814554/1‑150 (MQ=255)
                                                              ggCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCAtt                                                                                         <  2:1038578/150‑1 (MQ=255)
                                                                cGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTaa                                                                                       >  2:1987212/1‑150 (MQ=255)
                                                                 gTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTaaa                                                                                      >  2:180057/1‑150 (MQ=255)
                                                                 gTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTaaa                                                                                      >  2:73430/1‑150 (MQ=255)
                                                                    ttGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATa                                                                                   <  2:3100332/150‑1 (MQ=255)
                                                                    ttGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATa                                                                                   <  2:3250566/150‑1 (MQ=255)
                                                                    ttGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATa                                                                                   >  2:2262230/1‑150 (MQ=255)
                                                                     tGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAg                                                                                  <  1:3439728/150‑1 (MQ=255)
                                                                      gggTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAgg                                                                                 <  2:2210651/150‑1 (MQ=255)
                                                                       ggTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGt                                                                                <  2:36997/150‑1 (MQ=255)
                                                                        gTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTc                                                                               <  1:3019256/150‑1 (MQ=255)
                                                                        gTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTc                                                                               >  1:2224117/1‑150 (MQ=255)
                                                                        gTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTc                                                                               <  1:3034220/150‑1 (MQ=255)
                                                                        gTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTc                                                                               <  1:425723/150‑1 (MQ=255)
                                                                         tGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCa                                                                              <  1:347059/150‑1 (MQ=255)
                                                                           aTGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAgg                                                                            <  1:1202883/150‑1 (MQ=255)
                                                                           aTGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCAACTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAgg                                                                            <  1:3471595/150‑1 (MQ=255)
                                                                            tGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCTAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAggg                                                                           >  2:395206/1‑150 (MQ=255)
                                                                            tGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCTAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAggg                                                                           >  2:395158/1‑150 (MQ=255)
                                                                             ggCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGa                                                                          <  1:2565550/150‑1 (MQ=255)
                                                                              gCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGaa                                                                         <  1:687136/150‑1 (MQ=255)
                                                                              gCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGaa                                                                         <  2:4160630/150‑1 (MQ=255)
                                                                              gCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGaa                                                                         <  1:959119/150‑1 (MQ=255)
                                                                                gTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAAtt                                                                       <  2:1447008/150‑1 (MQ=255)
                                                                                   gAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTaaa                                                                    >  1:163553/1‑150 (MQ=255)
                                                                                     aGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTaaaaa                                                                  <  2:1008200/150‑1 (MQ=255)
                                                                                      gtggtgCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAAc                                                                 >  2:532636/1‑150 (MQ=255)
                                                                                       tggtgCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACt                                                                <  1:1064206/150‑1 (MQ=255)
                                                                                         gtgCACCAGGACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTaa                                                              <  2:3825952/150‑1 (MQ=255)
                                                                                          tgCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAAc                                                             <  1:2184970/150‑1 (MQ=255)
                                                                                              ccAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCtt                                                         <  1:1731810/150‑1 (MQ=255)
                                                                                              ccAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCtt                                                         <  1:1257025/150‑1 (MQ=255)
                                                                                               cAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTa                                                        >  2:641364/1‑150 (MQ=255)
                                                                                                aGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTAc                                                       <  2:779230/150‑1 (MQ=255)
                                                                                                   aCGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCtt                                                    >  2:1809860/1‑150 (MQ=255)
                                                                                                    cGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCttt                                                   <  2:1520998/150‑1 (MQ=255)
                                                                                                    cGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCttt                                                   <  1:1532397/150‑1 (MQ=255)
                                                                                                     ggCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTa                                                  <  2:3056922/150‑1 (MQ=255)
                                                                                                        aTTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAg                                               <  1:3057068/150‑1 (MQ=255)
                                                                                                        aTTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAg                                               >  1:4114428/1‑150 (MQ=255)
                                                                                                        aTTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAAACCTTACCTTTAAAg                                               >  1:3072332/1‑150 (MQ=255)
                                                                                                          tGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCg                                             <  2:3377405/150‑1 (MQ=255)
                                                                                                          tGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTGAAAGCg                                             >  1:202989/1‑150 (MQ=255)
                                                                                                                ggTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCg                                       <  1:841482/150‑1 (MQ=255)
                                                                                                                ggTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCg                                       >  2:4113638/1‑150 (MQ=255)
                                                                                                                ggTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCg                                       >  2:3218097/1‑150 (MQ=255)
                                                                                                                    ttCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGAttt                                   >  2:2656091/1‑150 (MQ=255)
                                                                                                                     tCAGCATTGTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTta                                  <  2:1157091/150‑1 (MQ=255)
                                                                                                                      cAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTtat                                 <  2:716905/150‑1 (MQ=255)
                                                                                                                       aGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTtata                                <  1:3903538/150‑1 (MQ=255)
                                                                                                                         cATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATc                              >  1:66151/1‑150 (MQ=255)
                                                                                                                         cATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATc                              >  1:1552191/1‑150 (MQ=255)
                                                                                                                          aTTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCg                             >  1:3399631/1‑150 (MQ=255)
                                                                                                                          aTTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCg                             >  1:2717632/1‑150 (MQ=255)
                                                                                                                             ttGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCAt                          <  2:3745966/150‑1 (MQ=255)
                                                                                                                               gCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATaa                        >  2:737502/1‑150 (MQ=255)
                                                                                                                               gCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATaa                        >  1:873407/1‑150 (MQ=255)
                                                                                                                                  tCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAagag                     <  2:1897855/150‑1 (MQ=255)
                                                                                                                                   cGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGt                    <  2:756891/150‑1 (MQ=255)
                                                                                                                                     tCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGc                  >  1:3728347/1‑150 (MQ=255)
                                                                                                                                             cTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCa          <  2:3913094/150‑1 (MQ=255)
                                                                                                                                              tttAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCAt         <  1:3798822/150‑1 (MQ=255)
                                                                                                                                              tttAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCAt         <  1:1441574/150‑1 (MQ=255)
                                                                                                                                                tAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGc       >  1:2643132/1‑150 (MQ=255)
                                                                                                                                                tAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGc       >  1:532632/1‑150 (MQ=255)
                                                                                                                                                tAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGc       <  2:1591940/150‑1 (MQ=255)
                                                                                                                                                 aGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGcc      <  2:2965769/150‑1 (MQ=255)
                                                                                                                                                 aGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGcc      >  1:3000897/1‑150 (MQ=255)
                                                                                                                                                 aGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGcc      <  1:3244093/150‑1 (MQ=255)
                                                                                                                                                   ccGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCgg    >  2:3877561/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGAT  >  CP009273/538251‑538549

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: ‑