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,0320CA99.5% 630.2 / ‑4.6 189intergenic (+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 (0/1);  new base A (85/103);  total (85/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 189 total aligned reads displayed.

CTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAAT  >  CP009273/3574887‑3575178
                                                                                                                                                 |                                                                                                                                                  
  ataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATa                                                                                                                                              <  1:5097119/150‑1 (MQ=255)
   taAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAc                                                                                                                                             >  2:3206772/1‑150 (MQ=255)
      aaaaTACGTCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcctt                                                                                                                                          <  1:4793409/150‑1 (MQ=255)
      aaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCTCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcctt                                                                                                                                          >  1:2594841/1‑150 (MQ=255)
      aaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAcctt                                                                                                                                          <  1:1753211/150‑1 (MQ=255)
       aaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttc                                                                                                                                         >  1:15730/1‑150 (MQ=255)
       aaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttc                                                                                                                                         >  1:690608/1‑150 (MQ=255)
         aTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcct                                                                                                                                       >  2:4389774/1‑150 (MQ=255)
         aTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcct                                                                                                                                       <  1:4117920/150‑1 (MQ=255)
         aTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATAccttcct                                                                                                                                       >  1:2645028/1‑150 (MQ=255)
                tttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGt                                                                                                                                >  2:4873151/1‑150 (MQ=255)
                 ttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGta                                                                                                                               <  2:3181536/150‑1 (MQ=255)
                   gAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatt                                                                                                                             <  1:334234/150‑1 (MQ=255)
                   gAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatt                                                                                                                             <  2:2557124/150‑1 (MQ=255)
                   gAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtata                                                                                                                             <  2:1725101/150‑2 (MQ=255)
                    aaaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtattt                                                                                                                            >  2:2799116/1‑150 (MQ=255)
                    aaaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtattt                                                                                                                            >  2:3280816/1‑150 (MQ=255)
                    aaaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtattt                                                                                                                            >  1:1231668/1‑150 (MQ=255)
                     aaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtattta                                                                                                                           <  1:220160/150‑1 (MQ=255)
                       ttGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatttatt                                                                                                                         <  1:2473330/150‑1 (MQ=255)
                        tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        <  1:971132/150‑1 (MQ=255)
                        tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        >  1:1974487/1‑150 (MQ=255)
                        tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        >  1:4096731/1‑150 (MQ=255)
                        tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTatttattt                                                                                                                        <  1:119261/150‑1 (MQ=255)
                            gggACTGTTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGccc                                                                                                                    <  2:1213636/150‑1 (MQ=255)
                            gggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGccc                                                                                                                    <  2:4532917/150‑1 (MQ=255)
                             ggACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                                   >  2:2351970/1‑150 (MQ=255)
                                cTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATa                                                                                                                <  2:3787125/150‑1 (MQ=255)
                                 tGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATAc                                                                                                               <  2:1811788/150‑1 (MQ=255)
                                  ggTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACa                                                                                                              <  2:4104885/150‑1 (MQ=255)
                                    tAAGAAGTACTCATTGATTAATTTTTAATATTTCGATGAATATTATATGGCGCGGGTGAATTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGTGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACata                                                                                                            <  2:2624427/150‑1 (MQ=255)
                                    tAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATTGCGCGGGTGAAGTCAGATAGTGCCCTTTTTCTGAGCTGGAAAAAGGGTGTCATTAAAAATAGACATACCTTCCTATAAGGTAGTTATTTGCCCAATACata                                                                                                            >  2:3141120/1‑150 (MQ=255)
                                    tAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACata                                                                                                            >  2:2007500/1‑150 (MQ=255)
                                      agaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                                          <  1:395556/150‑1 (MQ=255)
                                        tagTACGCATTGCGTATTTGTGATTATGTCGCTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCACTTTTTCTGAGATGGAAATAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATAtt                                                                                                        <  2:5094551/149‑1 (MQ=255)
                                              gCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGAtt                                                                                                  >  1:4890159/1‑150 (MQ=255)
                                                aTTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTa                                                                                                >  2:4500460/1‑150 (MQ=255)
                                                aTTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTa                                                                                                <  1:4781283/150‑1 (MQ=255)
                                                 ttGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTaa                                                                                               >  1:2978444/1‑150 (MQ=255)
                                                   gCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACt                                                                                             <  2:3348942/150‑1 (MQ=255)
                                                      tAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAt                                                                                          >  2:4679769/1‑150 (MQ=255)
                                                       aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATa                                                                                         <  1:3373438/150‑1 (MQ=255)
                                                       aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATa                                                                                         <  1:3213594/150‑1 (MQ=255)
                                                        aTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                                        >  1:3686056/1‑150 (MQ=255)
                                                          ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                                      <  1:1093461/150‑1 (MQ=255)
                                                          ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAAtt                                                                                      <  1:1093444/150‑1 (MQ=255)
                                                            tGAATATGTCGAGGCATATTATAGGGCGCGGGTGAAGTAAGTTAGTGCCCTTTTTCTGAGATGGAATAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACGCATAATTAc                                                                                    <  2:3369779/150‑1 (MQ=255)
                                                             gAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACt                                                                                   <  2:1860470/150‑1 (MQ=255)
                                                                 atGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCt                                                                               <  2:3086098/150‑1 (MQ=255)
                                                                   gTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaa                                                                             <  2:1850744/150‑1 (MQ=255)
                                                                    tCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                            <  1:3415666/150‑1 (MQ=255)
                                                                      gATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                          >  1:3289657/1‑150 (MQ=255)
                                                                       aTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                         <  2:3961167/150‑1 (MQ=255)
                                                                        tGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtgt                                                                        <  1:2484734/150‑1 (MQ=255)
                                                                           atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                     <  2:3869275/150‑1 (MQ=255)
                                                                           atatTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                     >  1:2833825/1‑150 (MQ=255)
                                                                             attatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAatat                                                                   >  1:2409837/1‑150 (MQ=255)
                                                                                atatGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtatttt                                                                >  1:1388263/1‑150 (MQ=255)
                                                                                 tatGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtatttta                                                               >  2:2844163/1‑150 (MQ=255)
                                                                                   tGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttatt                                                             >  1:2133204/1‑150 (MQ=255)
                                                                                     gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                           <  2:2897200/150‑1 (MQ=255)
                                                                                        cgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTa                                                        >  2:1539844/1‑150 (MQ=255)
                                                                                             gAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatata                                                   <  1:3206694/150‑1 (MQ=255)
                                                                                               aGTAAGATAGTGCCCCTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAAtatata                                                 >  1:4904098/1‑150 (MQ=255)
                                                                                                gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATAc                                                >  2:762049/1‑150 (MQ=255)
                                                                                                 tAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACg                                               >  2:1391829/1‑150 (MQ=255)
                                                                                                 tAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACg                                               <  2:4489898/150‑1 (MQ=255)
                                                                                                 tAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACg                                               >  1:4658305/1‑150 (MQ=255)
                                                                                                    gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCt                                            >  1:3563151/1‑150 (MQ=255)
                                                                                                     aTAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGAATTTATATGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATAAGCCTa                                           >  1:2237241/1‑150 (MQ=255)
                                                                                                       aGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACa                                         <  1:4389774/150‑1 (MQ=255)
                                                                                                           cccTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAttt                                     <  2:4284985/150‑1 (MQ=255)
                                                                                                           cccTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAttt                                     <  2:798262/150‑1 (MQ=255)
                                                                                                           cccTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAttt                                     <  1:2560466/150‑1 (MQ=255)
                                                                                                             cTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTcc                                   >  2:1674554/1‑150 (MQ=255)
                                                                                                             cTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTcc                                   <  2:4904157/150‑1 (MQ=255)
                                                                                                             cTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTcc                                   <  1:3017017/150‑1 (MQ=255)
                                                                                                             cTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTcc                                   <  1:4873151/150‑1 (MQ=255)
                                                                                                                 ttCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                               <  1:1326715/150‑1 (MQ=255)
                                                                                                                 ttCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                               <  1:4619505/150‑1 (MQ=255)
                                                                                                                  tCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAg                              >  1:4937147/1‑150 (MQ=255)
                                                                                                                  tCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAg                              <  1:2799116/150‑1 (MQ=255)
                                                                                                                     gagaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaa                           >  1:2112709/1‑150 (MQ=255)
                                                                                                                     gagaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaa                           >  1:392433/1‑150 (MQ=255)
                                                                                                                       taTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaat                         <  1:2597390/149‑1 (MQ=255)
                                                                                                                           gAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat                     >  1:2330237/1‑150 (MQ=255)
                                                                                                                             aaaaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAataaataa                   >  1:3253359/1‑150 (MQ=255)
                                                                                                                              aaaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAAc                  >  1:3392187/1‑150 (MQ=255)
                                                                                                                               aaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                 >  1:539060/1‑150 (MQ=255)
                                                                                                                                   gtgtCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTa             >  1:5133605/1‑150 (MQ=255)
                                                                                                                                   gtgtCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTa             >  2:1345164/1‑150 (MQ=255)
                                                                                                                                   gtgtCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTa             >  2:458418/1‑150 (MQ=255)
                                                                                                                                       cATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACaa         >  1:983308/1‑150 (MQ=255)
                                                                                                                                         ttCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGc       <  1:3095597/150‑1 (MQ=255)
                                                                                                                                           cAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat     <  2:1022533/150‑1 (MQ=255)
                                                                                                                                           cAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat     <  2:3605898/150‑1 (MQ=255)
                                                                                                                                           cAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCat     <  2:3480680/150‑1 (MQ=255)
                                                                                                                                            aaaaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCata    <  2:1348205/150‑1 (MQ=255)
                                                                                                                                             aaaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataa   <  2:144505/150‑1 (MQ=255)
                                                                                                                                              aaTAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataat  >  2:2576649/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAAT  >  CP009273/3574887‑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: