Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,168,773 A→T 36.8% A133A (GCA→GCT rplK → 50S ribosomal subunit protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,168,7730AT36.8% 266.4 / 298.3 247A133A (GCA→GCTrplK50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base A (77/79);  new base T (42/49);  total (119/128)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

TGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCG  >  CP009273/4168626‑4168921
                                                                                                                                                   |                                                                                                                                                    
  gTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGtt                                                                                                                                                  <  2:1931136/150‑1 (MQ=255)
  gTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGtt                                                                                                                                                  <  2:122853/150‑1 (MQ=255)
         gTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATggg                                                                                                                                           >  2:1473262/1‑150 (MQ=255)
           cTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGcc                                                                                                                                         >  1:4130227/1‑150 (MQ=255)
                tCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTa                                                                                                                                    <  1:506894/150‑1 (MQ=255)
                 ccGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAg                                                                                                                                   <  1:2321702/150‑1 (MQ=255)
                 ccGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAg                                                                                                                                   <  1:3691671/150‑1 (MQ=255)
                  cGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGt                                                                                                                                  >  2:1152762/1‑150 (MQ=255)
                  cGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGt                                                                                                                                  >  2:3792669/1‑150 (MQ=255)
                   ggTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTg                                                                                                                                 >  1:1934497/1‑150 (MQ=255)
                    gTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTgg                                                                                                                                <  2:3178864/150‑1 (MQ=255)
                    gTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTgg                                                                                                                                >  2:3230010/1‑150 (MQ=255)
                     tAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTgga                                                                                                                               <  1:2345066/150‑1 (MQ=255)
                      aaGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTggag                                                                                                                              >  2:2170657/1‑150 (MQ=255)
                       aGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTggagg                                                                                                                             <  1:1357547/150‑1 (MQ=255)
                           gAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTa                                                                                                                         >  1:400872/1‑150 (MQ=255)
                           gAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTa                                                                                                                         >  1:1669544/1‑150 (MQ=255)
                              caaagacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGa                                                                                                                      <  1:1244426/150‑1 (MQ=255)
                               aaagacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGaa                                                                                                                     <  2:3078924/150‑1 (MQ=255)
                               aaagacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGaa                                                                                                                     <  2:2081784/150‑1 (MQ=255)
                               aaagacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGaa                                                                                                                     <  1:3002465/150‑1 (MQ=255)
                                 agacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAAt                                                                                                                   >  1:189855/1‑150 (MQ=255)
                                  gacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATg                                                                                                                  <  1:2110823/150‑1 (MQ=255)
                                     aaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCt                                                                                                               <  2:170769/150‑1 (MQ=255)
                                        gTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTaaa                                                                                                            >  1:255002/1‑150 (MQ=255)
                                          gggTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACt                                                                                                          <  1:656041/150‑1 (MQ=255)
                                          gggTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACt                                                                                                          <  1:2728123/150‑1 (MQ=255)
                                          gggTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACt                                                                                                          <  1:362607/150‑1 (MQ=255)
                                            gTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGa                                                                                                        <  1:2371446/150‑1 (MQ=255)
                                             tAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGAc                                                                                                       <  1:2705598/150‑1 (MQ=255)
                                             tAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGAc                                                                                                       <  1:2705800/150‑1 (MQ=255)
                                             tAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGAc                                                                                                       <  1:2874714/150‑1 (MQ=255)
                                              aaaaTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGAcc                                                                                                      >  1:3583337/1‑150 (MQ=255)
                                              aaaaTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGAcc                                                                                                      >  2:4117273/1‑150 (MQ=255)
                                              aaaaTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGAcc                                                                                                      <  1:2315124/150‑1 (MQ=255)
                                                 aTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAg                                                                                                   <  1:1485372/150‑1 (MQ=255)
                                                 aTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAg                                                                                                   <  2:1174097/150‑1 (MQ=255)
                                                 aTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAg                                                                                                   <  2:786973/150‑1 (MQ=255)
                                                   ttCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGACATGGCTAAACTGACCAAgcg                                                                                                 >  2:1418355/1‑150 (MQ=255)
                                                    tCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAgcgc                                                                                                >  1:574650/1‑150 (MQ=255)
                                                    tCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAgcgc                                                                                                <  2:2406028/150‑1 (MQ=255)
                                                     cccGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCa                                                                                               >  2:846314/1‑150 (MQ=255)
                                                      ccGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCAt                                                                                              >  2:521441/1‑150 (MQ=255)
                                                         cgcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCg                                                                                           >  2:2803601/1‑150 (MQ=255)
                                                          gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCgt                                                                                          <  2:1210985/150‑1 (MQ=255)
                                                          gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCgt                                                                                          <  2:3303261/150‑1 (MQ=255)
                                                          gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCgt                                                                                          <  2:3427293/150‑1 (MQ=255)
                                                              aGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTa                                                                                      >  1:2209241/1‑150 (MQ=255)
                                                              aGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTa                                                                                      <  1:1264755/150‑1 (MQ=255)
                                                                cTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATc                                                                                    <  1:2204242/150‑1 (MQ=255)
                                                                 tGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATcc                                                                                   <  2:2703125/150‑1 (MQ=255)
                                                                 tGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATcc                                                                                   >  1:332745/1‑150 (MQ=255)
                                                                 tGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATcc                                                                                   <  1:873241/150‑1 (MQ=255)
                                                                        aaTCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGaa                                                                            <  1:716577/150‑1 (MQ=255)
                                                                         aTCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGaaa                                                                           <  1:2040732/150‑1 (MQ=255)
                                                                         aTCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGaaa                                                                           <  1:3062233/150‑1 (MQ=255)
                                                                           cgcgCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGt                                                                         <  1:3370052/150‑1 (MQ=255)
                                                                           cgcgCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGt                                                                         <  1:1313647/150‑1 (MQ=255)
                                                                             cgcAGACCAAAGCTGCCGACATGACTGGTGCCGNCATTGAAGCNATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTg                                                                       <  1:2946932/150‑1 (MQ=255)
                                                                             cgcAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTg                                                                       <  1:4155960/150‑1 (MQ=255)
                                                                             cgcAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTg                                                                       <  2:1995706/150‑1 (MQ=255)
                                                                              gcAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGa                                                                      >  2:1685021/1‑150 (MQ=255)
                                                                              gcAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGa                                                                      <  1:2766891/150‑1 (MQ=255)
                                                                                aGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATg                                                                    >  2:4044610/1‑150 (MQ=255)
                                                                                 gACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGc                                                                   >  2:1274422/1‑150 (MQ=255)
                                                                                 gACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGc                                                                   >  1:1754403/1‑150 (MQ=255)
                                                                                  aCCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCa                                                                  >  1:779535/1‑150 (MQ=255)
                                                                                      aaGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCa                                                              >  1:3328909/1‑150 (MQ=255)
                                                                                         cTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAAc                                                           <  1:2309446/150‑1 (MQ=255)
                                                                                         cTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAAc                                                           >  1:1115732/1‑150 (MQ=255)
                                                                                         cTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAAc                                                           <  1:3414471/150‑1 (MQ=255)
                                                                                          tGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACa                                                          >  1:1406462/1‑150 (MQ=255)
                                                                                             cGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTa                                                       >  2:4022404/1‑150 (MQ=255)
                                                                                               aCATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACg                                                     >  2:2283463/1‑150 (MQ=255)
                                                                                                    aCTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATc                                                >  1:3210667/1‑150 (MQ=255)
                                                                                                      tGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCaa                                              >  2:2918936/1‑150 (MQ=255)
                                                                                                       ggTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAAc                                             <  1:2403185/150‑1 (MQ=255)
                                                                                                       ggTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGGTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAAc                                             <  2:1577232/150‑1 (MQ=255)
                                                                                                             gACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCCAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCt                                       <  2:2217337/150‑1 (MQ=255)
                                                                                                                aTTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATc                                    <  1:3499916/150‑1 (MQ=255)
                                                                                                                aTTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATc                                    <  1:3258784/150‑1 (MQ=255)
                                                                                                                         aTGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCActgctg                           <  1:1028659/150‑1 (MQ=255)
                                                                                                                          tGACTCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGa                          <  1:1127157/150‑1 (MQ=255)
                                                                                                                             cTCGCTCCATCGAAGGTACTGCTAGTTCCAAGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCCAACAGTACGACATCAACTAAGCTATCGCACTGCTGAAag                       >  2:1454699/1‑150 (MQ=255)
                                                                                                                              tCGCTCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAaga                      >  1:3048304/1‑150 (MQ=255)
                                                                                                                              tCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAaga                      <  2:2955708/150‑1 (MQ=255)
                                                                                                                                gCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGc                    >  1:3962582/1‑150 (MQ=255)
                                                                                                                                 cTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCt                   <  1:3152408/150‑1 (MQ=255)
                                                                                                                                  tCCATCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTg                  <  2:4102273/150‑1 (MQ=255)
                                                                                                                                   ccATGGTAGGTTCTGCTCGTTCCATGGGCCTGGTAGTGGATGACTTAGAAATGGCTAAACTGACCAAGCGCATGAGTGTTATCCGCGTGAAAGTTGATGCATCCATTCTGTACGACTTCAACGAAGCTATCGCACTGCTGAAAGAGCTgg                 <  2:1942595/150‑1 (MQ=255)
                                                                                                                                    cATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGc                >  1:3776184/1‑150 (MQ=255)
                                                                                                                                     aTCGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCg               >  2:2320822/1‑150 (MQ=255)
                                                                                                                                       cGAAGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGAc             <  1:3266237/150‑1 (MQ=255)
                                                                                                                                        gAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACt            >  1:3815237/1‑150 (MQ=255)
                                                                                                                                        gAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACt            >  2:524128/1‑150 (MQ=255)
                                                                                                                                         aaGGTACTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTg           <  1:2605924/150‑1 (MQ=255)
                                                                                                                                            gTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTa        >  1:423531/1‑150 (MQ=255)
                                                                                                                                               cTGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAAt     <  2:2417538/150‑1 (MQ=255)
                                                                                                                                                tGCTCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAAtt    >  1:1609462/1‑150 (MQ=255)
                                                                                                                                                  cACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCg  >  2:383835/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCG  >  CP009273/4168626‑4168921

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: