Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,872,054 A→T 100% E28V (GAG→GTG)  yeaP → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,872,0540AT98.2% 900.1 / 3.2 280E28V (GAG→GTG) yeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (133/142);  total (136/144)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 280 total aligned reads displayed.

TGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTAT  >  CP009273/1871908‑1872200
                                                                                                                                                  |                                                                                                                                                  
tgTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGAt                                                                                                                                                 <  1:2480325/150‑1 (MQ=255)
 gTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATg                                                                                                                                                <  2:3114886/150‑1 (MQ=255)
   tttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCt                                                                                                                                              >  2:1443141/1‑150 (MQ=255)
   tttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCt                                                                                                                                              <  2:633147/150‑1 (MQ=255)
      aTCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATAATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGa                                                                                                                                           >  2:5124939/1‑150 (MQ=255)
       tCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGaa                                                                                                                                          >  1:142224/1‑150 (MQ=255)
        cATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGaaa                                                                                                                                         >  2:2054347/1‑150 (MQ=255)
             cTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTc                                                                                                                                    <  2:4846324/150‑1 (MQ=255)
                cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACt                                                                                                                                 >  1:1617553/1‑150 (MQ=255)
                 aaaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                                >  2:1766021/1‑150 (MQ=255)
                  aaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGa                                                                                                                               <  1:4116536/150‑1 (MQ=255)
                   aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGat                                                                                                                              <  2:3935997/150‑1 (MQ=255)
                   aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGat                                                                                                                              <  1:2324150/150‑1 (MQ=255)
                   aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGat                                                                                                                              <  1:2287890/150‑1 (MQ=255)
                    gATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGata                                                                                                                             <  2:775903/150‑1 (MQ=255)
                       tAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATgg                                                                                                                          <  1:3737546/150‑1 (MQ=255)
                           taatCAGCATTTCCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATc                                                                                                                      <  1:3906120/150‑1 (MQ=255)
                           taatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATc                                                                                                                      >  2:1021683/1‑150 (MQ=255)
                           taatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATc                                                                                                                      <  1:5099681/150‑1 (MQ=255)
                                aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacct                                                                                                                 <  2:1145219/150‑1 (MQ=255)
                                aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacct                                                                                                                 >  1:282034/1‑150 (MQ=255)
                                aGCATATACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGACCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATAAGGAATCAacct                                                                                                                 >  2:1262285/1‑150 (MQ=255)
                                   aTTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAACGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCTCTGATATGGAATCAacctacc                                                                                                              <  1:1089388/150‑1 (MQ=255)
                                   aTTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctacc                                                                                                              >  1:4228591/1‑150 (MQ=255)
                                     ttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTg                                                                                                            >  1:4653258/1‑150 (MQ=255)
                                      tACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGa                                                                                                           >  1:2567508/1‑150 (MQ=255)
                                      tACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGa                                                                                                           >  1:348094/1‑150 (MQ=255)
                                       aCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGAc                                                                                                          <  1:2713935/150‑1 (MQ=255)
                                        cccGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGAcc                                                                                                         >  1:4368775/1‑150 (MQ=255)
                                         ccgccgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCa                                                                                                        >  1:2418059/1‑150 (MQ=255)
                                             cgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGt                                                                                                    >  1:2109483/1‑150 (MQ=255)
                                              gTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTg                                                                                                   >  1:2168050/1‑150 (MQ=255)
                                                aTCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGa                                                                                                 >  1:442783/1‑150 (MQ=255)
                                                 tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGAt                                                                                                <  1:2982900/150‑1 (MQ=255)
                                                  ccTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATg                                                                                               >  2:1195492/1‑150 (MQ=255)
                                                  ccTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATg                                                                                               <  1:1723177/150‑1 (MQ=255)
                                                   cTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGt                                                                                              <  1:4173743/150‑1 (MQ=255)
                                                        gttgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAg                                                                                         >  1:758918/1‑150 (MQ=255)
                                                         ttgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgc                                                                                        <  1:5263284/150‑1 (MQ=255)
                                                         ttgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgc                                                                                        <  2:3007193/150‑1 (MQ=255)
                                                           gttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCGACCTACCTGACCAAAGTGGATGTCGAAgcgc                                                                                      >  1:1868683/1‑150 (MQ=255)
                                                           gttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgc                                                                                      >  2:1514860/1‑150 (MQ=255)
                                                           gttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgc                                                                                      >  2:3154493/1‑150 (MQ=255)
                                                           gttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgc                                                                                      <  1:1166984/150‑1 (MQ=255)
                                                            ttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcg                                                                                     >  2:1687344/1‑150 (MQ=255)
                                                            ttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcg                                                                                     >  2:414961/1‑150 (MQ=255)
                                                              ccGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGcc                                                                                   >  1:4994446/1‑150 (MQ=255)
                                                               cGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCt                                                                                  <  2:2066007/150‑1 (MQ=255)
                                                                    cagatcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgcagc                                                                             <  1:4551286/150‑1 (MQ=255)
                                                                       atcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCata                                                                          >  1:620263/1‑150 (MQ=255)
                                                                       atcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCata                                                                          <  1:1076662/150‑1 (MQ=255)
                                                                        tcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCatat                                                                         <  1:2434748/150‑1 (MQ=255)
                                                                        tcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCatat                                                                         <  2:433319/150‑1 (MQ=255)
                                                                            atTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATg                                                                     <  2:4495296/150‑1 (MQ=255)
                                                                             ttATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGt                                                                    <  1:1563048/150‑1 (MQ=255)
                                                                               aTCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGttt                                                                  <  2:3821894/150‑1 (MQ=255)
                                                                                      gcgTCTCGCAATCCATTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTa                                                           <  2:138906/150‑1 (MQ=255)
                                                                                       cgTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTaa                                                          >  2:5147236/1‑150 (MQ=255)
                                                                                        gTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                         >  2:959173/1‑150 (MQ=255)
                                                                                        gTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                         <  1:1668164/150‑1 (MQ=255)
                                                                                        gTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                         >  1:1097646/1‑150 (MQ=255)
                                                                                         tctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACa                                                        <  1:3102837/150‑1 (MQ=255)
                                                                                         tctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACa                                                        <  2:218865/150‑1 (MQ=255)
                                                                                          ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAg                                                       >  1:1919055/1‑150 (MQ=255)
                                                                                              cAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                   <  2:3411341/150‑1 (MQ=255)
                                                                                               aaTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGa                                                  >  2:672238/1‑150 (MQ=255)
                                                                                               aaTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGa                                                  <  2:78815/150‑1 (MQ=255)
                                                                                                aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGaa                                                 <  2:2975965/150‑1 (MQ=255)
                                                                                                   ccTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAAt                                              <  1:1713723/150‑1 (MQ=255)
                                                                                                     ttGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                            <  2:4653640/150‑1 (MQ=255)
                                                                                                      tGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTa                                           <  2:2345934/150‑1 (MQ=255)
                                                                                                       gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTAc                                          <  1:4892049/150‑1 (MQ=255)
                                                                                                         cAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                        >  1:4756758/1‑150 (MQ=255)
                                                                                                          aaaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATc                                       <  2:141794/150‑1 (MQ=255)
                                                                                                           aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATcc                                      <  1:1312284/150‑1 (MQ=255)
                                                                                                           aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATcc                                      <  1:3354588/150‑1 (MQ=255)
                                                                                                           aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATcc                                      >  2:1429953/1‑150 (MQ=255)
                                                                                                            aGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                     <  1:4455224/150‑1 (MQ=255)
                                                                                                             gAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCg                                    >  1:3030062/1‑150 (MQ=255)
                                                                                                                  gTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGaa                               <  1:1443141/150‑1 (MQ=255)
                                                                                                                   tCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAt                              <  1:865781/150‑1 (MQ=255)
                                                                                                                    cACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtt                             <  1:9262/150‑1 (MQ=255)
                                                                                                                     aCTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAttt                            <  2:481011/150‑1 (MQ=255)
                                                                                                                       tGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                          >  1:5033412/1‑150 (MQ=255)
                                                                                                                        ggAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAc                         <  2:2567508/150‑1 (MQ=255)
                                                                                                                         gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAcc                        >  2:5265450/1‑150 (MQ=255)
                                                                                                                         gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAcc                        >  2:2513966/1‑150 (MQ=255)
                                                                                                                          aaaGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCg                       >  1:2961205/1‑150 (MQ=255)
                                                                                                                            aGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGtc                     >  2:1873424/1‑150 (MQ=255)
                                                                                                                                  gTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTgg               >  2:3579059/1‑150 (MQ=255)
                                                                                                                                  gTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTgg               >  2:269092/1‑150 (MQ=255)
                                                                                                                                   tCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACATCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTggg              >  2:1665966/1‑150 (MQ=255)
                                                                                                                                   tCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTggg              <  1:3501709/150‑1 (MQ=255)
                                                                                                                                   tCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTggg              <  1:1392994/150‑1 (MQ=255)
                                                                                                                                      gACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTAAGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAtt           >  2:2806309/1‑150 (MQ=255)
                                                                                                                                       aCAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGACTTTTACCGTCTCGTGGGATTa          >  2:4140268/1‑150 (MQ=255)
                                                                                                                                            ttCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt     >  1:399263/1‑150 (MQ=255)
                                                                                                                                             tCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGtt    <  2:421150/150‑1 (MQ=255)
                                                                                                                                             tCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGtt    <  2:3591390/150‑1 (MQ=255)
                                                                                                                                              cTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTa   <  1:4081720/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTAT  >  CP009273/1871908‑1872200

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: