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,0540AT96.6% 733.8 / 11.7 232E28V (GAG→GTG) yeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base T (110/114);  total (113/119)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 235 total aligned reads displayed.

GTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTT  >  CP009273/1871909‑1872198
                                                                                                                                                 |                                                                                                                                                
  tttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATg                                                                                                                                              >  1:339025/1‑148 (MQ=255)
      tCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTgg                                                                                                                                          >  2:176754/1‑148 (MQ=255)
       cATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGa                                                                                                                                         >  1:859074/1‑148 (MQ=255)
       cATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGa                                                                                                                                         >  1:543110/1‑148 (MQ=255)
        aTGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGaa                                                                                                                                        <  1:237703/148‑1 (MQ=255)
          gCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAAt                                                                                                                                      <  2:2045228/148‑1 (MQ=255)
           ccTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATg                                                                                                                                     >  2:3445878/1‑148 (MQ=255)
           ccTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATg                                                                                                                                     <  2:2119411/148‑1 (MQ=255)
             tGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGt                                                                                                                                   <  2:3533347/148‑1 (MQ=255)
               cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCa                                                                                                                                 <  1:3543216/148‑1 (MQ=255)
               cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCa                                                                                                                                 <  1:3324608/148‑1 (MQ=255)
                aaaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCAc                                                                                                                                >  1:2930425/1‑148 (MQ=255)
                 aaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACt                                                                                                                               >  1:2216150/1‑148 (MQ=255)
                  aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                              >  1:3430323/1‑148 (MQ=255)
                  aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                              >  2:1729888/1‑148 (MQ=255)
                  aGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTg                                                                                                                              <  1:2501205/148‑1 (MQ=255)
                   gATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGa                                                                                                                             <  2:335619/148‑1 (MQ=255)
                    aTTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGat                                                                                                                            <  1:2593271/148‑1 (MQ=255)
                         ataatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGa                                                                                                                       >  1:3222651/1‑148 (MQ=255)
                          taatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGaa                                                                                                                      <  2:683653/148‑1 (MQ=255)
                           aatCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAAt                                                                                                                     <  1:1250551/148‑1 (MQ=255)
                            atCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATc                                                                                                                    >  1:916875/1‑148 (MQ=255)
                             tCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCa                                                                                                                   >  2:319176/1‑148 (MQ=255)
                              cAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCaa                                                                                                                  <  1:1525034/148‑1 (MQ=255)
                               aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAac                                                                                                                 <  2:3147201/148‑1 (MQ=255)
                               aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAac                                                                                                                 <  1:2888282/148‑1 (MQ=255)
                                 cATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacct                                                                                                               >  1:2466556/1‑148 (MQ=255)
                                 cATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacct                                                                                                               >  1:2679542/1‑148 (MQ=255)
                                   tttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctac                                                                                                             >  2:3339990/1‑148 (MQ=255)
                                    ttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAacctacc                                                                                                            >  2:516130/1‑148 (MQ=255)
                                      aCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAGCCTACCTg                                                                                                          <  2:1288073/148‑1 (MQ=255)
                                      aCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTg                                                                                                          <  1:2525003/148‑1 (MQ=255)
                                            cgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCaaa                                                                                                    >  2:1118927/1‑148 (MQ=255)
                                             gTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAg                                                                                                   <  2:482188/148‑1 (MQ=255)
                                              tATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGt                                                                                                  <  1:1289459/148‑1 (MQ=255)
                                               aTCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTg                                                                                                 >  2:2376556/1‑148 (MQ=255)
                                                tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                                                                <  2:3262059/148‑1 (MQ=255)
                                                tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                                                                >  1:1674811/1‑148 (MQ=255)
                                                 ccTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGa                                                                                               >  1:3055485/1‑148 (MQ=255)
                                                 ccTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGa                                                                                               >  2:1366858/1‑148 (MQ=255)
                                                 ccTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGa                                                                                               >  1:435056/1‑148 (MQ=255)
                                                    ggAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGt                                                                                            >  2:911585/1‑148 (MQ=255)
                                                     gAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTc                                                                                           >  2:1594765/1‑148 (MQ=255)
                                                      aGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCg                                                                                          >  1:72785/1‑148 (MQ=255)
                                                        ttgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGaa                                                                                        >  1:968876/1‑148 (MQ=255)
                                                        ttgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGaa                                                                                        >  1:941752/1‑148 (MQ=255)
                                                         tgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAg                                                                                       <  1:1308763/148‑1 (MQ=255)
                                                         tgttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAg                                                                                       <  1:596550/148‑1 (MQ=255)
                                                                tgtCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCt                                                                                >  1:1172579/1‑148 (MQ=255)
                                                                      atcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgcagca                                                                          <  2:933093/148‑1 (MQ=255)
                                                                        cagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCata                                                                        <  2:1155175/148‑1 (MQ=255)
                                                                          gatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCAtatt                                                                      <  2:273301/148‑1 (MQ=255)
                                                                           atTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCAtatta                                                                     >  2:1248316/1‑148 (MQ=255)
                                                                           atTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCAtatta                                                                     >  2:1249724/1‑148 (MQ=255)
                                                                                  cccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGc                                                              <  1:1596576/148‑1 (MQ=255)
                                                                                      cgTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGt                                                          >  2:526604/1‑148 (MQ=255)
                                                                                        tctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTaa                                                        <  1:3244100/148‑1 (MQ=255)
                                                                                         ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCTAAGCGCGCCTGCAGCAGATTATGTTTGCCCGTAAc                                                       >  2:2322272/1‑148 (MQ=255)
                                                                                         ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                       <  2:2355909/148‑1 (MQ=255)
                                                                                         ctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAAc                                                       >  1:242894/1‑148 (MQ=255)
                                                                                             cAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGcc                                                   <  2:296841/148‑1 (MQ=255)
                                                                                             cAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGca                                                   <  1:1366858/148‑2 (MQ=255)
                                                                                               aTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                                 >  1:2884034/1‑148 (MQ=255)
                                                                                                tCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGa                                                >  1:705879/1‑148 (MQ=255)
                                                                                                 cccTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGaa                                               <  2:72785/148‑1 (MQ=255)
                                                                                                   cTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGaaaa                                             <  1:1058555/148‑1 (MQ=255)
                                                                                                     tGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATg                                           <  2:1545792/148‑1 (MQ=255)
                                                                                                      gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                          <  2:1674811/148‑1 (MQ=255)
                                                                                                      gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                          >  2:596007/1‑148 (MQ=255)
                                                                                                      gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                          <  1:2041820/148‑1 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      <  2:69227/148‑1 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      >  2:717775/1‑148 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      >  2:2680754/1‑148 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      >  1:2407775/1‑148 (MQ=255)
                                                                                                          aaGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACAt                                      <  1:1482215/148‑1 (MQ=255)
                                                                                                           aGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATc                                     >  2:668375/1‑148 (MQ=255)
                                                                                                           aGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATc                                     <  1:1690966/148‑1 (MQ=255)
                                                                                                             aaCAGTCACTTGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATTGAATCAAACTACCTGAACAAACTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   <  2:2489259/148‑1 (MQ=255)
                                                                                                             aaCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   >  1:1902871/1‑148 (MQ=255)
                                                                                                             aaCAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATccc                                   >  1:1729885/1‑148 (MQ=255)
                                                                                                                 gTCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGgag                               >  2:2642161/1‑148 (MQ=255)
                                                                                                                  tCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGgaga                              <  2:2774902/148‑1 (MQ=255)
                                                                                                                  tCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGgaga                              <  1:1564633/148‑1 (MQ=255)
                                                                                                                  tCACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGgaga                              <  1:3277505/148‑1 (MQ=255)
                                                                                                                   cACTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGaa                             <  2:1920289/148‑1 (MQ=255)
                                                                                                                    aCTGGAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAt                            <  2:1124416/148‑1 (MQ=255)
                                                                                                                        gAAAGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAATCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                        >  1:2601054/1‑148 (MQ=255)
                                                                                                                           aGTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCg                     >  1:3165993/1‑148 (MQ=255)
                                                                                                                            gTCTGGTCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGt                    <  1:180128/148‑1 (MQ=255)
                                                                                                                                ggTCCGACAGCTTCTGGTGATGCTGGAAATGGTCAATGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCg                <  2:2034470/148‑1 (MQ=255)
                                                                                                                                  tCCGACAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTg              <  1:1345225/148‑1 (MQ=255)
                                                                                                                                   ccGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTgg             <  2:334350/148‑1 (MQ=255)
                                                                                                                                       cAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAtt         >  1:1224899/1‑148 (MQ=255)
                                                                                                                                       cAGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAtt         >  1:569220/1‑148 (MQ=255)
                                                                                                                                        aGCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTa        <  2:1410099/148‑1 (MQ=255)
                                                                                                                                         gCTTCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTAc       >  1:1019122/1‑148 (MQ=255)
                                                                                                                                           ttCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTc     >  2:3030221/1‑148 (MQ=255)
                                                                                                                                            tCTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCg    >  1:3328166/1‑148 (MQ=255)
                                                                                                                                             cTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt   <  2:246519/148‑1 (MQ=255)
                                                                                                                                             cTGGTGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAA‑TTTACCGTCTCGTGGGATTACTCGtt  <  2:2884034/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTT  >  CP009273/1871909‑1872198

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