Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 109,050 C→T 9.2% E13K (GAA→AAA)  zapD ← FtsZ stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273109,0500CT 9.2% 1242.4 / 86.2 364E13K (GAA→AAA) zapDFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base C (172/157);  new base T (24/9);  total (197/167)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 367 total aligned reads displayed.

CTCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAT  >  CP009273/108901‑109198
                                                                                                                                                     |                                                                                                                                                    
     tGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTa                                                                                                                                                 <  1:4525091/150‑1 (MQ=255)
         gACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTgg                                                                                                                                             >  1:212743/1‑150 (MQ=255)
              cgccgcGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGACATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGtt                                                                                                                                        <  1:1067106/150‑1 (MQ=255)
               gccgcGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTc                                                                                                                                       >  2:3794264/1‑150 (MQ=255)
                 cgcgCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaa                                                                                                                                     <  1:4066965/150‑1 (MQ=255)
                 cgcgCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaa                                                                                                                                     >  1:4337379/1‑150 (MQ=255)
                  gcgcTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaaa                                                                                                                                    >  1:422781/1‑150 (MQ=255)
                         aaaaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGAcc                                                                                                                             <  1:316316/150‑1 (MQ=255)
                            aCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTgg                                                                                                                          >  1:963399/1‑150 (MQ=255)
                              aTCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGt                                                                                                                        <  1:2526941/150‑1 (MQ=255)
                              aTCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTc                                                                                                                       <  2:715099/150‑1 (MQ=255)
                                ccAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCt                                                                                                                      <  1:904504/150‑1 (MQ=255)
                                ccAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCt                                                                                                                      <  2:1947606/150‑1 (MQ=255)
                                 cAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTg                                                                                                                     >  2:1429659/1‑150 (MQ=255)
                                  aGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGc                                                                                                                    >  2:4168000/1‑150 (MQ=255)
                                  aGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGc                                                                                                                    >  2:1261808/1‑150 (MQ=255)
                                  aGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGc                                                                                                                    <  2:2601786/150‑1 (MQ=255)
                                    tAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCAt                                                                                                                  >  1:3358005/1‑150 (MQ=255)
                                      aTTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTa                                                                                                                <  1:4755805/150‑1 (MQ=255)
                                      aTTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTa                                                                                                                <  2:957291/150‑1 (MQ=255)
                                       ttCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCACTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTAc                                                                                                               <  2:4139403/150‑1 (MQ=255)
                                       ttCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTAc                                                                                                               >  2:2547625/1‑150 (MQ=255)
                                        tCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACg                                                                                                              <  1:5140501/150‑1 (MQ=255)
                                            tGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGttt                                                                                                          >  2:1646346/1‑150 (MQ=255)
                                             gACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGtttt                                                                                                         <  1:4171809/150‑1 (MQ=255)
                                                     ggAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTga                                                                                                 <  2:2952585/150‑1 (MQ=255)
                                                       aagaaATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTgaga                                                                                               >  1:2937385/1‑150 (MQ=255)
                                                         gaaATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACa                                                                                             <  1:873928/150‑1 (MQ=255)
                                                           aaTGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACaaa                                                                                           <  1:203016/150‑1 (MQ=255)
                                                            aTGCAGCGCGCCAGCGTGGTCAACGATGGGTTAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACaaag                                                                                          <  2:3561358/150‑2 (MQ=255)
                                                             tGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACt                                                                                         <  2:4666692/150‑1 (MQ=255)
                                                               cAGCTCTCCAGCGTGGTCAACTTTGGGTAAATTAACTGTGAGTTGCTTAATAAAATACTCAATGCGCAGACATGTACGCATTTTTTAATTTATTTGATGTTCAAAAACGAACTGGGTCTGCATTACGGTTTTTCCTGTGAGACACACTGc                                                                                       <  2:1832551/150‑1 (MQ=255)
                                                               cAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGc                                                                                       >  1:1005609/1‑150 (MQ=255)
                                                                aGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCg                                                                                      >  1:1723176/1‑150 (MQ=255)
                                                                  cgcgcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGAc                                                                                    >  2:2455109/1‑150 (MQ=255)
                                                                      ccAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCaa                                                                                <  2:153624/150‑1 (MQ=255)
                                                                      ccAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCaa                                                                                <  2:153504/150‑1 (MQ=255)
                                                                        aGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGc                                                                              >  1:2657083/1‑150 (MQ=255)
                                                                        aGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGc                                                                              >  1:1408565/1‑150 (MQ=255)
                                                                         gCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCt                                                                             >  2:1969607/1‑150 (MQ=255)
                                                                            tGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCa                                                                          >  1:610213/1‑150 (MQ=255)
                                                                             ggTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCaa                                                                         >  2:640701/1‑150 (MQ=255)
                                                                             ggTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCaa                                                                         >  1:532782/1‑150 (MQ=255)
                                                                               tCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAAt                                                                       >  1:3450338/1‑150 (MQ=255)
                                                                               tCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAAt                                                                       <  1:3049653/150‑1 (MQ=255)
                                                                               tCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAAt                                                                       >  1:3331138/1‑150 (MQ=255)
                                                                                 aaCGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAg                                                                     <  1:3599680/150‑1 (MQ=255)
                                                                                      tGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTgtgc                                                                >  1:4606608/1‑150 (MQ=255)
                                                                                       gggTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTgtgcg                                                               <  2:3981052/150‑1 (MQ=255)
                                                                                         gTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTgtgcgtg                                                             >  2:3439080/1‑150 (MQ=255)
                                                                                           aaaTTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCa                                                           >  2:4333401/1‑150 (MQ=255)
                                                                                               tAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCg                                                       <  1:2352096/150‑1 (MQ=255)
                                                                                                aaCGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCgg                                                      >  2:1608713/1‑150 (MQ=255)
                                                                                                aaCGGTGAGTTGCTGAATCAAAAAATCAATGCGCAGCAATGTACGCATTTTTTTATTTATTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCgg                                                      >  1:2014148/1‑150 (MQ=255)
                                                                                                 aCGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGACATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCAGGGTCTGCATTAAGGATATTACTGTGAGACAAACTGCTACCCAAGCTGCAAATAGTGTGCGTGCAGGCggg                                                     >  2:2252283/1‑150 (MQ=255)
                                                                                                    gTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCat                                                  <  2:4588519/150‑2 (MQ=255)
                                                                                                     tGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAAc                                                 >  2:1133494/1‑150 (MQ=255)
                                                                                                     tGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAAc                                                 >  2:2692736/1‑150 (MQ=255)
                                                                                                      gAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACa                                                <  1:5097646/150‑1 (MQ=255)
                                                                                                       aGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACAt                                               >  2:3212363/1‑150 (MQ=255)
                                                                                                         ttGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATcc                                             >  2:5020879/1‑150 (MQ=255)
                                                                                                         ttGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATcc                                             <  1:1300715/150‑1 (MQ=255)
                                                                                                            cTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGAt                                          >  2:1780660/1‑150 (MQ=255)
                                                                                                              gAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGc                                        >  2:2001614/1‑150 (MQ=255)
                                                                                                              gAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGc                                        <  2:3358005/150‑1 (MQ=255)
                                                                                                              gAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGc                                        <  2:368292/150‑1 (MQ=255)
                                                                                                              gAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGc                                        <  1:809712/150‑1 (MQ=255)
                                                                                                               aaTCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCg                                       >  1:1566363/1‑150 (MQ=255)
                                                                                                                aTCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGa                                      >  1:2294197/1‑150 (MQ=255)
                                                                                                                  cAAAAACTCAATGCGCAGACATGTACTCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTATGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATa                                    <  1:3526565/150‑1 (MQ=255)
                                                                                                                     aaaCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCa                                 <  1:2960495/150‑1 (MQ=255)
                                                                                                                        cTCAATGCGCAGCCATGTACGCATTTTTTCAGTGAGTGGATGGTCAAAAAGGACCTGGGTCTGCATTACGGTGGTTCCTGTGAGACAAACTGCGACGCAAGCTGCCAATAGTGTGCGTGCAGGCGGGCAACCTCCGATGCGATCGCAGcc                              >  1:3080830/1‑150 (MQ=255)
                                                                                                                         tCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCg                             <  2:817867/150‑1 (MQ=255)
                                                                                                                         tCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCg                             >  1:383200/1‑150 (MQ=255)
                                                                                                                         tCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCg                             >  1:383013/1‑150 (MQ=255)
                                                                                                                           aaTGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGt                           >  1:1987375/1‑150 (MQ=255)
                                                                                                                            aTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTg                          >  1:160314/1‑150 (MQ=255)
                                                                                                                              gcgcAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTgcg                        <  2:4404770/150‑1 (MQ=255)
                                                                                                                               cgcAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGAATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCAGGTgcgc                       >  1:729427/1‑150 (MQ=255)
                                                                                                                               cgcAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTgcgc                       >  2:2086687/1‑150 (MQ=255)
                                                                                                                                gcAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGcc                      >  2:3102256/1‑150 (MQ=255)
                                                                                                                                  aGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGt                    >  1:4480951/1‑150 (MQ=255)
                                                                                                                                     cATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttat                 >  1:4870083/1‑150 (MQ=255)
                                                                                                                                      aTGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttaat                >  2:240552/1‑148 (MQ=255)
                                                                                                                                       tGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttatta               >  1:2729095/1‑150 (MQ=255)
                                                                                                                                          aCGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCa            >  2:2807576/1‑150 (MQ=255)
                                                                                                                                           cGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCaa           <  2:4755439/150‑1 (MQ=255)
                                                                                                                                             cATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATg         <  2:2704952/150‑1 (MQ=255)
                                                                                                                                               ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc       <  1:2824392/150‑1 (MQ=255)
                                                                                                                                               ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc       >  2:42567/1‑150 (MQ=255)
                                                                                                                                               ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc       <  2:4962184/150‑1 (MQ=255)
                                                                                                                                                tttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACg      >  1:646635/1‑150 (MQ=255)
                                                                                                                                                 ttttTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGt     >  2:1311421/1‑150 (MQ=255)
                                                                                                                                                 ttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGt     >  1:18723/1‑150 (MQ=255)
                                                                                                                                                  tttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTc    <  1:3971272/150‑1 (MQ=255)
                                                                                                                                                   ttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCa   >  1:5018458/1‑150 (MQ=255)
                                                                                                                                                    tCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAt  <  1:4739541/150‑1 (MQ=255)
                                                                                                                                                    tCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAt  <  2:2964942/150‑1 (MQ=255)
                                                                                                                                                    tCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAt  >  1:1158770/1‑150 (MQ=255)
                                                                                                                                                    tCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAt  >  2:846432/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCAT  >  CP009273/108901‑109198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: