Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 301,808 C→A 6.8% E241* (GAA→TAA)  ecpD ← polymerized tip adhesin of ECP fibers

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273301,8080CA 6.8% 1046.1 / 38.5 280E241* (GAA→TAA) ecpDpolymerized tip adhesin of ECP fibers
Reads supporting (aligned to +/- strand):  ref base C (118/143);  new base A (8/11);  total (126/154)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 284 total aligned reads displayed.

GGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGG  >  CP009273/301659‑301957
                                                                                                                                                     |                                                                                                                                                     
ggCCGAGGCTAACGACGAGTTGGCGATCGCCGGGACTATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTc                                                                                                                                                       >  2:3343887/1‑150 (MQ=255)
    gAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGccc                                                                                                                                                   <  2:1787998/150‑1 (MQ=255)
          aaCGACGAGTTGGCGATCGCCGGGAATATATTGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTc                                                                                                                                             <  1:5232751/150‑1 (MQ=255)
          aaCGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTc                                                                                                                                             <  2:577790/150‑1 (MQ=255)
          aaCGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTc                                                                                                                                             <  2:479103/150‑1 (MQ=255)
           acgacgAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCg                                                                                                                                            <  1:3950247/150‑1 (MQ=255)
           acgacgAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCg                                                                                                                                            >  2:2218794/1‑150 (MQ=255)
            cgacgaGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCgg                                                                                                                                           <  1:1023834/150‑1 (MQ=255)
             gacgaGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGt                                                                                                                                          >  2:3401583/1‑150 (MQ=255)
              tcttgTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGACCGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTa                                                                                                                                         <  2:3119676/146‑1 (MQ=255)
              acgaGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTa                                                                                                                                         <  2:3119523/150‑1 (MQ=255)
                gaGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACt                                                                                                                                       <  2:1768876/150‑1 (MQ=255)
                gaGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACt                                                                                                                                       <  1:1535060/150‑1 (MQ=255)
                       cGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGc                                                                                                                                >  1:1875627/1‑150 (MQ=255)
                       cGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTAGCCCAGAGTCGGTACTCCGTCGc                                                                                                                                >  2:3938011/1‑150 (MQ=255)
                         aTCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTg                                                                                                                              >  2:3603074/1‑150 (MQ=255)
                         aTCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTg                                                                                                                              <  1:708709/150‑1 (MQ=255)
                              cGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGAt                                                                                                                         <  1:1801919/150‑1 (MQ=255)
                                ggAATATAGGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATaa                                                                                                                       <  1:2719224/150‑1 (MQ=255)
                                 gAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATaaa                                                                                                                      >  1:3626585/1‑150 (MQ=255)
                                    tatatGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATAc                                                                                                                   <  1:5191135/150‑1 (MQ=255)
                                    tatatGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATAc                                                                                                                   <  1:4702943/150‑1 (MQ=255)
                                    tatatGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATAc                                                                                                                   <  2:4317005/150‑1 (MQ=255)
                                     atatGGATTGATGTGTTGCTTAGTCCGTTTGTGTGCAGGGTATAGTTAACCATCTTACAACTTATTCCTGTACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACt                                                                                                                  <  2:2824497/150‑1 (MQ=255)
                                     atatGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGACCGGCAGTCGGCGTTCCCAACGCCCAGAGTCGGTACTCCGACGCAGTTGATAAAAAcg                                                                                                                  >  2:4889261/1‑149 (MQ=255)
                                        tGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTcc                                                                                                               <  1:4845217/150‑1 (MQ=255)
                                         ggATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCg                                                                                                              <  2:40766/150‑1 (MQ=255)
                                         ggATTGAGGTGGTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCg                                                                                                              <  1:4366392/150‑1 (MQ=255)
                                          gATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGc                                                                                                             >  2:122800/1‑150 (MQ=255)
                                           aTTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGcc                                                                                                            >  2:2155325/1‑150 (MQ=255)
                                              gAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGc                                                                                                         >  1:4068434/1‑150 (MQ=255)
                                                ggtggtGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGa                                                                                                       <  1:3744892/150‑1 (MQ=255)
                                                  tgtTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGt                                                                                                     <  1:1516218/150‑1 (MQ=255)
                                                   gtTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTAGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTg                                                                                                    <  2:152709/150‑1 (MQ=255)
                                                     tGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTggg                                                                                                  >  2:4553427/1‑150 (MQ=255)
                                                      gCTGAGTCCGTTTGATTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGt                                                                                                 <  1:3554266/150‑1 (MQ=255)
                                                       cTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGTGTGGGTa                                                                                                <  2:2344965/150‑1 (MQ=255)
                                                        tGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTAt                                                                                               >  1:4483887/1‑150 (MQ=255)
                                                         gAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTAtt                                                                                              <  2:5088355/150‑1 (MQ=255)
                                                         gAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTAtt                                                                                              <  1:4559277/150‑1 (MQ=255)
                                                           gTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCATTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATAATTACGCCAGCGAGTGGGTATTTa                                                                                            >  2:1086117/1‑150 (MQ=255)
                                                           gTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTa                                                                                            <  1:948134/150‑1 (MQ=255)
                                                            tCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTAt                                                                                           <  2:4915798/150‑1 (MQ=255)
                                                              cGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCa                                                                                         <  1:2047964/150‑1 (MQ=255)
                                                               gtttgtttGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCaa                                                                                        >  1:4126371/1‑150 (MQ=255)
                                                                tttgtttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAAc                                                                                       <  1:3483220/150‑1 (MQ=255)
                                                                tttgtttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAAc                                                                                       <  1:2951500/150‑1 (MQ=255)
                                                                attgtttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAAc                                                                                       <  2:2765417/149‑1 (MQ=255)
                                                                 ttgtttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAAcc                                                                                      <  1:2413247/150‑1 (MQ=255)
                                                                  tgtttgCATGGTATAGTTAACCATCTTACAACTTAATCCTTCACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCg                                                                                     <  2:1267653/150‑1 (MQ=255)
                                                                  tgtttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCg                                                                                     <  2:3883216/150‑1 (MQ=255)
                                                                   gttttCAGTGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGt                                                                                    <  1:5169019/150‑1 (MQ=255)
                                                                    tttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTa                                                                                   >  1:2686166/1‑150 (MQ=255)
                                                                    tttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTa                                                                                   >  1:1435837/1‑150 (MQ=255)
                                                                       gCAGGGTATAGTTAACCCTCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGAt                                                                                <  2:2027108/150‑1 (MQ=255)
                                                                             tataGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTtgctg                                                                          >  1:379532/1‑150 (MQ=255)
                                                                                aGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGAAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCtt                                                                       <  2:1982653/150‑1 (MQ=255)
                                                                                 gTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTCTTTATCAACCGTAGATTTGCTGCttt                                                                      >  2:3730099/1‑150 (MQ=255)
                                                                                 gTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCttt                                                                      <  1:1625131/150‑1 (MQ=255)
                                                                                 gTTAACCATCCTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCttt                                                                      <  2:5089326/150‑1 (MQ=255)
                                                                                    aaCCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGt                                                                   >  1:3951029/1‑150 (MQ=255)
                                                                                     aCCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCCGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTa                                                                  >  1:4833616/1‑150 (MQ=255)
                                                                                     aCCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTa                                                                  >  1:2618673/1‑150 (MQ=255)
                                                                                     aCCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTa                                                                  >  1:2614877/1‑150 (MQ=255)
                                                                                      ccATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTAt                                                                 >  1:4108241/1‑150 (MQ=255)
                                                                                      ccATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTAt                                                                 >  1:4410614/1‑150 (MQ=255)
                                                                                      ccATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTAt                                                                 <  2:3715502/150‑1 (MQ=255)
                                                                                            tACAACTTAATCCTGAACGGCTGCCTTTTGTTTGCGTCCGGCTGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATATATACTTACGACAGCGAGTGGGTTTTTATCAACCGTAGATTTGCTGCTTTCGTTAGGGTGa                                                           <  2:1716078/150‑1 (MQ=255)
                                                                                            tACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGa                                                           <  1:495101/150‑1 (MQ=255)
                                                                                             aCAACTTAATCCTGACCGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGAc                                                          <  2:4704852/150‑1 (MQ=255)
                                                                                             aCAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGAc                                                          <  1:5164934/150‑1 (MQ=255)
                                                                                               aaCTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGt                                                        >  2:343894/1‑150 (MQ=255)
                                                                                                aCTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGtt                                                       >  2:909246/1‑150 (MQ=255)
                                                                                                 cTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTg                                                      >  2:4067195/1‑150 (MQ=255)
                                                                                                 cTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACGTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTg                                                      >  2:2579345/1‑150 (MQ=255)
                                                                                                     aTCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGaa                                                  >  2:3056431/1‑150 (MQ=255)
                                                                                                     aTCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTAGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCCGCGAGGGGGGATTTATCAACCGTAGAGTTGCTGCTTTCGTATGGGTGACGTTGCGaa                                                  >  1:2224603/1‑150 (MQ=255)
                                                                                                       ccTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACa                                                <  1:2095037/150‑1 (MQ=255)
                                                                                                       ccTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACa                                                <  1:381063/150‑1 (MQ=255)
                                                                                                        cTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACAt                                               <  1:2295724/150‑1 (MQ=255)
                                                                                                           aaCGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATAcc                                            >  2:2645968/1‑150 (MQ=255)
                                                                                                           aaCGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATAcc                                            >  1:3783529/1‑150 (MQ=255)
                                                                                                             cGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAg                                          <  1:1083962/150‑1 (MQ=255)
                                                                                                             cGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAg                                          >  1:5144589/1‑150 (MQ=255)
                                                                                                              ggCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGt                                         <  1:1472113/150‑1 (MQ=255)
                                                                                                                  gCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGcc                                     <  1:1267112/150‑1 (MQ=255)
                                                                                                                   ccGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCg                                    >  1:4249170/1‑150 (MQ=255)
                                                                                                                   ccGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCg                                    >  1:4603554/1‑150 (MQ=255)
                                                                                                                     gATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGa                                  >  1:266426/1‑150 (MQ=255)
                                                                                                                        ggTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGGTGCGAACATACCAGTTGCCGGAGGc                               >  1:1770163/1‑150 (MQ=255)
                                                                                                                           ttGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTg                            <  1:718507/150‑1 (MQ=255)
                                                                                                                               gTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATcc                        >  1:3298989/1‑150 (MQ=255)
                                                                                                                                tCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCt                       <  1:3156090/150‑1 (MQ=255)
                                                                                                                                  cGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTa                     <  1:1165636/150‑1 (MQ=255)
                                                                                                                                      aGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGc                 >  2:5133984/1‑150 (MQ=255)
                                                                                                                                         cGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTg              >  2:745340/1‑150 (MQ=255)
                                                                                                                                           gCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTgcg            >  1:3108178/1‑150 (MQ=255)
                                                                                                                                                 ccTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCgctg      <  2:5170170/150‑1 (MQ=255)
                                                                                                                                                  cTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCgctgc     <  1:668272/150‑1 (MQ=255)
                                                                                                                                                    tCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGctgctg   <  1:3343657/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGG  >  CP009273/301659‑301957

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: