Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,664,370 G→T 6.9% A473A (GCG→GCT treF → cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,664,3700GT 6.9% 1527.4 / 57.6 404A473A (GCG→GCTtreFcytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base G (162/214);  new base T (16/12);  total (178/226)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 408 total aligned reads displayed.

GTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGG  >  CP009273/3664221‑3664519
                                                                                                                                                     |                                                                                                                                                     
   ggaaTGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCTTGCGCAGACGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAAACGGTGAACAGTGGGATAAAACCAACGGCTGGGCACCGTTCCAAAGGAAGGAGAtt                                                                                                                                                    >  2:2322776/4‑150 (MQ=255)
   ggTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGAtt                                                                                                                                                    <  2:2440033/150‑1 (MQ=255)
   ggTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGAtt                                                                                                                                                    >  1:3008007/1‑150 (MQ=255)
     tATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCa                                                                                                                                                  >  2:3406563/1‑150 (MQ=255)
     tATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCa                                                                                                                                                  <  1:2344467/150‑1 (MQ=255)
      aTGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAg                                                                                                                                                 >  1:2824138/1‑150 (MQ=255)
        ggCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAggg                                                                                                                                               >  1:2513059/1‑150 (MQ=255)
         gCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGa                                                                                                                                              <  2:4603458/150‑1 (MQ=255)
          cGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGAt                                                                                                                                             >  2:1538030/1‑150 (MQ=255)
          cGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGAt                                                                                                                                             >  1:3098987/1‑150 (MQ=255)
          cGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGAt                                                                                                                                             >  1:2740716/1‑150 (MQ=255)
          cGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGAt                                                                                                                                             >  1:5549795/1‑150 (MQ=255)
            aaCCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGAttt                                                                                                                                           <  1:4243092/150‑1 (MQ=255)
                  gAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATg                                                                                                                                     >  2:3799139/1‑150 (MQ=255)
                      aGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACg                                                                                                                                 >  2:674736/1‑150 (MQ=255)
                       ggCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACgg                                                                                                                                >  1:5385090/1‑150 (MQ=255)
                       ggCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACgg                                                                                                                                <  1:1117388/150‑1 (MQ=255)
                         ccGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCg                                                                                                                              <  1:3588402/150‑1 (MQ=255)
                           gATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGAt                                                                                                                            >  1:3782333/1‑150 (MQ=255)
                           gATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGAt                                                                                                                            >  1:2159556/1‑150 (MQ=255)
                           gATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGAt                                                                                                                            >  2:1948329/1‑150 (MQ=255)
                              cGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAc                                                                                                                         >  1:855624/1‑150 (MQ=255)
                              cGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAc                                                                                                                         >  2:4328487/1‑150 (MQ=255)
                              cGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAc                                                                                                                         >  2:1647622/1‑150 (MQ=255)
                                tCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGAcat                                                                                                                       >  2:3018771/1‑148 (MQ=255)
                                       aaCGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGt                                                                                                                <  2:4537260/150‑1 (MQ=255)
                                         cGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGtga                                                                                                              <  1:2425100/150‑1 (MQ=255)
                                             gTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGaa                                                                                                          <  1:1199050/150‑1 (MQ=255)
                                             gTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGaa                                                                                                          >  1:2547062/1‑150 (MQ=255)
                                               gcgcAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAAt                                                                                                        <  1:4509386/150‑1 (MQ=255)
                                                cgcAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATc                                                                                                       >  2:2769325/1‑150 (MQ=255)
                                                     tCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATcgcgcg                                                                                                  >  1:605454/1‑150 (MQ=255)
                                                       ggTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAAACGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTACTGGGTGATGAAATAGCGCGaa                                                                                                >  2:5595715/1‑150 (MQ=255)
                                                           aCTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctg                                                                                            >  1:3743702/1‑150 (MQ=255)
                                                            cTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGACATCGCGCGAAgctgg                                                                                           >  1:5084885/1‑150 (MQ=255)
                                                            cTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctgg                                                                                           >  1:2845694/1‑150 (MQ=255)
                                                            cTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctgg                                                                                           >  1:1768233/1‑150 (MQ=255)
                                                              gACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctggct                                                                                         <  2:366190/150‑1 (MQ=255)
                                                              gACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctggct                                                                                         <  2:369054/150‑1 (MQ=255)
                                                               acacCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAgctggctg                                                                                        <  2:1591793/150‑1 (MQ=255)
                                                                  ccTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAg                                                                                     >  2:946752/1‑150 (MQ=255)
                                                                  ccTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAg                                                                                     <  1:322336/150‑1 (MQ=255)
                                                                  ccTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAg                                                                                     <  1:1652449/150‑1 (MQ=255)
                                                                   cTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGa                                                                                    >  1:1902624/1‑150 (MQ=255)
                                                                       cGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGt                                                                                <  1:4965462/150‑1 (MQ=255)
                                                                       cGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGt                                                                                >  2:2628642/1‑150 (MQ=255)
                                                                       cGGGAGTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGt                                                                                <  2:5569170/150‑1 (MQ=255)
                                                                        gggATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTg                                                                               <  1:3330467/150‑1 (MQ=255)
                                                                            ttCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATc                                                                           >  1:992439/1‑150 (MQ=255)
                                                                            ttCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATc                                                                           <  2:4127059/150‑1 (MQ=255)
                                                                             tCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCa                                                                          >  2:512297/1‑150 (MQ=255)
                                                                               tGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGt                                                                        <  2:2165300/150‑1 (MQ=255)
                                                                                    aGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTAt                                                                   >  1:2624168/1‑150 (MQ=255)
                                                                                    aGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTAt                                                                   <  1:2029157/150‑1 (MQ=255)
                                                                                      cGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCTGTTCTATCt                                                                 >  1:2965223/1‑150 (MQ=255)
                                                                                      cGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCt                                                                 >  1:4214659/1‑150 (MQ=255)
                                                                                       gAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTg                                                                <  2:4894745/150‑1 (MQ=255)
                                                                                       gAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTg                                                                <  2:853650/150‑1 (MQ=255)
                                                                                           aCGAAACCGGTGAACATTGGGATAACCCCAACGGCTGGGCACCGTTTCACTGGATGGCGATTCAGGGATTTAAAAGGTACGGAGATGACCTTCTGGGTGATGTAATCGCGCGAAGCTGGCTGAAGACGGTGAATCATTTGTATCTGGAAc                                                            <  2:1306459/150‑1 (MQ=255)
                                                                                            cGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACa                                                           <  2:4011409/150‑1 (MQ=255)
                                                                                            cGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACa                                                           <  2:628832/150‑1 (MQ=255)
                                                                                            cGAAACCGGTGAACAGTGGGATAAACCCAACGGATGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACa                                                           <  2:2082151/150‑1 (MQ=255)
                                                                                             gAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAg                                                          <  1:5147515/150‑1 (MQ=255)
                                                                                              aaaCCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGc                                                         <  1:432705/150‑1 (MQ=255)
                                                                                               aaCCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAACTCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGca                                                        >  1:4607796/1‑150 (MQ=255)
                                                                                               aaCCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGca                                                        <  1:750903/150‑1 (MQ=255)
                                                                                                  cGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACCATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGATGACGGTGAATCAGTGCGATCTGGAACAGCACaa                                                     >  2:5347287/1‑150 (MQ=255)
                                                                                                    gTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                   <  2:1239690/150‑1 (MQ=255)
                                                                                                    gTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAAc                                                   <  1:5043435/150‑1 (MQ=255)
                                                                                                             gggATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGaaa                                          <  2:2547062/150‑1 (MQ=255)
                                                                                                              ggATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGaaaa                                         <  1:1415560/150‑1 (MQ=255)
                                                                                                                aTAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAAt                                       >  2:2874590/1‑150 (MQ=255)
                                                                                                                 tAAACCCATCGGCTGGGCACCGTTACAATGGATGGCGAGTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATa                                      <  2:3848327/150‑1 (MQ=255)
                                                                                                                 tAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATa                                      <  2:3084672/150‑1 (MQ=255)
                                                                                                                 tAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATa                                      <  2:222643/150‑1 (MQ=255)
                                                                                                                   aaCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAcc                                    <  1:4883838/150‑1 (MQ=255)
                                                                                                                   aaCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAcc                                    >  1:4016596/1‑150 (MQ=255)
                                                                                                                   aaCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATAcc                                    >  2:2576622/1‑150 (MQ=255)
                                                                                                                    aCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGTTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCa                                   <  1:3797359/150‑1 (MQ=255)
                                                                                                                    aCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCa                                   <  1:2395506/150‑1 (MQ=255)
                                                                                                                    aCCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCa                                   <  1:1230858/150‑1 (MQ=255)
                                                                                                                          cGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGc                             >  1:1731395/1‑150 (MQ=255)
                                                                                                                           ggCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGcc                            >  2:723406/1‑150 (MQ=255)
                                                                                                                            gCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCg                           <  2:1198107/150‑1 (MQ=255)
                                                                                                                             cTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGa                          <  1:765300/150‑1 (MQ=255)
                                                                                                                               gggCACCGTTACAATGGATGGCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATg                        <  2:4730444/150‑1 (MQ=255)
                                                                                                                                  cACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGgtg                     >  1:3729678/1‑150 (MQ=255)
                                                                                                                                    ccGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCGTCTTGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGtt                   <  2:3576342/150‑1 (MQ=255)
                                                                                                                                     cGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTc                  <  2:2737144/150‑1 (MQ=255)
                                                                                                                                      gTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTcc                 >  1:271608/1‑150 (MQ=255)
                                                                                                                                          cAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCcgc             <  2:1472609/150‑1 (MQ=255)
                                                                                                                                          cAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCcgc             >  2:1903006/1‑150 (MQ=255)
                                                                                                                                           aaTGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACCAACTGATCGAAAAATACCATATAGCCGATGGTGTTCCCcgcg            >  2:5227737/1‑150 (MQ=255)
                                                                                                                                           aaTGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCcgcg            <  2:5385090/150‑1 (MQ=255)
                                                                                                                                               gatggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAgg        >  1:4485387/1‑150 (MQ=255)
                                                                                                                                                atggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGACggc       >  2:5067856/1‑150 (MQ=255)
                                                                                                                                                atggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggc       <  1:667951/150‑1 (MQ=255)
                                                                                                                                                  ggCTATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggcgg     >  1:1020087/1‑150 (MQ=255)
                                                                                                                                                  ggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggcgg     <  2:4309609/150‑1 (MQ=255)
                                                                                                                                                  ggCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAggcgg     >  2:4217892/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCACCGTTACAATGGATGGCGATTCAGGGATTTAAAATGTACGGCGATGACCTTCTGGGTGATGAAATCGCGCGAAGCTGGCTGAAGACGGTGAATCAGTTCTATCTGGAACAGCACAAACTGATCGAAAAATACCATATTGCCGATGGTGTTCCCCGCGAAGGCGGCGG  >  CP009273/3664221‑3664519

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: