Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,575,032 C→A 41.8% 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,0320CA41.8% 313.2 / inf 306intergenic (+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 (87/91);  new base A (58/70);  total (145/161)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 309 total aligned reads displayed.

TCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATA  >  CP009273/3574886‑3575179
                                                                                                                                                  |                                                                                                                                                   
tCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCg                                                                                                                                                    <  2:2985836/148‑1 (MQ=255)
 cTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGa                                                                                                                                                   <  1:105877/148‑1 (MQ=255)
  tataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGAc                                                                                                                                                  >  2:2327675/1‑148 (MQ=255)
  tataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGAc                                                                                                                                                  <  2:3703541/148‑1 (MQ=255)
   ataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACa                                                                                                                                                 >  1:2649733/1‑148 (MQ=255)
   ataAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACa                                                                                                                                                 <  1:2841584/148‑1 (MQ=255)
    taAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACAt                                                                                                                                                <  2:2867551/148‑1 (MQ=255)
     aaaaaaTACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATa                                                                                                                                               >  1:2765604/1‑148 (MQ=255)
           tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcc                                                                                                                                         <  1:1124217/148‑1 (MQ=255)
           tACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATAccttcc                                                                                                                                         <  1:1926139/148‑1 (MQ=255)
                cTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTaa                                                                                                                                    >  1:3759247/1‑148 (MQ=255)
                 tttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAg                                                                                                                                   <  1:1373738/148‑1 (MQ=255)
                 tttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAg                                                                                                                                   >  1:3293539/1‑148 (MQ=255)
                 tttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAg                                                                                                                                   >  1:3293092/1‑148 (MQ=255)
                  ttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAgg                                                                                                                                  >  1:947605/1‑148 (MQ=255)
                  ttGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAgg                                                                                                                                  >  1:3653878/1‑148 (MQ=255)
                     aaaTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGtat                                                                                                                               <  2:1208531/148‑1 (MQ=255)
                       aTTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtattt                                                                                                                             >  1:3397029/1‑148 (MQ=255)
                         tGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGtatttat                                                                                                                           <  2:1819999/148‑1 (MQ=255)
                                aCTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                                    <  1:774457/148‑1 (MQ=255)
                                aCTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCa                                                                                                                    >  1:1553355/1‑148 (MQ=255)
                                  tGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAAt                                                                                                                  <  2:768444/148‑1 (MQ=255)
                                  tGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAAt                                                                                                                  <  2:2685416/148‑1 (MQ=255)
                                   ggTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATa                                                                                                                 >  1:596525/1‑148 (MQ=255)
                                      aagaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACat                                                                                                              >  1:486717/1‑148 (MQ=255)
                                      aagaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACat                                                                                                              <  1:1600537/148‑1 (MQ=255)
                                        gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatat                                                                                                            >  1:2523795/1‑148 (MQ=255)
                                        gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatat                                                                                                            <  2:583010/148‑1 (MQ=255)
                                        gaagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatat                                                                                                            <  1:399364/148‑1 (MQ=255)
                                         aagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                                           <  1:1407254/148‑1 (MQ=255)
                                         aagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                                           >  2:3616793/1‑148 (MQ=255)
                                         aagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                                           >  2:3232893/1‑148 (MQ=255)
                                         aagTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatata                                                                                                           <  1:118316/148‑1 (MQ=255)
                                          agTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACatatat                                                                                                          <  1:2501373/148‑1 (MQ=255)
                                           gTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATAtt                                                                                                         >  2:2145425/1‑148 (MQ=255)
                                           gTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATAtt                                                                                                         >  2:1003449/1‑148 (MQ=255)
                                              cGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTgag                                                                                                      >  1:1788731/1‑148 (MQ=255)
                                               gCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTgaga                                                                                                     <  1:236336/148‑1 (MQ=255)
                                                 aTTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGAtt                                                                                                   <  1:3516045/148‑1 (MQ=255)
                                                 aTTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGAtt                                                                                                   <  2:3521446/148‑1 (MQ=255)
                                                   tGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACATTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTa                                                                                                 >  2:2064711/1‑148 (MQ=255)
                                                        aaTGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCa                                                                                            >  2:305294/1‑148 (MQ=255)
                                                          tGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATa                                                                                          >  2:2037057/1‑148 (MQ=255)
                                                          tGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATa                                                                                          <  1:18941/148‑1 (MQ=255)
                                                           ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                                         <  1:730180/148‑1 (MQ=255)
                                                           ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                                         <  1:866942/148‑1 (MQ=255)
                                                           ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                                         >  2:680481/1‑148 (MQ=255)
                                                           ggTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATaa                                                                                         <  1:2032166/148‑1 (MQ=255)
                                                              gAATATGTCGATTCTTATTATTTGGGGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGAAATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTa                                                                                      <  1:2262885/148‑1 (MQ=255)
                                                                       gATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                             <  1:1342514/148‑1 (MQ=255)
                                                                        aTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAg                                                                            >  2:2475479/1‑148 (MQ=255)
                                                                          gCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                          <  2:709129/148‑1 (MQ=255)
                                                                             tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                                       <  2:1479817/148‑1 (MQ=255)
                                                                             tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                                       <  2:1936418/148‑1 (MQ=255)
                                                                              attatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                      >  2:25308/1‑148 (MQ=255)
                                                                                tataTGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAatat                                                                    <  2:169686/148‑1 (MQ=255)
                                                                                tataTGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAatat                                                                    <  1:3205432/148‑1 (MQ=255)
                                                                                      gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttatt                                                              <  2:2796767/148‑1 (MQ=255)
                                                                                      gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttatt                                                              <  2:1632860/148‑1 (MQ=255)
                                                                                      gcgcGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttatt                                                              <  2:3578540/148‑1 (MQ=255)
                                                                                       cgcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                             <  2:369545/148‑1 (MQ=255)
                                                                                       cgcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattttattt                                                             <  2:276763/148‑1 (MQ=255)
                                                                                        gcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                            <  1:3689358/148‑1 (MQ=255)
                                                                                        gcgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATattttatttt                                                            <  1:1010090/148‑1 (MQ=255)
                                                                                            gTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTaa                                                        >  1:2883698/1‑148 (MQ=255)
                                                                                                aGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatata                                                    >  2:2565/1‑148 (MQ=255)
                                                                                                     gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGc                                               >  2:2945200/1‑148 (MQ=255)
                                                                                                     gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGc                                               >  2:966589/1‑148 (MQ=255)
                                                                                                      attgTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTACTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTGATTATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGcc                                              <  1:1165507/145‑1 (MQ=255)
                                                                                                      aTAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGcc                                              >  2:2635837/1‑148 (MQ=255)
                                                                                                        aGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTa                                            <  2:594546/148‑1 (MQ=255)
                                                                                                           gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACaa                                         <  1:1375044/148‑1 (MQ=255)
                                                                                                           gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACaa                                         >  1:1906227/1‑148 (MQ=255)
                                                                                                           gCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACaa                                         >  2:42036/1‑148 (MQ=255)
                                                                                                              cTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAAttt                                      >  1:3238856/1‑148 (MQ=255)
                                                                                                               tttttCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTc                                     >  1:701770/1‑148 (MQ=255)
                                                                                                                   tCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTgg                                 <  1:2978485/148‑1 (MQ=255)
                                                                                                                   tCTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTgg                                 <  2:316862/148‑1 (MQ=255)
                                                                                                                    cTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                                <  1:2194989/148‑1 (MQ=255)
                                                                                                                    cTGAGATGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                                <  2:3041667/148‑1 (MQ=255)
                                                                                                                      gagaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGt                              <  1:2037057/148‑1 (MQ=255)
                                                                                                                       agaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGta                             >  2:1089324/1‑148 (MQ=255)
                                                                                                                        gaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGTCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaa                            <  2:92324/148‑1 (MQ=255)
                                                                                                                        gaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaa                            <  2:2234228/148‑1 (MQ=255)
                                                                                                                        gaTGGAAAAAGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaa                            >  2:2669943/1‑148 (MQ=255)
                                                                                                                            gAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataa                        >  1:570211/1‑148 (MQ=255)
                                                                                                                             aaaaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataaa                       >  2:3611514/1‑148 (MQ=255)
                                                                                                                              aaaaGGGTGTCATTCAAAATAGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTaaataaat                      <  2:8015/148‑1 (MQ=255)
                                                                                                                               aaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAaataaata                     <  2:334426/148‑1 (MQ=255)
                                                                                                                               aaaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAaataaata                     >  2:312987/1‑148 (MQ=255)
                                                                                                                                aaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGACTACAATTTCCTGGAGTAAataaataa                    >  2:3460259/1‑148 (MQ=255)
                                                                                                                                aaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATAAGCCTACAATTTCCTGGAGTAAataaataa                    >  2:3545337/1‑148 (MQ=255)
                                                                                                                                  gggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                  >  1:3559746/1‑148 (MQ=255)
                                                                                                                                  gggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACa                  >  2:2123771/1‑148 (MQ=255)
                                                                                                                                   ggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                 >  2:2324116/1‑148 (MQ=255)
                                                                                                                                      gtCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTa              <  2:396219/148‑1 (MQ=255)
                                                                                                                                                 tCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataat   <  1:3062188/148‑1 (MQ=255)
                                                                                                                                                 tCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataat   <  2:1553355/148‑1 (MQ=255)
                                                                                                                                                  cGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataata  <  1:165808/148‑1 (MQ=255)
                                                                                                                                                  cGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataata  >  2:725115/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCTATAAAAAATACGGCTTTGAAATTGAAGGGACTGGTAAGAAGTACGCATTGCGTAATGGTGAATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATA  >  CP009273/3574886‑3575179

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: