Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 302,428 G→T 29.6% P34Q (CCG→CAG)  ecpD ← polymerized tip adhesin of ECP fibers

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273302,4280GT29.6% 645.7 / 269.1 324P34Q (CCG→CAG) ecpDpolymerized tip adhesin of ECP fibers
Reads supporting (aligned to +/- strand):  ref base G (126/102);  new base T (40/56);  total (166/158)
Fisher's exact test for biased strand distribution p-value = 2.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 334 total aligned reads displayed.

CGAGGCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGC  >  CP009273/302279‑302576
                                                                                                                                                     |                                                                                                                                                    
cGAGGCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCg                                                                                                                                                      >  1:1069212/1‑150 (MQ=255)
   ggCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTg                                                                                                                                                   >  1:1777559/1‑150 (MQ=255)
   ggCTTTGCTGAGAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCtgtg                                                                                                                                                   <  2:4526441/150‑1 (MQ=255)
    gCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGc                                                                                                                                                  <  1:2205029/150‑1 (MQ=255)
      tttGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCAt                                                                                                                                                >  1:2271577/1‑150 (MQ=255)
      tttGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCAt                                                                                                                                                >  2:1621232/1‑150 (MQ=255)
       ttGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATc                                                                                                                                               >  2:2918775/1‑150 (MQ=255)
        tGCTGTTAGAGGGTTCCTGAACTTGTGTATTTTACACCGGTCGAGGGGTTGGCAGCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGTCGAATATGTTGTCGTCTGAGTTGTTTTCGACATACACAAACTCGCGTGCCGGTGCTTcc                                                                                                                                              <  1:42491/150‑1 (MQ=255)
        tGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGGCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATcc                                                                                                                                              <  2:4454704/150‑1 (MQ=255)
             aTAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCa                                                                                                                                         <  1:801966/150‑1 (MQ=255)
             aTAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCa                                                                                                                                         >  1:4625568/1‑150 (MQ=255)
               aGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                       >  1:4787803/1‑150 (MQ=255)
               aGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAgg                                                                                                                                       >  1:1872220/1‑150 (MQ=255)
                 aTGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGgtt                                                                                                                                     <  2:1219467/150‑1 (MQ=255)
                 aTGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGgtt                                                                                                                                     >  2:1155968/1‑150 (MQ=255)
                 aTGGTTCCTGAACCAGGGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgtt                                                                                                                                     >  2:3482295/1‑150 (MQ=255)
                   ggTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttgt                                                                                                                                   >  1:870273/1‑150 (MQ=255)
                   ggTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttgt                                                                                                                                   >  2:2498819/1‑150 (MQ=255)
                         tGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCCGGAGTGACGAAACAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCt                                                                                                                             <  1:4353722/150‑1 (MQ=255)
                           aaCCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCa                                                                                                                           >  1:4515035/1‑150 (MQ=255)
                             ccAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTAGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACt                                                                                                                         >  1:3378346/1‑150 (MQ=255)
                               aGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGc                                                                                                                       >  1:1100291/1‑150 (MQ=255)
                                  gtATTTTAAACCGGTCCAGCGGTTGGCCCCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGc                                                                                                                    >  2:2207544/1‑150 (MQ=255)
                                   tATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGct                                                                                                                   >  1:514587/1‑150 (MQ=255)
                                   tATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGct                                                                                                                   >  1:1739087/1‑150 (MQ=255)
                                     ttttAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGctct                                                                                                                 <  2:3918038/150‑1 (MQ=255)
                                       ttAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGa                                                                                                               <  2:816608/150‑1 (MQ=255)
                                        tAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAg                                                                                                              >  2:2767945/1‑150 (MQ=255)
                                        tAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAg                                                                                                              >  2:2475144/1‑150 (MQ=255)
                                            ccGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCa                                                                                                          >  2:5819/1‑150 (MQ=255)
                                              ggTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGt                                                                                                        >  1:2649380/1‑150 (MQ=255)
                                               gTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGtt                                                                                                       <  2:3233860/150‑1 (MQ=255)
                                               gTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGtt                                                                                                       <  1:2630588/150‑1 (MQ=255)
                                                tCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTa                                                                                                      <  2:3156507/150‑1 (MQ=255)
                                                     cGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGc                                                                                                 >  2:3136535/1‑150 (MQ=255)
                                                       gTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCa                                                                                               >  2:1721860/1‑150 (MQ=255)
                                                        ttGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAg                                                                                              <  2:3671032/150‑1 (MQ=255)
                                                           gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGAtt                                                                                           <  2:2804178/150‑1 (MQ=255)
                                                            cACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTCCTGGCCAGAGTa                                                                                          >  2:4150054/1‑150 (MQ=255)
                                                               cGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAgcg                                                                                       <  1:4592513/150‑1 (MQ=255)
                                                                  tCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGa                                                                                    <  1:218750/150‑1 (MQ=255)
                                                                  tCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGa                                                                                    <  1:2499432/150‑1 (MQ=255)
                                                                   cAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaa                                                                                   <  2:1074497/150‑1 (MQ=255)
                                                                   cAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaa                                                                                   >  1:4030485/1‑150 (MQ=255)
                                                                    aGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaaa                                                                                  >  1:359665/1‑150 (MQ=255)
                                                                       cgcgGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTGTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAAtt                                                                               >  2:428656/1‑150 (MQ=255)
                                                                          ggATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatc                                                                            >  1:925016/1‑150 (MQ=255)
                                                                          ggATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatc                                                                            <  2:1347454/150‑1 (MQ=255)
                                                                          ggATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatc                                                                            >  1:3443270/1‑150 (MQ=255)
                                                                                  cgcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGttat                                                                    <  2:2306778/150‑1 (MQ=255)
                                                                                   gcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGttatt                                                                   >  2:2642686/1‑150 (MQ=255)
                                                                                    ccccGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGttatta                                                                  <  1:4300302/150‑1 (MQ=255)
                                                                                    ccccGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGttatta                                                                  >  1:1299399/1‑150 (MQ=255)
                                                                                        gCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATtagtag                                                              >  1:273657/1‑150 (MQ=255)
                                                                                        gCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATtagtag                                                              <  1:2159865/150‑1 (MQ=255)
                                                                                         ccAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAgg                                                             <  1:695623/150‑1 (MQ=255)
                                                                                          cAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGt                                                            >  2:4065956/1‑150 (MQ=255)
                                                                                           aGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGtt                                                           >  2:1881670/1‑150 (MQ=255)
                                                                                            ggAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTa                                                          >  1:4785611/1‑150 (MQ=255)
                                                                                               gttaCGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAAct                                                       >  1:2044200/4‑150 (MQ=255)
                                                                                                tGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActc                                                      >  2:2858761/1‑150 (MQ=255)
                                                                                                tGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActc                                                      <  1:2992582/150‑1 (MQ=255)
                                                                                                 gACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActct                                                     >  1:4609146/1‑150 (MQ=255)
                                                                                                 gACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActct                                                     <  1:3962142/150‑1 (MQ=255)
                                                                                                 gACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActct                                                     >  2:47343/1‑150 (MQ=255)
                                                                                                 gACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActct                                                     <  1:4239112/150‑1 (MQ=255)
                                                                                                  aCGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACtctc                                                    <  1:1202163/150‑1 (MQ=255)
                                                                                                    gAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAt                                                  >  1:1513738/1‑150 (MQ=255)
                                                                                                     aaaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAtt                                                 >  2:375732/1‑150 (MQ=255)
                                                                                                      aaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttc                                                >  1:98983/1‑150 (MQ=255)
                                                                                                       aaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttct                                               <  1:4336626/150‑1 (MQ=255)
                                                                                                          gTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctctt                                            <  1:3179785/150‑1 (MQ=255)
                                                                                                          gTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctctt                                            <  1:651541/150‑1 (MQ=255)
                                                                                                          gTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctctt                                            <  2:4346024/150‑1 (MQ=255)
                                                                                                             gtcgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                         <  1:2476112/150‑1 (MQ=255)
                                                                                                                 tcTGAGTTGTTGTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGt                                     <  2:3583332/150‑1 (MQ=255)
                                                                                                                  cTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTc                                    >  2:2630805/1‑150 (MQ=255)
                                                                                                                  cTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTc                                    >  2:735791/1‑150 (MQ=255)
                                                                                                                   tGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCt                                   <  2:1119466/150‑1 (MQ=255)
                                                                                                                   tGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCt                                   >  2:4235945/1‑150 (MQ=255)
                                                                                                                    gAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCGGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTg                                  <  1:1549019/150‑1 (MQ=255)
                                                                                                                            ttCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAccg                          >  1:1028247/1‑150 (MQ=255)
                                                                                                                               gACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTCTTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAccgccg                       >  2:3096475/1‑150 (MQ=255)
                                                                                                                               gACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAccgccg                       >  2:1635048/1‑150 (MQ=255)
                                                                                                                                  aaacacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGccc                    >  2:4599078/1‑150 (MQ=255)
                                                                                                                                  aaacacaaacTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGccc                    >  1:2925001/1‑150 (MQ=255)
                                                                                                                                    acacaaacTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCAc                  >  2:3985527/1‑150 (MQ=255)
                                                                                                                                     cacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACg                 <  2:2389275/150‑1 (MQ=255)
                                                                                                                                     cacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACg                 <  2:4787803/150‑1 (MQ=255)
                                                                                                                                     cacaaacTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACg                 >  1:3624704/1‑150 (MQ=255)
                                                                                                                                        aaacTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATg              >  2:612092/1‑150 (MQ=255)
                                                                                                                                         aacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATga             <  2:1515467/150‑1 (MQ=255)
                                                                                                                                           cTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATgaga           <  2:3439963/150‑1 (MQ=255)
                                                                                                                                             cgcgTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAgg         >  1:4644819/1‑150 (MQ=255)
                                                                                                                                               cgTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGcc       <  2:4197569/150‑1 (MQ=255)
                                                                                                                                                gTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCg      >  1:2588314/1‑150 (MQ=255)
                                                                                                                                                  gCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCt    <  1:2100423/150‑1 (MQ=255)
                                                                                                                                                  gCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCt    <  1:4411350/150‑1 (MQ=255)
                                                                                                                                                   ccGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTg   >  2:1419300/1‑150 (MQ=255)
                                                                                                                                                    cGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGc  >  2:1191763/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGGCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGC  >  CP009273/302279‑302576

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: