Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,957,480 C→A 8.7% L19F (TTG→TTT ppdA ← putative prepilin peptidase‑dependent protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,957,4800CA 8.7% 640.5 / 24.6 183L19F (TTG→TTTppdAputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base C (53/114);  new base A (9/7);  total (62/121)
Fisher's exact test for biased strand distribution p-value = 5.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 184 total aligned reads displayed.

AACGCTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAG  >  CP009273/2957333‑2957623
                                                                                                                                                   |                                                                                                                                               
aaCGCTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGAcc                                                                                                                                                 >  2:2830496/1‑148 (MQ=255)
   gCTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCaaa                                                                                                                                              <  2:1536665/148‑1 (MQ=255)
    cTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCaaaa                                                                                                                                             >  2:2753570/1‑148 (MQ=255)
    cTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTAGTGCCAGCCATCGAGTCCACGTGCGCTTATCATGACCaaaa                                                                                                                                             <  2:2970581/148‑1 (MQ=255)
    cTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAAGAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCaaaa                                                                                                                                             <  2:3128092/148‑1 (MQ=255)
     tGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAAt                                                                                                                                            <  2:1288380/148‑1 (MQ=255)
      gATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATc                                                                                                                                           <  1:3144652/148‑1 (MQ=255)
      gATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATc                                                                                                                                           <  2:2417167/148‑1 (MQ=255)
       aTACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCa                                                                                                                                          <  1:2785481/148‑1 (MQ=255)
           tgtgGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCAt                                                                                                                                      <  2:3171244/148‑1 (MQ=255)
             tgGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATcg                                                                                                                                    >  2:3084415/1‑148 (MQ=255)
             tgGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATcg                                                                                                                                    >  2:1444360/1‑148 (MQ=255)
               gTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATcgcg                                                                                                                                  >  2:2748583/1‑148 (MQ=255)
                   cgGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCa                                                                                                                              >  1:1281623/1‑148 (MQ=255)
                   cgGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCa                                                                                                                              >  1:2937819/1‑148 (MQ=255)
                     gTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGc                                                                                                                            <  1:173746/148‑1 (MQ=255)
                      ttATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCg                                                                                                                           <  2:2558429/148‑1 (MQ=255)
                      ttATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCg                                                                                                                           <  1:2485627/148‑1 (MQ=255)
                      ttATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCg                                                                                                                           <  1:2483514/148‑1 (MQ=255)
                      ttATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCg                                                                                                                           >  1:1332234/1‑148 (MQ=255)
                         tGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGttt                                                                                                                        <  2:528628/148‑1 (MQ=255)
                            cAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGTCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCaa                                                                                                                     <  1:1387003/148‑1 (MQ=255)
                            cAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGTCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCaa                                                                                                                     <  1:3035363/148‑1 (MQ=255)
                             aGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAAt                                                                                                                    >  1:1922277/1‑148 (MQ=255)
                               ttGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCa                                                                                                                  <  2:2245910/148‑1 (MQ=255)
                                 ggCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACAAAAATCAGCATCGCGACCAGCGTTTCAATCAGc                                                                                                                >  1:1087015/1‑148 (MQ=255)
                                   cATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGt                                                                                                              <  2:403016/148‑1 (MQ=255)
                                    aTCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGta                                                                                                             <  2:2523662/148‑1 (MQ=255)
                                    aTCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGta                                                                                                             <  2:2523275/148‑1 (MQ=255)
                                    aTCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGta                                                                                                             >  2:2615148/1‑148 (MQ=255)
                                          aCGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCa                                                                                                       <  2:3193796/148‑1 (MQ=255)
                                            gTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACg                                                                                                     <  1:1503253/148‑1 (MQ=255)
                                            gTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACg                                                                                                     >  1:3119207/1‑148 (MQ=255)
                                            gTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACg                                                                                                     <  2:2667501/148‑1 (MQ=255)
                                            gTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACg                                                                                                     <  2:1309070/148‑1 (MQ=255)
                                             tAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGt                                                                                                    >  2:944889/1‑148 (MQ=255)
                                             tAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGt                                                                                                    <  2:2637332/148‑1 (MQ=255)
                                             tAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGt                                                                                                    >  1:3019185/1‑148 (MQ=255)
                                             tAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGt                                                                                                    <  1:264011/148‑1 (MQ=255)
                                               aaTAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTtg                                                                                                  <  1:2321688/148‑1 (MQ=255)
                                                  agagCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTgtgt                                                                                               <  1:2170037/148‑1 (MQ=255)
                                                  agagCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTgtgt                                                                                               <  1:2858845/148‑1 (MQ=255)
                                                  agagCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTgtgt                                                                                               <  1:2859301/148‑1 (MQ=255)
                                                   gagCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGtt                                                                                              >  1:2135888/1‑148 (MQ=255)
                                                   gagCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGtt                                                                                              <  1:1652410/148‑1 (MQ=255)
                                                    agCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGttt                                                                                             >  2:2125319/1‑148 (MQ=255)
                                                    agCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGttt                                                                                             <  2:1408709/148‑1 (MQ=255)
                                                     gCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGtttt                                                                                            <  2:2771015/148‑1 (MQ=255)
                                                     gCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGtttt                                                                                            >  1:3312540/1‑148 (MQ=255)
                                                      cAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTc                                                                                           <  2:2612345/148‑1 (MQ=255)
                                                      cAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTc                                                                                           <  1:2021297/148‑1 (MQ=255)
                                                      cAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTc                                                                                           <  2:1244408/148‑1 (MQ=255)
                                                       aaaTAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCa                                                                                          >  2:1337336/1‑148 (MQ=255)
                                                         aTAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATg                                                                                        <  1:1168966/148‑1 (MQ=255)
                                                             tCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACTAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCgg                                                                                    >  1:2204316/1‑148 (MQ=255)
                                                               ccGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCa                                                                                  <  1:107598/148‑1 (MQ=255)
                                                               ccGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTAGTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCa                                                                                  <  1:107829/148‑1 (MQ=255)
                                                                    ccTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATg                                                                             <  2:2964329/148‑1 (MQ=255)
                                                                     ctggctggCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATgg                                                                            <  1:2988400/148‑1 (MQ=255)
                                                                     ctggctggCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATgg                                                                            <  2:1207589/148‑1 (MQ=255)
                                                                      tggctggcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGa                                                                           <  2:2704847/148‑1 (MQ=255)
                                                                       ggctggcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAg                                                                          <  1:1960232/148‑1 (MQ=255)
                                                                       ggctggcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAg                                                                          <  2:1920775/148‑1 (MQ=255)
                                                                       ggctggcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAg                                                                          <  1:3043110/148‑1 (MQ=255)
                                                                       ggctggcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAg                                                                          <  1:52321/148‑1 (MQ=255)
                                                                          tggcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCga                                                                       <  1:232971/148‑1 (MQ=255)
                                                                            gcGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCgagg                                                                     <  2:1916708/148‑1 (MQ=255)
                                                                             cGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCgagga                                                                    <  2:781645/148‑1 (MQ=255)
                                                                                tttGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGaa                                                                 <  2:3030529/148‑1 (MQ=255)
                                                                                tttACCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGaa                                                                 >  2:582582/1‑148 (MQ=255)
                                                                                   gCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGaaaaa                                                              <  1:2942397/148‑1 (MQ=255)
                                                                                      aTAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAg                                                           <  1:1802875/148‑1 (MQ=255)
                                                                                      aTAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAg                                                           >  1:38614/1‑148 (MQ=255)
                                                                                            gTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAgg                                                     <  2:1918926/148‑1 (MQ=255)
                                                                                                 gACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAg                                                <  2:3247691/148‑1 (MQ=255)
                                                                                                         ccAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTAtt                                        >  1:859885/1‑148 (MQ=255)
                                                                                                         ccAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTAtt                                        <  1:654726/148‑1 (MQ=255)
                                                                                                          cAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCTTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTAttt                                       <  1:1431052/148‑1 (MQ=255)
                                                                                                          cAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTAttt                                       >  2:2005058/1‑148 (MQ=255)
                                                                                                              aTTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATttcttc                                   <  2:790316/148‑1 (MQ=255)
                                                                                                                tGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCgg                                 >  1:845774/1‑148 (MQ=255)
                                                                                                                tGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCgg                                 <  2:2946564/148‑1 (MQ=255)
                                                                                                                tGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCgg                                 >  1:3196749/1‑148 (MQ=255)
                                                                                                                 gccagccaTAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGgc                                <  1:2660075/148‑1 (MQ=255)
                                                                                                                      ccaTAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATc                           <  1:2273584/148‑1 (MQ=255)
                                                                                                                      ccaTAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATc                           >  2:513987/1‑148 (MQ=255)
                                                                                                                           gagTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCg                      >  1:3127410/1‑148 (MQ=255)
                                                                                                                            agTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCgg                     >  2:2403408/1‑148 (MQ=255)
                                                                                                                               ccACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACt                  >  2:446912/1‑148 (MQ=255)
                                                                                                                                cACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTa                 <  1:582582/148‑1 (MQ=255)
                                                                                                                                 ccTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTAt                <  1:1044070/147‑1 (MQ=255)
                                                                                                                                  cTTTCGCTTAGCATGTCCAAAATCAGCATCTCTACCAGCTTTTCAATCAGCTTATAACCACTTTGTGTTTTCATGCCGGAAGTATGGTGCGAGTAGAAAAATAGACGAGGGCCATTTTCTATTTCTTCGGCGCATCTTCCGGACTAtt               <  2:1514628/148‑1 (MQ=255)
                                                                                                                                  cTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTAtt               >  2:2491570/1‑148 (MQ=255)
                                                                                                                                  cTTGCGCTTAGCATGACAAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTAtt               >  1:1688877/1‑148 (MQ=255)
                                                                                                                                   ttGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTAttt              <  1:2753402/148‑1 (MQ=255)
                                                                                                                                      cgcTTAGCATGACAAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACg           <  2:1688877/148‑1 (MQ=255)
                                                                                                                                            gCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTg     >  1:1834911/1‑148 (MQ=255)
                                                                                                                                             cATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGc    <  1:2165971/148‑1 (MQ=255)
                                                                                                                                              aTGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCa   <  2:715361/148‑1 (MQ=255)
                                                                                                                                               tGACAAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAg  <  2:1087015/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AACGCTGATACTGTGGTCGCGGTTATGCCAGTTGGCATCTTCACGTAAATAGAGCAAATAGTCCCGCGCCTGGCTGGCGGTTTGCCATAGCCGTTGCGACTGCTGCCAGTATTGCCAGCCATAGAGTCCACTTGCGCTTAGCATGACCAAAATCAGCATCGCGACCAGCGTTTCAATCAGCGTATAACCACGTTGTGTTTTCATGCCGGCAGTATGGAGCGAGGAGAAAAAAAGACGAGGGCCAGTTTCTATTTCTTCGGCGCATCTTCCGGACTATTTACGCCGTTGCAG  >  CP009273/2957333‑2957623

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: