Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,575,032 C→A 100% 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,0320CA98.1% 1250.8 / 9.9 389intergenic (+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 (5/3);  new base A (169/212);  total (174/215)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 389 total aligned reads displayed.

GTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAAT  >  CP009273/3574885‑3575178
                                                                                                                                                   |                                                                                                                                                  
 tCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGAc                                                                                                                                                 <  2:5276323/150‑1 (MQ=255)
  cTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACa                                                                                                                                                <  1:3893532/150‑1 (MQ=255)
   tataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACAt                                                                                                                                               <  1:3922061/150‑1 (MQ=255)
     taAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAc                                                                                                                                             >  2:1786861/1‑150 (MQ=255)
     taAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAc                                                                                                                                             >  1:4319805/1‑150 (MQ=255)
      aaaaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcc                                                                                                                                            <  2:1920696/150‑1 (MQ=255)
       aaaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcct                                                                                                                                           <  2:4165768/150‑1 (MQ=255)
       aaaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcct                                                                                                                                           <  2:610439/150‑1 (MQ=255)
        aaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcctt                                                                                                                                          <  1:333313/150‑1 (MQ=255)
           aTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcct                                                                                                                                       <  2:5299743/150‑1 (MQ=255)
            tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcctt                                                                                                                                      >  1:2270278/1‑150 (MQ=255)
             aCGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCttt                                                                                                                                     <  1:4581949/150‑1 (MQ=255)
               ggCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTaa                                                                                                                                   >  2:3359283/1‑150 (MQ=255)
                    tGAAATTGAAGGGACTGGTAATAAGTACGCATTGAGTAATGGTGAATATGTCTATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtat                                                                                                                              <  1:5276410/150‑1 (MQ=255)
                      aaaTTGAAGTGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAACTTCCTTTAAGGtattt                                                                                                                            >  2:2989264/1‑150 (MQ=255)
                       aaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtattta                                                                                                                           >  2:2442349/1‑150 (MQ=255)
                        aTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatttat                                                                                                                          >  2:2920411/1‑150 (MQ=255)
                        aTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatttat                                                                                                                          >  2:3826237/1‑150 (MQ=255)
                         ttGTAGGGACTGGTTAGAAGTACGGAGTGCGTAATGGTGAAGCTGTCGATGCATATTATATGGAGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGTGCTGGCACAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatttatt                                                                                                                         <  2:2668623/150‑1 (MQ=255)
                          tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        <  2:131026/150‑1 (MQ=255)
                          tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        <  2:2243301/150‑1 (MQ=255)
                          tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        <  2:4070134/150‑1 (MQ=255)
                            aaGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGc                                                                                                                      >  2:4629189/1‑150 (MQ=255)
                               ggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                                   <  1:1092676/150‑1 (MQ=255)
                                 aCTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAAt                                                                                                                 <  2:4793585/150‑1 (MQ=255)
                                   tGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATAc                                                                                                               <  2:5023934/150‑1 (MQ=255)
                                      tAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACata                                                                                                            <  2:1706317/150‑1 (MQ=255)
                                         gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACAGACCTTCCTTTAAGGTATTTATTTGCCCAATACatatat                                                                                                         <  2:4872487/150‑1 (MQ=255)
                                           agTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTg                                                                                                       <  1:2685573/150‑1 (MQ=255)
                                           agTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACAGACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTg                                                                                                       <  1:880974/150‑1 (MQ=255)
                                             tACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTgag                                                                                                     <  1:737973/150‑1 (MQ=255)
                                                    tGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATANTGAGATTTAAc                                                                                              <  1:2766544/150‑1 (MQ=255)
                                                     gCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACt                                                                                             <  2:4976731/150‑1 (MQ=255)
                                                       gTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCa                                                                                           <  1:4280118/150‑1 (MQ=255)
                                                        tAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAt                                                                                          <  2:365071/150‑1 (MQ=255)
                                                           tGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAt                                                                                       >  2:3760621/1‑150 (MQ=255)
                                                            ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                                      <  1:1034864/150‑1 (MQ=255)
                                                            ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                                      <  1:2536343/150‑1 (MQ=255)
                                                            ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                                      >  2:4551729/1‑150 (MQ=255)
                                                             gTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTa                                                                                     >  1:1349916/1‑150 (MQ=255)
                                                               gAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACt                                                                                   >  2:2296254/1‑150 (MQ=255)
                                                                   atGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCt                                                                               <  2:2548705/150‑1 (MQ=255)
                                                                      tCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                            <  1:25643/150‑1 (MQ=255)
                                                                       cGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAg                                                                           <  2:4330484/150‑1 (MQ=255)
                                                                       cGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAg                                                                           >  2:4505013/1‑150 (MQ=255)
                                                                        gATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                          >  2:3098573/1‑150 (MQ=255)
                                                                         aTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                         <  2:3440415/150‑1 (MQ=255)
                                                                         aTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                         >  2:1031495/1‑150 (MQ=255)
                                                                            cATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                                      >  2:2155589/1‑150 (MQ=255)
                                                                             atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                     >  1:3525974/1‑150 (MQ=255)
                                                                             atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                     <  2:5121370/150‑1 (MQ=255)
                                                                                ttatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtatt                                                                  <  1:3196134/150‑1 (MQ=255)
                                                                                 tataTGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattt                                                                 >  1:1061868/1‑150 (MQ=255)
                                                                                    atGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttat                                                              >  2:2464313/1‑150 (MQ=255)
                                                                                    atGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttat                                                              >  2:1404866/1‑150 (MQ=255)
                                                                                      ggCGCGGGTGAAGTAAGATAGTGCACTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                            <  1:4158802/150‑1 (MQ=255)
                                                                                       gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                           <  2:820007/150‑1 (MQ=255)
                                                                                       gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                           <  1:5321054/150‑1 (MQ=255)
                                                                                            ggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAat                                                      <  2:4851548/150‑1 (MQ=255)
                                                                                              tGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatat                                                    <  1:3752797/150‑1 (MQ=255)
                                                                                               gAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatata                                                   >  1:2147599/1‑150 (MQ=255)
                                                                                                aaGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatatat                                                  >  1:1753581/1‑150 (MQ=255)
                                                                                                aaGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatatat                                                  <  1:4388046/150‑1 (MQ=255)
                                                                                                  gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATAc                                                >  1:1655268/1‑150 (MQ=255)
                                                                                                  gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATAc                                                <  1:4904699/150‑1 (MQ=255)
                                                                                                   tAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACg                                               <  2:4382764/150‑1 (MQ=255)
                                                                                                      gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATTTTTTATTTTTTAATATATACGCCt                                            <  1:1506280/150‑1 (MQ=255)
                                                                                                      gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCt                                            <  1:2076894/150‑1 (MQ=255)
                                                                                                        tAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTAc                                          >  1:972392/1‑150 (MQ=255)
                                                                                                         aGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTATGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACa                                         >  2:1568301/1‑150 (MQ=255)
                                                                                                         aGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACa                                         <  1:1837043/150‑1 (MQ=255)
                                                                                                            gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAtt                                      >  2:4190230/1‑150 (MQ=255)
                                                                                                             cccTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAttt                                     >  2:1999191/1‑150 (MQ=255)
                                                                                                                 ttttCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTg                                 >  1:3817726/1‑150 (MQ=255)
                                                                                                                   ttCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                               >  2:3150449/1‑150 (MQ=255)
                                                                                                                     cTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGt                             <  1:4468569/150‑1 (MQ=255)
                                                                                                                      tGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGta                            >  2:283421/1‑150 (MQ=255)
                                                                                                                        agaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaa                          <  1:40758/150‑1 (MQ=255)
                                                                                                                        agaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaa                          >  2:3905701/1‑150 (MQ=255)
                                                                                                                         gaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaat                         >  1:160306/1‑150 (MQ=255)
                                                                                                                          aTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaata                        >  1:3287849/1‑150 (MQ=255)
                                                                                                                           tGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataa                       >  2:2095161/1‑150 (MQ=255)
                                                                                                                             gAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat                     >  1:386216/1‑150 (MQ=255)
                                                                                                                               aaaaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAataaataa                   >  1:2011485/1‑150 (MQ=255)
                                                                                                                                 aaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                 <  2:1877536/150‑1 (MQ=255)
                                                                                                                                 aaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                 >  1:59095/1‑150 (MQ=255)
                                                                                                                                 aaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                 >  2:1520541/1‑150 (MQ=255)
                                                                                                                                  aGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                <  1:4766985/150‑1 (MQ=255)
                                                                                                                                  aGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                <  1:5143564/150‑1 (MQ=255)
                                                                                                                                  aGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                >  1:559701/1‑150 (MQ=255)
                                                                                                                                   gggTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAAt               <  1:5021880/150‑1 (MQ=255)
                                                                                                                                    ggTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAAtt              <  1:2798440/150‑1 (MQ=255)
                                                                                                                                    ggTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAAtt              <  1:3681984/150‑1 (MQ=255)
                                                                                                                                          aTTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAg        >  2:936246/1‑150 (MQ=255)
                                                                                                                                            tCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCa      <  2:3524352/150‑1 (MQ=255)
                                                                                                                                             cAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat     <  1:1201533/150‑1 (MQ=255)
                                                                                                                                             cAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat     <  2:661371/150‑1 (MQ=255)
                                                                                                                                                aaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataat  <  1:2989264/150‑1 (MQ=255)
                                                                                                                                                aaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataat  <  2:5082774/150‑1 (MQ=255)
                                                                                                                                                aaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataat  >  2:1673900/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAAT  >  CP009273/3574885‑3575178

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: