Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 961,946 T→A 8.7% F724Y (TTT→TAT)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273961,9460TA 8.7% 1331.0 / 92.0 392F724Y (TTT→TAT) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (177/181);  new base A (14/20);  total (191/201)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 426 total aligned reads displayed.

AATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATC  >  CP009273/961797‑962095
                                                                                                                                                     |                                                                                                                                                     
aaTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAttt                                                                                                                                                       <  1:2828416/150‑1 (MQ=255)
aaTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAttt                                                                                                                                                       <  2:241523/150‑1 (MQ=255)
aaTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAttt                                                                                                                                                       >  2:3576414/1‑150 (MQ=255)
aaTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAttt                                                                                                                                                       >  1:2410985/1‑150 (MQ=255)
 aTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAATCAACGGTCAAtttt                                                                                                                                                      >  2:2712439/1‑150 (MQ=255)
  ttCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAattatt                                                                                                                                                     <  1:2714618/150‑1 (MQ=255)
    cAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCa                                                                                                                                                   <  2:216070/150‑1 (MQ=255)
      gCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAgc                                                                                                                                                 <  1:2108597/150‑1 (MQ=255)
       cGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAgcg                                                                                                                                                >  1:966818/1‑150 (MQ=255)
        gagagaATGGAAAAGTGGCACTCGCATAGGCATCGCGCTATAACGCATGGCGACTGCCCTCAAACAAAGTAAAGCATCGCTATCAACTGCAAGTATATCAATGGATTGATACACCACATCAACGACAAACAACGGTCAATTTTTCAgcgc                                                                                                                                               >  2:3168551/1‑150 (MQ=255)
          gTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCaa                                                                                                                                             >  1:4756145/1‑150 (MQ=255)
             aaTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGc                                                                                                                                          >  1:3453729/1‑150 (MQ=255)
              aTGGAAACGTGGCTCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCACTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCTAACAACGGTCAATTTTTCAGCGCAAGGCt                                                                                                                                         <  2:5265049/150‑1 (MQ=255)
                 gAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGc                                                                                                                                      >  1:5161230/1‑150 (MQ=255)
                  aaaaGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCg                                                                                                                                     <  2:2833735/150‑1 (MQ=255)
                     aGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGAt                                                                                                                                  >  1:3186222/1‑150 (MQ=255)
                       tGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAATGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTa                                                                                                                                <  2:5572391/150‑1 (MQ=255)
                       tGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTa                                                                                                                                <  2:5442398/150‑1 (MQ=255)
                        ggCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTag                                                                                                                               <  2:3386392/150‑1 (MQ=255)
                         gCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCATCTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTagc                                                                                                                              <  2:4289341/150‑1 (MQ=255)
                         gCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTagc                                                                                                                              <  2:3108107/150‑1 (MQ=255)
                           aCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTagcag                                                                                                                            >  1:1442318/1‑150 (MQ=255)
                                 atcggcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAg                                                                                                                      >  1:3580690/1‑150 (MQ=255)
                                    ggcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGa                                                                                                                   <  2:3768752/150‑1 (MQ=255)
                                     gcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAg                                                                                                                  <  2:818673/150‑1 (MQ=255)
                                     gcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAg                                                                                                                  >  1:1853629/1‑150 (MQ=255)
                                      catcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGc                                                                                                                 >  1:3567964/1‑150 (MQ=255)
                                         cgcgCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCaaa                                                                                                              >  2:586009/1‑150 (MQ=255)
                                           cgcTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAAtt                                                                                                            <  1:290044/150‑1 (MQ=255)
                                                taACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTAcc                                                                                                       <  1:574229/150‑1 (MQ=255)
                                                 aaCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCt                                                                                                      >  1:1318027/1‑150 (MQ=255)
                                                 aaCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCt                                                                                                      <  1:2151494/150‑1 (MQ=255)
                                                 aaCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCt                                                                                                      >  2:4899132/1‑150 (MQ=255)
                                                 aaCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCt                                                                                                      >  2:4112004/1‑150 (MQ=255)
                                                  aCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCATTTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTc                                                                                                     <  2:5394216/150‑1 (MQ=255)
                                                  aCGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTc                                                                                                     <  1:4897083/150‑1 (MQ=255)
                                                   cGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCg                                                                                                    <  1:1772284/150‑1 (MQ=255)
                                                          cGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTat                                                                                             >  1:5592581/1‑150 (MQ=255)
                                                              tGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTatcatc                                                                                         <  1:1000478/150‑1 (MQ=255)
                                                               gCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTAtcatca                                                                                        <  2:2361984/150‑1 (MQ=255)
                                                                  ctctAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTg                                                                                     >  2:3916227/1‑150 (MQ=255)
                                                                    ctAACAAAGTTAAGCATCGCTATCAACGGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGt                                                                                   >  1:2593713/1‑150 (MQ=255)
                                                                         aaaGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGc                                                                              <  2:4033329/150‑1 (MQ=255)
                                                                                aGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAg                                                                       <  1:4025769/150‑1 (MQ=255)
                                                                                 gCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGt                                                                      <  1:5301244/150‑1 (MQ=255)
                                                                                    tCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGa                                                                   <  1:1080804/150‑1 (MQ=255)
                                                                                     cGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGAt                                                                  >  1:5324793/1‑150 (MQ=255)
                                                                                      gCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATa                                                                 <  2:4405558/150‑1 (MQ=255)
                                                                                      gCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATa                                                                 >  2:5278908/1‑150 (MQ=255)
                                                                                      gCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATa                                                                 <  1:5433938/150‑1 (MQ=255)
                                                                                         aTCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACg                                                              >  1:1615749/1‑150 (MQ=255)
                                                                                          tCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCAGGCGGATTAGCAGCCACAGGGAGCAAATTTTACCACGTTGGTATCATCAGTGGGATGGCGTGCCAGTGGATAACgg                                                             >  2:1635162/1‑150 (MQ=255)
                                                                                           cAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACggg                                                            >  2:4170858/1‑150 (MQ=255)
                                                                                           cAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACggg                                                            >  1:2231133/1‑150 (MQ=255)
                                                                                           cAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACggg                                                            <  1:665482/150‑1 (MQ=255)
                                                                                           cAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACggg                                                            >  1:4499370/1‑150 (MQ=255)
                                                                                            aaCTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGt                                                           >  2:1351581/1‑150 (MQ=255)
                                                                                            aaCTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGt                                                           >  2:5069569/1‑150 (MQ=255)
                                                                                             aCTGCAAGGATATCAATGGATTGATACTCCACATCCAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTa                                                          <  1:2648689/150‑1 (MQ=255)
                                                                                              cTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAg                                                         <  1:2809860/150‑1 (MQ=255)
                                                                                               tGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGa                                                        <  2:5132394/150‑1 (MQ=255)
                                                                                               tGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGa                                                        <  2:2662408/150‑1 (MQ=255)
                                                                                               tGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGa                                                        >  2:5589096/1‑150 (MQ=255)
                                                                                                gCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGaa                                                       <  1:2025714/150‑1 (MQ=255)
                                                                                                 cAAGTATATCAAGTGAGTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAat                                                      <  2:3966412/150‑1 (MQ=255)
                                                                                                  aaGTATATCAATGTAATGATACTCCACATCAAGGTCAAACAACGGTCAATTATTAAGCGCAAGGCTGGCGGATTAGCAGCCTCATGGAGCAAATTTTACCTCGTTTGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAata                                                     <  2:5531325/150‑1 (MQ=255)
                                                                                                  aaGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAata                                                     <  1:3262630/150‑1 (MQ=255)
                                                                                                   aGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAatat                                                    >  2:940801/1‑150 (MQ=255)
                                                                                                      atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCg                                                 <  2:294263/150‑1 (MQ=255)
                                                                                                      atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCg                                                 <  1:879981/150‑1 (MQ=255)
                                                                                                      atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCg                                                 >  2:1539635/1‑150 (MQ=255)
                                                                                                        atCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGc                                               <  2:1947919/150‑1 (MQ=255)
                                                                                                            aTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGTATATGCGGCTAtt                                           >  2:3859521/1‑150 (MQ=255)
                                                                                                             tGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAttt                                          >  1:2319979/1‑150 (MQ=255)
                                                                                                             tGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAttt                                          <  2:4495605/150‑1 (MQ=255)
                                                                                                             tGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAttt                                          <  2:1841204/150‑1 (MQ=255)
                                                                                                               gcttGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTca                                        <  1:991106/148‑1 (MQ=255)
                                                                                                               gATTGATACTCCACCTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGAGTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTca                                        <  2:2116005/150‑1 (MQ=255)
                                                                                                                 ttGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaac                                      <  2:5161230/150‑1 (MQ=255)
                                                                                                                   gATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaacaa                                    <  1:4405234/150‑1 (MQ=255)
                                                                                                                   gATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaacaa                                    <  2:2234335/150‑1 (MQ=255)
                                                                                                                     tACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAAt                                  <  1:2996015/150‑1 (MQ=255)
                                                                                                                        tCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCt                               <  1:4359458/150‑1 (MQ=255)
                                                                                                                          cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTgg                             <  1:774800/150‑1 (MQ=255)
                                                                                                                          cacaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTgg                             <  1:4115244/150‑1 (MQ=255)
                                                                                                                           acaTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGt                            >  2:4352984/1‑150 (MQ=255)
                                                                                                                             aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTAGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGttt                          >  2:1263950/1‑150 (MQ=255)
                                                                                                                              tCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGtttt                         >  1:3959394/1‑150 (MQ=255)
                                                                                                                                  ggTCAAACAGCGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                     >  1:3797803/1‑150 (MQ=255)
                                                                                                                                  ggTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                     <  1:3168551/150‑1 (MQ=255)
                                                                                                                                  ggTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                     >  2:1439565/1‑150 (MQ=255)
                                                                                                                                   gTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACACATGCTGGTTTTTTGaa                    >  2:5416525/1‑150 (MQ=255)
                                                                                                                                         caacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATa              <  2:4704057/150‑1 (MQ=255)
                                                                                                                                           acGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAAc            <  1:5106/150‑1 (MQ=255)
                                                                                                                                              gTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGAc         >  2:3739324/1‑150 (MQ=255)
                                                                                                                                                  attattCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAg     >  1:4650778/1‑150 (MQ=255)
                                                                                                                                                  aTTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAg     >  1:4298364/1‑150 (MQ=255)
                                                                                                                                                   tttttCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGa    <  1:483557/150‑1 (MQ=255)
                                                                                                                                                     tttCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGAtc  <  1:3140041/150‑1 (MQ=255)
                                                                                                                                                     tttCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGAtc  <  2:2259875/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAGATC  >  CP009273/961797‑962095

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: