Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,351,022 G→T 100% R116S (CGT→AGT)  sapA ← antimicrobial peptide transport ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,351,0220GT97.6% 782.4 / 6.7 255R116S (CGT→AGT) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (124/125);  total (126/129)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 256 total aligned reads displayed.

ATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCAT  >  CP009273/1350876‑1351171
                                                                                                                                                  |                                                                                                                                                     
aTCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGt                                                                                                                                                    <  2:3718746/150‑1 (MQ=255)
     gTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAg                                                                                                                                               >  2:1170949/1‑150 (MQ=255)
       ttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTa                                                                                                                                             >  1:5144099/1‑150 (MQ=255)
         gCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTaaa                                                                                                                                           >  1:4337024/1‑150 (MQ=255)
          cGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAc                                                                                                                                          >  1:4000362/1‑150 (MQ=255)
           ggACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAcc                                                                                                                                         >  1:2776169/1‑150 (MQ=255)
              ctcTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGt                                                                                                                                      >  2:5143328/3‑150 (MQ=255)
               tcTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTc                                                                                                                                     <  1:1352968/149‑1 (MQ=255)
                cTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCg                                                                                                                                    <  1:3420998/150‑1 (MQ=255)
                    tCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGttt                                                                                                                                >  2:4499561/1‑150 (MQ=255)
                      aCGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGttttt                                                                                                                              <  2:4158990/150‑1 (MQ=255)
                      aCGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGttttt                                                                                                                              <  2:3448730/150‑1 (MQ=255)
                      aCGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGttttt                                                                                                                              <  2:1551978/150‑1 (MQ=255)
                      aCGTTATCGGCAAATTGCAGGCTGTCGAAGGAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGttttt                                                                                                                              <  2:1383098/150‑1 (MQ=255)
                       cGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGtttttt                                                                                                                             >  1:287511/1‑150 (MQ=255)
                         ttATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGa                                                                                                                           >  1:1088397/1‑150 (MQ=255)
                             cGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACg                                                                                                                       <  2:2377774/150‑1 (MQ=255)
                                 aaaTTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAAc                                                                                                                   >  1:1721834/1‑150 (MQ=255)
                                     tGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATcg                                                                                                               <  2:4100541/150‑1 (MQ=255)
                                       cAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATcgcg                                                                                                             >  2:1937622/1‑150 (MQ=255)
                                        aGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGa                                                                                                            >  1:330996/1‑150 (MQ=255)
                                         ggCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGAc                                                                                                           <  1:4399635/150‑1 (MQ=255)
                                         ggCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGAc                                                                                                           >  2:857989/1‑150 (MQ=255)
                                          gCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACg                                                                                                          <  2:4258887/150‑1 (MQ=255)
                                               cGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGt                                                                                                     <  1:1028863/150‑1 (MQ=255)
                                                gAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGAc                                                                                                           <  1:2049292/143‑1 (MQ=255)
                                                 aaGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTgg                                                                                                   >  1:2073164/1‑150 (MQ=255)
                                                  aGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGa                                                                                                  >  2:2422135/1‑150 (MQ=255)
                                                  aGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGa                                                                                                  >  2:562901/1‑150 (MQ=255)
                                                   gTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGaa                                                                                                 >  1:1243794/1‑150 (MQ=255)
                                                    tAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAg                                                                                                >  2:3662226/1‑150 (MQ=255)
                                                      tgggAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCg                                                                                              >  2:4934733/2‑150 (MQ=255)
                                                      ggggAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCg                                                                                              <  1:197605/150‑1 (MQ=255)
                                                         gAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATa                                                                                           <  1:656048/150‑1 (MQ=255)
                                                          aaGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAg                                                                                          <  1:300778/150‑1 (MQ=255)
                                                              tgctgcCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTcg                                                                                      <  2:1441512/150‑1 (MQ=255)
                                                              tgctgcCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTcg                                                                                      >  2:2279649/1‑150 (MQ=255)
                                                                ctgcCGTGGACGTTATGCCACTGGTTGTTGCGGTCAAAAATTCGCTGGAATGTAAACACCACATCGTCGGCATTAATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTcgcg                                                                                    <  2:4151477/150‑1 (MQ=255)
                                                                  gcCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGcc                                                                                  <  1:1511137/150‑1 (MQ=255)
                                                                   ccGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCg                                                                                 >  2:1829614/1‑150 (MQ=255)
                                                                     gTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGtt                                                                               <  1:762238/150‑1 (MQ=255)
                                                                     gTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGtt                                                                               <  1:1677287/150‑1 (MQ=255)
                                                                        gACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTc                                                                            <  1:471331/150‑1 (MQ=255)
                                                                        gACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTc                                                                            <  1:989009/150‑1 (MQ=255)
                                                                          cGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGa                                                                          >  1:888152/1‑150 (MQ=255)
                                                                          cGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGa                                                                          <  2:1691882/150‑1 (MQ=255)
                                                                            ttATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGt                                                                        >  2:2062434/1‑150 (MQ=255)
                                                                              aTGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTAc                                                                      >  2:3297341/1‑150 (MQ=255)
                                                                               tGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGTATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACt                                                                     <  1:4597817/150‑1 (MQ=255)
                                                                                 ccACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTc                                                                   <  1:3726521/150‑1 (MQ=255)
                                                                                  cACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTcc                                                                  >  2:4895316/1‑150 (MQ=255)
                                                                                  cACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTcc                                                                  <  2:3901608/150‑1 (MQ=255)
                                                                                  cACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTcc                                                                  <  2:3868539/150‑1 (MQ=255)
                                                                                  cACGGGTTGTTGCGGTCAAAAATTCGCGGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTcc                                                                  >  2:447030/1‑150 (MQ=255)
                                                                                    cGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCa                                                                <  1:2612381/150‑1 (MQ=255)
                                                                                    cGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCa                                                                <  2:795200/150‑1 (MQ=255)
                                                                                       gttgttGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCt                                                             <  2:1448740/150‑1 (MQ=255)
                                                                                        ttgttgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCtt                                                            <  2:1292155/150‑1 (MQ=255)
                                                                                            tgCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCgg                                                        <  1:850043/150‑1 (MQ=255)
                                                                                               ggTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCaa                                                     <  1:177337/150‑1 (MQ=255)
                                                                                               ggTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCaa                                                     <  2:3786130/150‑1 (MQ=255)
                                                                                               ggTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCaa                                                     >  2:28058/1‑150 (MQ=255)
                                                                                               ggTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCaa                                                     <  2:3580640/150‑1 (MQ=255)
                                                                                                gTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAg                                                    <  1:5147635/150‑1 (MQ=255)
                                                                                                    aaaaTTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTcc                                                <  1:4881567/150‑1 (MQ=255)
                                                                                                        ttCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCa                                            <  2:3786032/150‑1 (MQ=255)
                                                                                                        ttCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCa                                            <  1:5013270/150‑1 (MQ=255)
                                                                                                         tCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCAt                                           <  1:3521706/150‑1 (MQ=255)
                                                                                                         tCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCAt                                           <  1:513496/150‑1 (MQ=255)
                                                                                                            cTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAg                                        <  1:3175002/150‑1 (MQ=255)
                                                                                                             tGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAgg                                       <  2:2880096/150‑1 (MQ=255)
                                                                                                              ggAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGc                                      >  1:131932/1‑150 (MQ=255)
                                                                                                              ggAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCGTTCGGCAAGTTCCGGCCTCAGGc                                      >  2:4349132/1‑150 (MQ=255)
                                                                                                                 atgTAAACACCACATCGTCGGAATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAAAGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTACCATCTTTCGGCAAGTTCCGGCATCAGGCGAt                                   >  2:4723115/3‑150 (MQ=255)
                                                                                                                  ggTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATa                                  >  2:4230793/1‑150 (MQ=255)
                                                                                                                    tAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAgg                                <  2:5056890/150‑1 (MQ=255)
                                                                                                                       acacCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGtat                             <  2:3757249/150‑1 (MQ=255)
                                                                                                                       acacCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGtat                             <  1:3180942/150‑1 (MQ=255)
                                                                                                                        caccacATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGtata                            <  1:2598536/150‑1 (MQ=255)
                                                                                                                        caccacATCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGtata                            >  1:3165516/1‑150 (MQ=255)
                                                                                                                              aTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATc                      >  1:3874491/1‑150 (MQ=255)
                                                                                                                              aTCGTCGGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGGTCCGGCATCAGGCGATAGGTATAGGGATc                      <  2:3122464/150‑1 (MQ=255)
                                                                                                                               tcgtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCg                     <  1:2279651/150‑1 (MQ=255)
                                                                                                                                cgtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATAGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGa                    >  1:2140358/1‑150 (MQ=255)
                                                                                                                                 gtcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGAc                   >  2:2595865/1‑150 (MQ=255)
                                                                                                                                  tcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACa                  >  2:2947352/1‑150 (MQ=255)
                                                                                                                                  tcgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACa                  >  2:2811178/1‑150 (MQ=255)
                                                                                                                                   cgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACAt                 >  2:4555099/1‑150 (MQ=255)
                                                                                                                                   cgGCATTCATTTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACAt                 >  2:2153772/1‑150 (MQ=255)
                                                                                                                                      cattcattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATccc              >  1:3384781/1‑149 (MQ=255)
                                                                                                                                       attcattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcag             <  2:3196075/150‑1 (MQ=255)
                                                                                                                                         tcattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagca           <  2:4715385/150‑1 (MQ=255)
                                                                                                                                         tcattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagca           <  2:2008967/150‑1 (MQ=255)
                                                                                                                                          cattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagcag          >  2:2785876/1‑150 (MQ=255)
                                                                                                                                          cattTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagcag          <  2:692997/150‑1 (MQ=255)
                                                                                                                                           attTTACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGt         >  1:706771/1‑150 (MQ=255)
                                                                                                                                            ttttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTc        <  1:3221279/150‑1 (MQ=255)
                                                                                                                                            ttttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTc        >  2:5085835/1‑150 (MQ=255)
                                                                                                                                              ttACTAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGa      <  1:4126902/150‑1 (MQ=255)
                                                                                                                                                  gAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCAt  >  2:3983394/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCAT  >  CP009273/1350876‑1351171

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: