Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,474,418 G→T 89.9% S413* (TCA→TAA)  emrY ← putative multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,474,4180GT89.9% 672.9 / 68.4 237S413* (TCA→TAA) emrYputative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (16/8);  new base T (95/118);  total (111/126)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 239 total aligned reads displayed.

CTGAAAGCGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAATGATGCCTGTAAAATCA  >  CP009273/2474269‑2474567
                                                                                                                                                     |                                                                                                                                                     
cTGAAAGCGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTg                                                                                                                                                       <  1:2195395/150‑1 (MQ=255)
    aaGCGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACaa                                                                                                                                                   >  1:1264865/1‑150 (MQ=255)
     aGCGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAg                                                                                                                                                  >  2:611238/1‑150 (MQ=255)
      gCGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAgg                                                                                                                                                 >  2:77551/1‑150 (MQ=255)
       cGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGt                                                                                                                                                >  1:4061336/1‑150 (MQ=255)
        gAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTa                                                                                                                                               >  1:1144113/1‑150 (MQ=255)
        gAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTa                                                                                                                                               <  2:2142999/150‑1 (MQ=255)
        gAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTa                                                                                                                                               >  2:2154808/1‑150 (MQ=255)
          aCCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACg                                                                                                                                             <  1:3553498/150‑1 (MQ=255)
           ccATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGa                                                                                                                                            >  1:1005665/1‑150 (MQ=255)
               aataTTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGaaaaa                                                                                                                                        <  1:3816384/150‑1 (MQ=255)
                atatTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGaaaaaa                                                                                                                                       >  2:4115030/1‑150 (MQ=255)
                 tatTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGaaaaaaa                                                                                                                                      >  1:3529328/1‑150 (MQ=255)
                  atTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAAt                                                                                                                                     <  1:2067623/150‑1 (MQ=255)
                    ttATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTAc                                                                                                                                              <  2:2342916/139‑1 (MQ=255)
                    ttATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTa                                                                                                                                   >  2:105496/1‑150 (MQ=255)
                    ttATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTa                                                                                                                                   >  2:4182784/1‑150 (MQ=255)
                    ttATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTa                                                                                                                                   >  1:2284874/1‑150 (MQ=255)
                     tATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTAc                                                                                                                                  <  2:2297977/150‑1 (MQ=255)
                     tATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTAc                                                                                                                                  >  1:173011/1‑150 (MQ=255)
                       tCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTc                                                                                                                                <  1:538628/150‑1 (MQ=255)
                       tCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTc                                                                                                                                >  1:916457/1‑150 (MQ=255)
                          aTAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATa                                                                                                                             >  2:4041248/1‑150 (MQ=255)
                           tAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTTACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATAc                                                                                                                            >  1:3783589/1‑150 (MQ=255)
                             aTTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTc                                                                                                                          >  2:1493096/1‑150 (MQ=255)
                             aTTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTc                                                                                                                          <  2:2203373/150‑1 (MQ=255)
                             aTTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTc                                                                                                                          <  2:2203347/150‑1 (MQ=255)
                              tttGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCg                                                                                                                         >  1:3834950/1‑150 (MQ=255)
                              tttGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCg                                                                                                                         <  1:2431158/150‑1 (MQ=255)
                                  tgACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGc                                                                                                                     >  1:997310/1‑150 (MQ=255)
                                   gACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCa                                                                                                                    <  1:695424/150‑1 (MQ=255)
                                    aCGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCAt                                                                                                                   >  2:2488258/1‑150 (MQ=255)
                                     cGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCAtt                                                                                                                  >  1:2480472/1‑150 (MQ=255)
                                       aTGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTgg                                                                                                                <  2:94158/150‑1 (MQ=255)
                                         gAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCa                                                                                                              >  1:325043/1‑150 (MQ=255)
                                          aaTTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCaa                                                                                                             >  2:3074333/1‑150 (MQ=255)
                                            ttAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAAt                                                                                                           >  2:1178487/1‑150 (MQ=255)
                                            ttAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAAt                                                                                                           >  2:3376274/1‑150 (MQ=255)
                                            ttAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAATTACTCATACTCGAGGCATTGGCAAAtt                                                                                                          >  1:1357959/1‑150 (MQ=255)
                                               aaCACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATtta                                                                                                        >  2:3363602/1‑150 (MQ=255)
                                                 cacGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATttatt                                                                                                      >  2:2904004/1‑150 (MQ=255)
                                                       gTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTgg                                                                                                >  1:885417/1‑150 (MQ=255)
                                                          aaaTTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCaa                                                                                             >  1:3927738/1‑150 (MQ=255)
                                                           aaTTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAg                                                                                            >  2:3949352/1‑150 (MQ=255)
                                                             ttGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGcc                                                                                          >  2:472939/1‑150 (MQ=255)
                                                             ttGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGACATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGcc                                                                                          <  2:2590552/150‑1 (MQ=255)
                                                              tGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCt                                                                                         <  2:188848/150‑1 (MQ=255)
                                                              tGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCt                                                                                         <  1:3365108/150‑1 (MQ=255)
                                                              tGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCt                                                                                         >  1:3199075/1‑150 (MQ=255)
                                                               gATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTg                                                                                        >  1:769506/1‑150 (MQ=255)
                                                                  cGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGaaa                                                                                     >  1:4081318/1‑150 (MQ=255)
                                                                  cGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGaaa                                                                                     <  2:4304974/150‑1 (MQ=255)
                                                                   gATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGaaaa                                                                                    <  2:1852100/150‑1 (MQ=255)
                                                                   gATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGaaaa                                                                                    >  2:825330/1‑150 (MQ=255)
                                                                    aTGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAAc                                                                                   >  2:551771/1‑150 (MQ=255)
                                                                      ggTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGa                                                                                 <  1:3508114/150‑1 (MQ=255)
                                                                       gTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGaa                                                                                >  1:668259/1‑150 (MQ=255)
                                                                       gTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGaa                                                                                >  1:3069628/1‑150 (MQ=255)
                                                                         tGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAAt                                                                              <  1:2617105/150‑1 (MQ=255)
                                                                          gCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATc                                                                             <  1:131004/150‑1 (MQ=255)
                                                                           cTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCg                                                                            <  1:3021498/150‑1 (MQ=255)
                                                                           cTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCg                                                                            <  1:1255764/150‑1 (MQ=255)
                                                                            tGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCgt                                                                           <  1:291243/150‑1 (MQ=255)
                                                                                aaCTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCgttgtt                                                                       >  2:3764571/1‑150 (MQ=255)
                                                                                aaCTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCgttgtt                                                                       <  1:457038/150‑1 (MQ=255)
                                                                                      cTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGt                                                                 >  1:1714823/1‑150 (MQ=255)
                                                                                        agggTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTaa                                                               <  2:3992518/148‑1 (MQ=255)
                                                                                        aTGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTaa                                                               <  1:419355/150‑1 (MQ=255)
                                                                                           gtgtAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTaaaaa                                                            <  2:1264865/150‑1 (MQ=255)
                                                                                            tgtAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAg                                                           <  2:3156744/150‑1 (MQ=255)
                                                                                             gtAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTCAAAAGa                                                          >  1:4051359/1‑150 (MQ=255)
                                                                                             gtAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGa                                                          <  1:1216731/150‑1 (MQ=255)
                                                                                                   aTTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAgg                                                    >  2:1991598/1‑150 (MQ=255)
                                                                                                   aTTCGCGTCGTCCCCACAGCGTCAGTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATGGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTCAAAAGAAACAgg                                                    >  2:2935491/1‑150 (MQ=255)
                                                                                                    ttCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGc                                                   >  1:3751332/1‑150 (MQ=255)
                                                                                                     tCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCa                                                  >  2:1843344/1‑150 (MQ=255)
                                                                                                      cgcgTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCaa                                                 <  1:2260889/150‑1 (MQ=255)
                                                                                                      cgcgTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCaa                                                 >  1:1765096/1‑150 (MQ=255)
                                                                                                      cgcgTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCaa                                                 <  1:2155863/150‑1 (MQ=255)
                                                                                                      cgcgTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCaa                                                 >  1:2384922/1‑150 (MQ=255)
                                                                                                      cgcgTCGTCCACACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCATATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCaa                                                 <  2:963783/150‑1 (MQ=255)
                                                                                                        cgtcgtCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACg                                               <  2:2803200/150‑1 (MQ=255)
                                                                                                           cgtcCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCg                                            <  2:4267036/150‑1 (MQ=255)
                                                                                                           cgtcCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCg                                            <  2:125895/150‑1 (MQ=255)
                                                                                                            gtcCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTCAGGGTAAAAAGAAACAGGCAACGGCGa                                           >  1:2602456/1‑150 (MQ=255)
                                                                                                              ccccACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAAt                                         >  2:1790725/1‑150 (MQ=255)
                                                                                                              ccccACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAAt                                         >  2:518916/1‑150 (MQ=255)
                                                                                                               cccACAGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATc                                        >  1:1968505/1‑150 (MQ=255)
                                                                                                                  acaGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCt                                     <  1:553966/150‑1 (MQ=255)
                                                                                                                    aGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGa                                   >  2:3384294/1‑150 (MQ=255)
                                                                                                                    aGCGTCATTGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGa                                   >  2:2262519/1‑150 (MQ=255)
                                                                                                                           ttGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACt                            <  2:3751332/150‑1 (MQ=255)
                                                                                                                           ttGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACt                            <  1:3353190/150‑1 (MQ=255)
                                                                                                                           ttGTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACt                            <  2:1690764/150‑1 (MQ=255)
                                                                                                                             gTCAACGACGTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGc                          <  2:3836635/150‑1 (MQ=255)
                                                                                                                                     cgTACCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAat                  <  2:3515312/150‑1 (MQ=255)
                                                                                                                                      gTACCAACTGATCACTAAAAGCTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGTTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAatg                 <  2:3799441/150‑1 (MQ=255)
                                                                                                                                        aCCAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAatgat               >  1:222114/1‑150 (MQ=255)
                                                                                                                                          cAACTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAATGATGc             <  1:2736829/150‑1 (MQ=255)
                                                                                                                                            aCTGATCCTTACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAATGATGCCt           <  1:1541254/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAAAGCGAACCATAATATTTATCCATAATTTGTGACGATGAATTAAACACGGGGTTAAATTGATCGATGGTTGCTGTCAACTGACTATGGTGTAACGATTCGCGTCGTCCCCACAGCGTCATTGTCAACGACGTACCAACTGATCCTGACAAGGTACGAAAAAAATTACTCATACTCGAGGCATTGGCAAATTTATTATCTGGCAAGCCTGAAAACGAAATCGTTGTTAAGGGTAAAAAGAAACAGGCAACGGCGAATCCCTGAAAAAACTGCGGCAAAATGATGCCTGTAAAATCA  >  CP009273/2474269‑2474567

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: