Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,575,020 G→A 9.0% intergenic (+34/‑203) yhhY → / → yhhZ putative acetyltransferase/putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,575,0200GA 9.0% 1050.9 / 63.5 310intergenic (+34/‑203)yhhY/yhhZputative acetyltransferase/putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
Reads supporting (aligned to +/- strand):  ref base G (150/132);  new base A (12/16);  total (162/148)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

CCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTA  >  CP009273/3574872‑3575167
                                                                                                                                                    |                                                                                                                                                   
 ccGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGgtgt                                                                                                                                                    >  2:470259/1‑149 (MQ=255)
  cGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTc                                                                                                                                                   >  2:1753937/1‑149 (MQ=255)
    gCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCAt                                                                                                                                                 <  1:3142942/149‑1 (MQ=255)
    gCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGAGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCAt                                                                                                                                                 >  2:1833452/1‑149 (MQ=255)
        ttAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCaa                                                                                                                                             >  1:5136649/1‑149 (MQ=255)
         tAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCaaa                                                                                                                                            >  1:2268409/1‑149 (MQ=255)
         tAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCaaa                                                                                                                                            <  1:3257940/149‑1 (MQ=255)
          aaGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTTGAAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGAGAAGTAAGATAGTGCCCTTTATCTGAGATGGAAAAAGGGTGTCATTCaaaa                                                                                                                                           >  1:3628844/1‑149 (MQ=255)
          aaGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCaaaa                                                                                                                                           <  1:4946501/149‑1 (MQ=255)
          aaGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCaaaa                                                                                                                                           >  1:3925352/1‑149 (MQ=255)
          aaGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTATCATTCaaaa                                                                                                                                           >  2:866520/1‑149 (MQ=255)
           aGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAAt                                                                                                                                          <  1:1399927/149‑1 (MQ=255)
           aGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAAt                                                                                                                                          <  1:4203488/149‑1 (MQ=255)
            ggTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATc                                                                                                                                         >  2:2605249/1‑149 (MQ=255)
              tCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGa                                                                                                                                       >  1:2620090/1‑149 (MQ=255)
              tCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGa                                                                                                                                       >  2:636779/1‑149 (MQ=255)
                tataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACa                                                                                                                                     >  2:4844814/1‑149 (MQ=255)
                     aaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTAATTCAAAATCGACATAcct                                                                                                                                >  2:2447460/1‑149 (MQ=255)
                      aaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAcctt                                                                                                                               >  1:4409318/1‑149 (MQ=255)
                        aTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcc                                                                                                                             >  2:4647244/1‑149 (MQ=255)
                        aTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcc                                                                                                                             <  2:451257/149‑1 (MQ=255)
                         tACGGATTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcct                                                                                                                            <  1:1440938/149‑1 (MQ=255)
                          aCGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcctt                                                                                                                           <  2:4219986/149‑1 (MQ=255)
                            ggCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTa                                                                                                                         >  1:1479410/1‑149 (MQ=255)
                             gCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTaa                                                                                                                        >  1:1040120/1‑149 (MQ=255)
                             gCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTaa                                                                                                                        >  1:2827563/1‑149 (MQ=255)
                             gCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTaa                                                                                                                        >  2:2946447/1‑149 (MQ=255)
                               tttGAAATTGAATGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTATCATTCAAAATCGACATACCTTCCTTTAAgg                                                                                                                      <  2:3977952/149‑1 (MQ=255)
                               tttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAgg                                                                                                                      >  2:1560356/1‑149 (MQ=255)
                                  gAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtat                                                                                                                   <  1:3545295/149‑1 (MQ=255)
                                    aaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtattt                                                                                                                 <  2:2829424/149‑1 (MQ=255)
                                     aTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtattta                                                                                                                <  1:2331203/149‑1 (MQ=255)
                                     aTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTATCATTCAAAATCGACATACCTTCCTTTAAGGtattta                                                                                                                >  1:4807448/1‑149 (MQ=255)
                                          aGGGACTGGTAAGAAGTACGCATTGCGTATTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTATCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGc                                                                                                           <  2:2330660/149‑1 (MQ=255)
                                          aGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTATCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGc                                                                                                           <  2:2272392/149‑1 (MQ=255)
                                           gggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGcc                                                                                                          >  1:4239635/1‑149 (MQ=255)
                                           gggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGcc                                                                                                          <  1:4488399/149‑1 (MQ=255)
                                            ggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGccc                                                                                                         <  2:3423803/149‑1 (MQ=255)
                                             gACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                        >  1:2113973/1‑149 (MQ=255)
                                               cTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAAt                                                                                                      >  2:731683/1‑149 (MQ=255)
                                                 ggTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATAc                                                                                                    >  1:503046/1‑149 (MQ=255)
                                                  gTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACa                                                                                                   >  2:2663253/1‑149 (MQ=255)
                                                      gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                               <  1:1753449/149‑1 (MQ=255)
                                                      gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                               <  1:4122299/149‑1 (MQ=255)
                                                      gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                               <  2:1461221/149‑1 (MQ=255)
                                                          tACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTga                                                                                           <  1:369411/149‑1 (MQ=255)
                                                          tACGCATTGCGTAATGGTGAATAGGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTga                                                                                           <  1:2925024/149‑1 (MQ=255)
                                                                 tGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTaa                                                                                    <  2:2113370/149‑1 (MQ=255)
                                                                    gTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTc                                                                                 >  2:2517274/1‑149 (MQ=255)
                                                                    gTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTc                                                                                 >  1:2630918/1‑149 (MQ=255)
                                                                     tAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCa                                                                                >  1:4180197/1‑149 (MQ=255)
                                                                      aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAt                                                                               >  2:263796/1‑149 (MQ=255)
                                                                      aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAt                                                                               >  1:367511/1‑149 (MQ=255)
                                                                       aTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATa                                                                              >  2:3804805/1‑149 (MQ=255)
                                                                        tGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                             <  2:3178773/149‑1 (MQ=255)
                                                                        tGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                             <  1:285235/149‑1 (MQ=255)
                                                                         ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAt                                                                            <  1:2605249/149‑1 (MQ=255)
                                                                          gTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                           >  1:3343925/1‑149 (MQ=255)
                                                                           tGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTa                                                                          >  2:4903195/1‑149 (MQ=255)
                                                                              atatGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACtt                                                                       <  1:801918/149‑1 (MQ=255)
                                                                              atatGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACtt                                                                       <  2:1802851/149‑1 (MQ=255)
                                                                                atGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTcc                                                                     >  1:2129879/1‑149 (MQ=255)
                                                                                atGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTcc                                                                     <  1:4873831/149‑1 (MQ=255)
                                                                                  gTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTa                                                                   >  1:40433/1‑149 (MQ=255)
                                                                                      aTGCATATTATATGGCGCGGGTGAAGTAAGATTGTGCCCTTTGTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                               <  2:2719627/149‑1 (MQ=255)
                                                                                      aTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                               >  1:2730650/1‑149 (MQ=255)
                                                                                          atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                           >  2:690085/1‑149 (MQ=255)
                                                                                           tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                          >  2:2187797/1‑149 (MQ=255)
                                                                                             ttatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAatat                                                        >  2:4382913/1‑149 (MQ=255)
                                                                                             ttatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAatat                                                        <  1:636779/149‑1 (MQ=255)
                                                                                             ttatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAatat                                                        >  2:2378772/1‑149 (MQ=255)
                                                                                                   ggCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttatt                                                  >  1:2187661/1‑149 (MQ=255)
                                                                                                    gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                 <  1:3069726/149‑1 (MQ=255)
                                                                                                     cgcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                <  2:3343925/149‑1 (MQ=255)
                                                                                                       cgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTAtttttt                                              <  2:4150639/149‑1 (MQ=255)
                                                                                                        gggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTa                                             >  2:2880572/1‑149 (MQ=255)
                                                                                                           tGAAGTAAGATAGTGCCCTTTGTCTGAGATTGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAata                                          <  2:3275841/149‑1 (MQ=255)
                                                                                                            gAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatat                                         >  1:1325511/1‑149 (MQ=255)
                                                                                                               gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTGAAGTGTAATATTTTATTTTTTAAtatata                                      <  2:3395750/149‑1 (MQ=255)
                                                                                                               gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAAtatata                                      >  2:529939/1‑149 (MQ=255)
                                                                                                                    aTAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCt                                 <  1:892828/149‑1 (MQ=255)
                                                                                                                      aGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTAc                               <  1:342405/149‑1 (MQ=255)
                                                                                                                      aGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTAc                               <  1:2663253/149‑1 (MQ=255)
                                                                                                                          cccTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAtt                           >  2:3033162/1‑149 (MQ=255)
                                                                                                                            ctatttCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTc                         <  2:1191145/146‑1 (MQ=255)
                                                                                                                            cTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAAACGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTc                         >  1:4651058/1‑149 (MQ=255)
                                                                                                                             tttttCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTcc                        >  1:3281006/1‑149 (MQ=255)
                                                                                                                                 tCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                    >  1:1957952/1‑149 (MQ=255)
                                                                                                                                      gaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaa               >  1:724983/1‑149 (MQ=255)
                                                                                                                                      gaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaa               <  1:2485777/149‑1 (MQ=255)
                                                                                                                                       aTGGAAAAAAGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaat              <  1:1215817/149‑1 (MQ=255)
                                                                                                                                        tGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaata             >  2:44270/1‑149 (MQ=255)
                                                                                                                                        tGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaata             <  1:2695301/149‑1 (MQ=255)
                                                                                                                                        tGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaata             <  1:2185091/149‑1 (MQ=255)
                                                                                                                                         ggAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataa            >  2:4094705/1‑149 (MQ=255)
                                                                                                                                          gAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataaa           >  1:3993468/1‑149 (MQ=255)
                                                                                                                                           aaaaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat          <  2:3656183/149‑1 (MQ=255)
                                                                                                                                           aaaaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat          <  2:1058477/149‑1 (MQ=255)
                                                                                                                                           aaaaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat          >  2:5016479/1‑149 (MQ=255)
                                                                                                                                            aaaaGGGTATCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAaataaata         <  2:1446523/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTA  >  CP009273/3574872‑3575167

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: