Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,575,032 C→A 33.7% intergenic (+46/‑191) 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,0320CA33.7% 406.8 / 221.8 258intergenic (+46/‑191)yhhY/yhhZputative acetyltransferase/putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
Reads supporting (aligned to +/- strand):  ref base C (67/104);  new base A (50/37);  total (117/141)
Fisher's exact test for biased strand distribution p-value = 5.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

Only 100 of 262 total aligned reads displayed.

GGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATA  >  CP009273/3574884‑3575179
                                                                                                                                                    |                                                                                                                                                   
ggTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATc                                                                                                                                                     <  1:3507367/149‑1 (MQ=255)
 gTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCg                                                                                                                                                    >  1:3718222/1‑149 (MQ=255)
  tCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGa                                                                                                                                                   >  2:1287090/1‑149 (MQ=255)
    tataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACa                                                                                                                                                 >  2:786915/1‑149 (MQ=255)
      taAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATa                                                                                                                                               >  2:1088945/1‑149 (MQ=255)
       aaaaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAc                                                                                                                                              >  1:828388/1‑149 (MQ=255)
        aaaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAcc                                                                                                                                             >  1:2835808/1‑149 (MQ=255)
        aaaaaTACGGATTTTAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAcc                                                                                                                                             <  1:1774174/149‑1 (MQ=255)
         aaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAcct                                                                                                                                            <  2:3262439/149‑1 (MQ=255)
             tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcct                                                                                                                                        <  2:1899659/149‑1 (MQ=255)
             tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcct                                                                                                                                        <  2:3602192/149‑1 (MQ=255)
             tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcct                                                                                                                                        <  1:2913346/149‑1 (MQ=255)
             tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcct                                                                                                                                        <  1:281197/149‑1 (MQ=255)
              aCGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcctt                                                                                                                                       <  1:2889471/149‑1 (MQ=255)
                    ttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACAGACCTTCCTTTAAGGt                                                                                                                                 <  2:1348648/149‑1 (MQ=255)
                     tGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGta                                                                                                                                <  2:146468/149‑1 (MQ=255)
                       aaaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatt                                                                                                                              >  1:2004294/1‑149 (MQ=255)
                        aaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtattt                                                                                                                             <  2:2534948/149‑1 (MQ=255)
                        aaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtattt                                                                                                                             <  1:136349/149‑1 (MQ=255)
                         aTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtattta                                                                                                                            >  1:258416/1‑149 (MQ=255)
                            gAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTatttattt                                                                                                                         >  1:3650584/1‑149 (MQ=255)
                            gAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                         <  2:3964705/149‑1 (MQ=255)
                             aaGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTg                                                                                                                        <  2:1500288/149‑1 (MQ=255)
                             aaGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTg                                                                                                                        >  1:2293668/1‑149 (MQ=255)
                              aGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGc                                                                                                                       <  1:610565/149‑1 (MQ=255)
                               gggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGcc                                                                                                                      >  2:1534965/1‑149 (MQ=255)
                                ggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGccc                                                                                                                     <  2:3225069/149‑1 (MQ=255)
                                ggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGccc                                                                                                                     <  1:35679/149‑1 (MQ=255)
                                 gACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                                    <  1:460855/149‑1 (MQ=255)
                                 gACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                                    <  1:791015/149‑1 (MQ=255)
                                   cTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAAt                                                                                                                  <  1:2213254/149‑1 (MQ=255)
                                      gTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACa                                                                                                               >  2:3345904/1‑149 (MQ=255)
                                       tAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACat                                                                                                              <  1:319582/149‑1 (MQ=255)
                                        aagaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACata                                                                                                             <  1:358694/149‑1 (MQ=255)
                                         agaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatat                                                                                                            <  2:2659547/149‑1 (MQ=255)
                                          gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                                           <  1:3318578/149‑1 (MQ=255)
                                           aagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatatat                                                                                                          <  1:216416/149‑1 (MQ=255)
                                            agTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATAtt                                                                                                         >  2:3525795/1‑149 (MQ=255)
                                             gTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTg                                                                                                        >  1:3927658/1‑149 (MQ=255)
                                                  cATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGAtt                                                                                                   <  1:3773730/149‑1 (MQ=255)
                                                   aTTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTCTTTATTTGCCCAATACATATATTGAGAttt                                                                                                  >  1:1542673/1‑149 (MQ=255)
                                                   aTTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGAttt                                                                                                  <  2:2140800/149‑1 (MQ=255)
                                                          aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAt                                                                                           <  2:3834956/149‑1 (MQ=255)
                                                          aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAt                                                                                           <  2:1791781/149‑1 (MQ=255)
                                                            nggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                                         >  1:3429465/2‑149 (MQ=255)
                                                              gTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                                       >  2:680039/1‑149 (MQ=255)
                                                                    atGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTcc                                                                                 <  2:1357670/149‑1 (MQ=255)
                                                                     tgtataTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCt                                                                                <  2:2376225/144‑1 (MQ=255)
                                                                      gTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTa                                                                               >  2:3393709/1‑149 (MQ=255)
                                                                        cGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                             <  1:1016760/149‑1 (MQ=255)
                                                                        cGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                             <  2:3927658/149‑1 (MQ=255)
                                                                         gATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAg                                                                            >  2:4024395/1‑149 (MQ=255)
                                                                          aTGCATATTATATGGCGCGGGTGAAGTAAGTTTGTGCCCTTTTTCTGAGAGGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                           <  2:2156627/149‑1 (MQ=255)
                                                                           tGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                          >  2:1567266/1‑149 (MQ=255)
                                                                             cATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTa                                                                        <  1:3845083/149‑1 (MQ=255)
                                                                             cATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTa                                                                        <  1:3997278/149‑1 (MQ=255)
                                                                             cATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTa                                                                        >  2:2029275/1‑149 (MQ=255)
                                                                              atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCCATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                                       <  2:1972764/149‑1 (MQ=255)
                                                                              atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                                       <  1:3875335/149‑1 (MQ=255)
                                                                               tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                      >  1:2704582/1‑149 (MQ=255)
                                                                                     atGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtagtttg                                                                >  2:1879202/1‑148 (MQ=255)
                                                                                      tGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttat                                                               <  2:2740015/149‑1 (MQ=255)
                                                                                      tGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttat                                                               <  2:2647973/149‑1 (MQ=255)
                                                                                        gcgcGGGTGAAGTGAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                             >  2:529981/1‑149 (MQ=255)
                                                                                        gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                             >  2:1474827/1‑149 (MQ=255)
                                                                                        gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                             <  1:1534965/149‑1 (MQ=255)
                                                                                        gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                             <  1:4046084/149‑1 (MQ=255)
                                                                                        gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATAAATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtagtttattt                                                             >  2:3457756/1‑149 (MQ=255)
                                                                                         cgcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                            <  2:2763273/149‑1 (MQ=255)
                                                                                          gcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTAttttt                                                           <  1:305827/149‑1 (MQ=255)
                                                                                          gcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTAttttt                                                           <  1:2377262/149‑1 (MQ=255)
                                                                                            gggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTa                                                         >  2:3573455/1‑149 (MQ=255)
                                                                                             ggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTaa                                                        >  1:3891872/1‑149 (MQ=255)
                                                                                             ggTGAACTAAGATACGGCCCTTTTTCTGATCTTGAAAAAGGGTGTCAATCAAAATCCTCATTCCTTACTTTAATGTATTTATTTGCCCAATACAAATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTaa                                                        <  1:1693203/149‑1 (MQ=255)
                                                                                              gTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCCATACATATATTGAGATTTAACTCATAATTACTTCATAAAGTGTAATATTTTATTTTTTAat                                                       >  2:1532920/1‑149 (MQ=255)
                                                                                                   gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAAtatata                                                  <  1:1088945/149‑1 (MQ=255)
                                                                                                     aaGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACg                                                >  1:2218126/1‑149 (MQ=255)
                                                                                                        aTAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCt                                             >  1:3138346/1‑149 (MQ=255)
                                                                                                        aTAGTGCCCTTTTTCGGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCt                                             <  2:2215517/149‑1 (MQ=255)
                                                                                                           gTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACa                                          >  1:1803447/1‑149 (MQ=255)
                                                                                                             gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAt                                        >  1:3001329/1‑149 (MQ=255)
                                                                                                              cccTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAtt                                       <  2:233650/149‑1 (MQ=255)
                                                                                                                    ttCTGAGATGGAAAAAGGGTGACATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTgg                                 <  1:2835052/149‑1 (MQ=255)
                                                                                                                     tCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                                <  1:3573455/149‑1 (MQ=255)
                                                                                                                       tGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGt                              <  1:3369466/149‑1 (MQ=255)
                                                                                                                          gaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaa                           <  2:3768763/149‑1 (MQ=255)
                                                                                                                           aTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaat                          <  2:1020629/149‑1 (MQ=255)
                                                                                                                            tGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaata                         >  1:1182190/1‑149 (MQ=255)
                                                                                                                            tGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaata                         >  1:649750/1‑149 (MQ=255)
                                                                                                                               aaaaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat                      <  1:1195504/149‑1 (MQ=255)
                                                                                                                                 aaaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACNCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAataaataa                    >  1:2084731/1‑149 (MQ=255)
                                                                                                                                   aGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                  <  2:1704688/149‑1 (MQ=255)
                                                                                                                                      gtgtCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAAtt               >  1:977596/1‑149 (MQ=255)
                                                                                                                                         tCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAAc            >  2:83374/1‑149 (MQ=255)
                                                                                                                                             tCAAAATCGACATACCTTCATTTAAGGTATTTATTTGCCAAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGc        <  1:3251564/149‑1 (MQ=255)
                                                                                                                                              cAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCa       <  2:2004294/149‑1 (MQ=255)
                                                                                                                                               aaaaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat      >  1:2609903/1‑149 (MQ=255)
                                                                                                                                               aaaaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat      <  2:286020/149‑1 (MQ=255)
                                                                                                                                                aaaTCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCata     >  2:2871250/1‑149 (MQ=255)
                                                                                                                                                   tCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataata  <  2:2834236/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATA  >  CP009273/3574884‑3575179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: