Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 961,946 T→A 26.5% 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,9460TA26.5% 484.3 / 216.8 249F724Y (TTT→TAT) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (81/102);  new base A (25/41);  total (106/143)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 254 total aligned reads displayed.

ATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAG  >  CP009273/961798‑962092
                                                                                                                                                    |                                                                                                                                                  
aTTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAtttt                                                                                                                                                   >  1:3019412/1‑150 (MQ=255)
 ttCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAAttttt                                                                                                                                                  <  1:3788538/150‑1 (MQ=255)
   cAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCa                                                                                                                                                <  2:1877718/150‑1 (MQ=255)
   cAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCa                                                                                                                                                >  1:713214/1‑150 (MQ=255)
   cAGCGAGTGAATGGAAAAGTGGCACTCGCATAGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCa                                                                                                                                                >  1:2469966/1‑150 (MQ=255)
    aGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAg                                                                                                                                               >  2:1269908/1‑150 (MQ=255)
     gCGAGTGAATTGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAgc                                                                                                                                              <  1:4639353/150‑1 (MQ=255)
       gAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAgcgc                                                                                                                                            >  1:227810/1‑150 (MQ=255)
        aGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCa                                                                                                                                           <  1:1288587/150‑1 (MQ=255)
           gAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAgg                                                                                                                                        >  2:1535260/1‑150 (MQ=255)
             aTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCt                                                                                                                                      <  1:1332166/150‑1 (MQ=255)
             aTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCt                                                                                                                                      <  1:3480705/150‑1 (MQ=255)
             aTGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCt                                                                                                                                      <  2:3538550/150‑1 (MQ=255)
              tGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTg                                                                                                                                     <  2:910020/150‑1 (MQ=255)
                gAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGc                                                                                                                                   <  1:1714011/150‑1 (MQ=255)
                 aaaaGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCg                                                                                                                                  >  1:659937/1‑150 (MQ=255)
                     gTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAAGGGTCAATTTTTCAGCGCAAGGCTGGCGGAtt                                                                                                                              <  1:76834/150‑1 (MQ=255)
                      tGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTa                                                                                                                             >  2:2193049/1‑150 (MQ=255)
                          aCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTagcag                                                                                                                         <  1:1489906/150‑1 (MQ=255)
                            tCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGcc                                                                                                                       <  1:1014008/150‑1 (MQ=255)
                             cGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCt                                                                                                                      <  2:2736778/150‑1 (MQ=255)
                              gcatcggcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTc                                                                                                                     >  2:1742965/1‑150 (MQ=255)
                                 tcggcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAgg                                                                                                                  >  2:705823/1‑150 (MQ=255)
                                    gcatcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAg                                                                                                               <  1:1398451/150‑1 (MQ=255)
                                     catcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGc                                                                                                              <  2:4231426/150‑1 (MQ=255)
                                      atcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCa                                                                                                             >  2:1339414/1‑150 (MQ=255)
                                      atcgCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAACTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCa                                                                                                             <  1:235357/150‑1 (MQ=255)
                                        cgcgCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCaaa                                                                                                           >  2:2356498/1‑150 (MQ=255)
                                           gcTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAAttt                                                                                                        <  2:3196666/150‑1 (MQ=255)
                                            cTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAAtttt                                                                                                       >  2:2111495/1‑150 (MQ=255)
                                             tataACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTa                                                                                                      >  2:3063136/1‑150 (MQ=255)
                                              ataACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTAc                                                                                                     >  1:2021548/1‑150 (MQ=255)
                                                  cGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCg                                                                                                 >  1:1384682/1‑150 (MQ=255)
                                                      tGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTgg                                                                                             <  2:4853384/150‑1 (MQ=255)
                                                        gCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTa                                                                                           <  2:140477/150‑1 (MQ=255)
                                                        gCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTa                                                                                           <  2:3498746/150‑1 (MQ=255)
                                                           aCTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTatca                                                                                        <  1:3546423/150‑1 (MQ=255)
                                                           aCTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATCCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTatca                                                                                        <  2:2521253/150‑1 (MQ=255)
                                                            cTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTatcat                                                                                       <  1:4313291/150‑1 (MQ=255)
                                                            cTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTatcat                                                                                       <  1:4313397/150‑1 (MQ=255)
                                                               cccTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAg                                                                                    >  2:2387950/1‑150 (MQ=255)
                                                               cccTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAg                                                                                    >  2:1251511/1‑150 (MQ=255)
                                                                 ctctAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTg                                                                                  >  2:2653930/1‑150 (MQ=255)
                                                                   ctAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGt                                                                                <  1:2730778/150‑1 (MQ=255)
                                                                       cAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTgg                                                                            <  2:2865565/150‑1 (MQ=255)
                                                                       cAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTgg                                                                            >  2:1824966/1‑150 (MQ=255)
                                                                        aaaGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGc                                                                           <  2:761500/150‑1 (MQ=255)
                                                                         aaGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCg                                                                          <  1:1646225/150‑1 (MQ=255)
                                                                          aGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGt                                                                         >  1:669430/1‑150 (MQ=255)
                                                                                     gCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATa                                                              >  1:4532435/1‑150 (MQ=255)
                                                                                       tATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAc                                                            >  2:3551240/1‑150 (MQ=255)
                                                                                       tATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAAc                                                            <  2:3019412/150‑1 (MQ=255)
                                                                                        aTCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACg                                                           >  2:3932238/1‑150 (MQ=255)
                                                                                         tCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACgg                                                          <  1:1900698/150‑1 (MQ=255)
                                                                                           aaCTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGt                                                        >  2:3415470/1‑150 (MQ=255)
                                                                                           aaCTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGt                                                        <  1:3258661/150‑1 (MQ=255)
                                                                                             cTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAg                                                      <  2:3493971/150‑1 (MQ=255)
                                                                                              tGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGa                                                     <  2:1942026/150‑1 (MQ=255)
                                                                                                cAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAat                                                   >  1:4704021/1‑150 (MQ=255)
                                                                                                cAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAat                                                   <  2:1701758/150‑1 (MQ=255)
                                                                                                  aGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAatat                                                 <  2:2213511/150‑1 (MQ=255)
                                                                                                     atatCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCg                                              <  1:2239885/150‑1 (MQ=255)
                                                                                                      tatCAATGGATTGATACTCCACATCAATGTCTAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCgg                                             <  2:2891820/150‑1 (MQ=255)
                                                                                                          aaTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAt                                         <  2:3597309/150‑1 (MQ=255)
                                                                                                          aaTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAt                                         >  2:3825098/1‑150 (MQ=255)
                                                                                                           aTGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAtt                                        >  2:2739915/1‑150 (MQ=255)
                                                                                                            tGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTAttt                                       >  2:3529293/1‑150 (MQ=255)
                                                                                                                 tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaaca                                  <  2:3732816/150‑1 (MQ=255)
                                                                                                                 tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaaca                                  <  1:4789374/150‑1 (MQ=255)
                                                                                                                 tGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaaca                                  <  2:115584/150‑1 (MQ=255)
                                                                                                                 tGATACTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTcaaca                                  <  1:2156367/150‑1 (MQ=255)
                                                                                                                    tACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAAt                               >  2:492202/1‑150 (MQ=255)
                                                                                                                      cTCCACATCAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGc                             <  1:1157578/150‑1 (MQ=255)
                                                                                                                        ccACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTg                           <  1:1596119/150‑1 (MQ=255)
                                                                                                                           caTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTCTCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGtt                        >  2:1943824/1‑150 (MQ=255)
                                                                                                                            aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGttt                       <  2:3212328/150‑1 (MQ=255)
                                                                                                                            aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGttt                       >  1:1259267/1‑150 (MQ=255)
                                                                                                                            aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGttt                       <  2:4626483/150‑1 (MQ=255)
                                                                                                                            aTCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGttt                       <  2:4381540/150‑1 (MQ=255)
                                                                                                                             tCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGtttt                      <  1:1520535/150‑1 (MQ=255)
                                                                                                                              cAAGGTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGttttt                     >  1:119082/1‑150 (MQ=255)
                                                                                                                               aaGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGtttttt                    <  2:1389967/150‑1 (MQ=255)
                                                                                                                               aaGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGtttttt                    <  2:452828/150‑1 (MQ=255)
                                                                                                                                 ggTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                  >  2:3787033/1‑150 (MQ=255)
                                                                                                                                 ggTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                  <  1:2597400/150‑1 (MQ=255)
                                                                                                                                 ggTCAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGa                  <  2:3996169/150‑1 (MQ=255)
                                                                                                                                    cAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATg               >  2:4284864/1‑150 (MQ=255)
                                                                                                                                    cAAACAACGGTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATg               >  2:4284771/1‑150 (MQ=255)
                                                                                                                                     aaaCAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACACAGGCTGGTTTTTTGGATGc              >  1:3483175/1‑150 (MQ=255)
                                                                                                                                     aaaCAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGc              >  1:1815968/1‑150 (MQ=255)
                                                                                                                                       acaacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCAt            <  1:3732973/150‑1 (MQ=255)
                                                                                                                                        caacGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATa           >  1:956010/1‑150 (MQ=255)
                                                                                                                                          acGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAAc         >  2:4141571/1‑150 (MQ=255)
                                                                                                                                          acGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAAc         <  1:341189/150‑1 (MQ=255)
                                                                                                                                             gTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGAc      <  1:3136462/150‑1 (MQ=255)
                                                                                                                                             gTCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGAc      <  1:891036/150‑1 (MQ=255)
                                                                                                                                              tCAATTATTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACa     <  2:1056167/150‑1 (MQ=255)
                                                                                                                                               cAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaa    >  2:1766556/1‑150 (MQ=255)
                                                                                                                                                aaTTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACaaa   >  1:519052/1‑150 (MQ=255)
                                                                                                                                                 aTTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAg  <  2:648157/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCCACATCAAGGTCAAACAACGGTCAATTTTTCAGCGCAAGGCTGGCGGATTAGCAGCCTCAGGGAGCAAATTTTACCTCGTTGGTATCATCAGTGGTTTGGCGTGCCAGTGGATAACGGGTAGAATATGCGGCTATTTCAACAAATGCTGGTTTTTTGAATGCATAACGACAAAG  >  CP009273/961798‑962092

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: