Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,866,685 T→G 37.4% F130V (TTT→GTT)  yeaJ → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,866,6850TG37.4% 171.2 / inf 337F130V (TTT→GTT) yeaJputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (95/116);  new base G (66/59);  total (162/175)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 343 total aligned reads displayed.

GCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGA  >  CP009273/1866537‑1866832
                                                                                                                                                    |                                                                                                                                                   
 cgcgcCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCtttt                                                                                                                                                   >  2:2443139/1‑150 (MQ=255)
   cgcCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGa                                                                                                                                                 <  1:4525745/150‑1 (MQ=255)
     ccATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATc                                                                                                                                               >  2:4003456/1‑150 (MQ=255)
         cTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGATTAAACCTGACTGCCCATAAATACCCGGGCTTACGAGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGAACTATTACCCGCTTTTGATCAGGc                                                                                                                                           <  2:1828031/150‑1 (MQ=255)
            cGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGt                                                                                                                                        >  1:2589260/1‑150 (MQ=255)
            cGAGGCTTTTGAGTCTGCCAATAACAACCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTTCGCGGCACACTACAAACCGCATACTCTGACTGCGACACAATTGTGGAAGCTGCAGCACGATTTCCCGCTGTTGATCAGGCAGt                                                                                                                                        <  1:666540/150‑1 (MQ=255)
             gAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTg                                                                                                                                       >  1:1167779/1‑150 (MQ=255)
                 cTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAg                                                                                                                                   >  1:3982817/1‑150 (MQ=255)
                       aGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAAcc                                                                                                                             <  1:844695/150‑1 (MQ=255)
                         tCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGc                                                                                                                           <  2:1473390/150‑1 (MQ=255)
                           tGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGcca                                                                                                                         <  2:2370146/150‑1 (MQ=255)
                               aaTAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAg                                                                                                                     <  2:4414999/150‑1 (MQ=255)
                                aTAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTAGTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAgg                                                                                                                    >  2:2062213/1‑150 (MQ=255)
                                  aaCACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGgat                                                                                                                  >  1:1402073/1‑150 (MQ=255)
                                     acCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGgatgat                                                                                                               <  2:3064577/150‑1 (MQ=255)
                                     acCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGgatgat                                                                                                               >  1:2976740/1‑150 (MQ=255)
                                       ccATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTa                                                                                                             >  1:3348073/1‑150 (MQ=255)
                                        cATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTAc                                                                                                            >  1:1736667/1‑150 (MQ=255)
                                        cATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTAc                                                                                                            <  2:3071100/150‑1 (MQ=255)
                                        cATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATACACTGACTGCGACACAAATGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTAc                                                                                                            >  2:3909722/1‑150 (MQ=255)
                                           ggCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTTCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGt                                                                                                         >  2:2522123/1‑150 (MQ=255)
                                            gCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGtt                                                                                                        >  2:1966281/1‑150 (MQ=255)
                                             cTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTc                                                                                                       <  1:934490/150‑1 (MQ=255)
                                                 aaCCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGt                                                                                                   >  1:212358/1‑150 (MQ=255)
                                                  aCCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTc                                                                                                  <  2:866040/150‑1 (MQ=255)
                                                      gACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTgg                                                                                              >  2:200921/1‑150 (MQ=255)
                                                      gACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTgg                                                                                              >  1:1350957/1‑150 (MQ=255)
                                                       aCTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTAGTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTggg                                                                                             <  2:3848689/150‑1 (MQ=255)
                                                         tGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGAt                                                                                           <  2:3327148/150‑1 (MQ=255)
                                                          gCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATg                                                                                          <  2:2611845/150‑1 (MQ=255)
                                                           cccATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGACGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATgg                                                                                         >  1:215857/1‑150 (MQ=255)
                                                             cATAAATACCCGGGCTTACGCGGCACACTACAAACCTCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCTCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGcc                                                                                       <  2:825327/150‑1 (MQ=255)
                                                             cATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGcc                                                                                       <  2:3170134/150‑1 (MQ=255)
                                                             cATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGcc                                                                                       >  1:2087250/1‑150 (MQ=255)
                                                             cATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGcc                                                                                       >  1:4152885/1‑150 (MQ=255)
                                                              aTAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCg                                                                                      <  1:2032789/150‑1 (MQ=255)
                                                              aTAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCg                                                                                      <  1:4360119/150‑1 (MQ=255)
                                                                aaaTACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                                                                    >  2:2535224/1‑150 (MQ=255)
                                                                     cccGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGaa                                                                               <  1:3312391/150‑1 (MQ=255)
                                                                      ccGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGaaa                                                                              <  1:541975/150‑1 (MQ=255)
                                                                       cGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAAt                                                                             >  1:360439/1‑150 (MQ=255)
                                                                        gggCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAAtt                                                                            >  2:1549726/1‑150 (MQ=255)
                                                                         ggCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAAttt                                                                           <  2:2931924/150‑1 (MQ=255)
                                                                            ttACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCAc                                                                        <  2:740313/150‑1 (MQ=255)
                                                                             tACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACt                                                                       >  2:608808/1‑150 (MQ=255)
                                                                             tACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACt                                                                       >  2:3751739/1‑150 (MQ=255)
                                                                              aCGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACta                                                                      <  1:1207617/150‑1 (MQ=255)
                                                                               cgcgGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACtat                                                                     >  2:470871/1‑150 (MQ=255)
                                                                                gcgGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACtatt                                                                    <  2:4061985/150‑1 (MQ=255)
                                                                                  ggCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACtattat                                                                  <  1:716019/150‑1 (MQ=255)
                                                                                    cacaCTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTAtct                                                                <  2:2772955/150‑1 (MQ=255)
                                                                                     acacTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTAtctc                                                               <  2:3155000/150‑1 (MQ=255)
                                                                                     acacTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTAtctc                                                               >  2:221993/1‑150 (MQ=255)
                                                                                      cacTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCg                                                              >  2:1770361/1‑150 (MQ=255)
                                                                                       acTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGa                                                             <  1:4133538/150‑1 (MQ=255)
                                                                                         tACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                           <  1:2614161/150‑1 (MQ=255)
                                                                                         tACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                           >  1:2883521/1‑150 (MQ=255)
                                                                                           cAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtctc                                                         >  1:3069383/1‑150 (MQ=255)
                                                                                            aaaCCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCa                                                        <  2:4239670/150‑1 (MQ=255)
                                                                                            aaaCCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCa                                                        <  1:4194810/150‑1 (MQ=255)
                                                                                            aaaCCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCa                                                        <  1:3842363/150‑1 (MQ=255)
                                                                                              aCCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTAGTATCTCGATCTCAAt                                                      <  2:417649/150‑1 (MQ=255)
                                                                                               ccGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATg                                                     <  1:1879818/150‑1 (MQ=255)
                                                                                                  cATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGAcc                                                  >  2:1163510/1‑150 (MQ=255)
                                                                                                  cATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGAcc                                                  <  1:1000422/150‑1 (MQ=255)
                                                                                                   aTCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCg                                                 >  1:3641829/1‑150 (MQ=255)
                                                                                                    tCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGc                                                <  1:161576/150‑1 (MQ=255)
                                                                                                       actgactgCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTAt                                             <  1:4595735/150‑1 (MQ=255)
                                                                                                        ctgactgCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATg                                            >  1:739765/1‑150 (MQ=255)
                                                                                                          gactgCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTc                                          >  2:633905/1‑150 (MQ=255)
                                                                                                             tgCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAt                                       <  2:2639100/150‑1 (MQ=255)
                                                                                                              gCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAtt                                      <  2:3421912/150‑1 (MQ=255)
                                                                                                               cGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAttt                                     >  1:3130582/1‑150 (MQ=255)
                                                                                                               cGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAttt                                     <  2:3232156/150‑1 (MQ=255)
                                                                                                               cGACACAATTGTGGAAGCTGCAGCACTAATACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGAATACGGATCAGGTCTTGGGATGGCCGAAGAGAAAATTCACTATTATCTCGATCTCAAAGACCGCTAAGTCAAtat                                     >  1:3495511/1‑148 (MQ=255)
                                                                                                                gACACAATTGTGGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAtttt                                    >  2:3532990/1‑150 (MQ=255)
                                                                                                                   acaATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTAt                                 <  2:4561306/150‑1 (MQ=255)
                                                                                                                       ttGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGc                             <  1:2681117/150‑1 (MQ=255)
                                                                                                                        tgtgGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTCTGTCTATTTTTATGAGcc                            >  2:4557237/1‑150 (MQ=255)
                                                                                                                        tgtgGAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGcc                            >  1:1528175/1‑150 (MQ=255)
                                                                                                                         gtgGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCg                           >  1:2425391/1‑150 (MQ=255)
                                                                                                                            gAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGtt                        <  2:4588643/150‑1 (MQ=255)
                                                                                                                            gAAGCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGtt                        >  1:4201815/1‑150 (MQ=255)
                                                                                                                             aaGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTa                       >  2:3861246/1‑150 (MQ=255)
                                                                                                                               gCTGCAGCACTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAAt                     >  1:2370103/1‑150 (MQ=255)
                                                                                                                                cTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATg                    <  1:3481349/150‑1 (MQ=255)
                                                                                                                                 tGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGt                   >  2:120342/1‑150 (MQ=255)
                                                                                                                                   cagcaCTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTg                 >  1:426199/1‑150 (MQ=255)
                                                                                                                                   cagcaCCATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTg                 >  2:442236/1‑150 (MQ=255)
                                                                                                                                     gcaCTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGaa               <  2:167883/150‑1 (MQ=255)
                                                                                                                                      caCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAAt              >  2:1923581/1‑150 (MQ=255)
                                                                                                                                      caCTATTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAAt              <  1:3198903/150‑1 (MQ=255)
                                                                                                                                       aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATa             >  2:2346790/1‑150 (MQ=255)
                                                                                                                                          aTTACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACtt          >  1:4533018/1‑150 (MQ=255)
                                                                                                                                           ttACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACttt         >  2:2318861/1‑150 (MQ=255)
                                                                                                                                           ttACCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACttt         <  2:1284058/150‑1 (MQ=255)
                                                                                                                                            tACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATGATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTg        <  2:2899678/150‑1 (MQ=255)
                                                                                                                                             aCCCGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACGAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGc       <  1:3236486/150‑1 (MQ=255)
                                                                                                                                               ccGCTGTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGa     >  1:1546485/1‑150 (MQ=255)
                                                                                                                                                cGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGAt    >  1:4175549/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAATTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGA  >  CP009273/1866537‑1866832

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: