Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 302,428 G→T 89.7% 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,4280GT89.7% 670.6 / 74.3 254P34Q (CCG→CAG) ecpDpolymerized tip adhesin of ECP fibers
Reads supporting (aligned to +/- strand):  ref base G (18/8);  new base T (113/115);  total (131/123)
Fisher's exact test for biased strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 258 total aligned reads displayed.

AGGCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGC  >  CP009273/302281‑302573
                                                                                                                                                   |                                                                                                                                                 
 ggCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGt                                                                                                                                                 <  2:873178/149‑1 (MQ=255)
  gCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCtgtg                                                                                                                                                >  2:4023410/1‑149 (MQ=255)
       gCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATcc                                                                                                                                           >  1:3061028/1‑149 (MQ=255)
         tGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGc                                                                                                                                         >  2:2089820/1‑149 (MQ=255)
           aTAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGccc                                                                                                                                       >  2:2473322/1‑149 (MQ=255)
           aTAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGccc                                                                                                                                       >  2:2921765/1‑149 (MQ=255)
           aTAGATGGTTCCTGAACCAGTGTAATTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGAGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGccc                                                                                                                                       >  2:1179103/1‑149 (MQ=255)
            tAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCa                                                                                                                                      >  1:1819792/1‑149 (MQ=255)
             aGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAg                                                                                                                                     >  2:166734/1‑149 (MQ=255)
             aGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAg                                                                                                                                     >  2:4202875/1‑149 (MQ=255)
             aGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAg                                                                                                                                     >  2:1764939/1‑149 (MQ=255)
              gATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                    <  1:157859/149‑1 (MQ=255)
              gATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                    >  1:1694915/1‑149 (MQ=255)
                tGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgtt                                                                                                                                  >  1:1312039/1‑149 (MQ=255)
                 ggTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttg                                                                                                                                 >  1:3382687/1‑149 (MQ=255)
                  ggtCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttgt                                                                                                                                <  2:1704947/147‑1 (MQ=255)
                  ggtCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttgt                                                                                                                                <  1:1324715/147‑1 (MQ=255)
                   ttCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttgtt                                                                                                                               <  2:3253855/149‑1 (MQ=255)
                   ttCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGgttgtt                                                                                                                               <  1:68196/149‑1 (MQ=255)
                    tCCTGATCCTGTGTATTTTAATCCGGGCCAGGGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGttt                                                                                                                              <  1:54523/149‑1 (MQ=255)
                     ccTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGtttt                                                                                                                             <  1:2121338/149‑1 (MQ=255)
                      cTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTg                                                                                                                            >  2:540773/1‑149 (MQ=255)
                      cTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTg                                                                                                                            >  1:2700382/1‑149 (MQ=255)
                      cTGAACCAGTGTATTTTAAACCGGTCCAGCGGTGGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTg                                                                                                                            <  1:898321/149‑1 (MQ=255)
                       tGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGc                                                                                                                           >  2:2053676/1‑149 (MQ=255)
                       tGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGc                                                                                                                           >  1:4252600/1‑149 (MQ=255)
                       tGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGc                                                                                                                           <  2:3799734/149‑1 (MQ=255)
                          aCCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCa                                                                                                                        >  1:1043292/1‑149 (MQ=255)
                          aCCAGGGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCCGGTTGTTTTGATCa                                                                                                                        >  1:2064813/1‑149 (MQ=255)
                            cAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACt                                                                                                                      <  1:3145971/149‑1 (MQ=255)
                            cAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACt                                                                                                                      >  2:4569993/1‑149 (MQ=255)
                              gtgtATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGc                                                                                                                    >  2:2929296/1‑149 (MQ=255)
                               tgtATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGcc                                                                                                                   >  2:3185433/1‑149 (MQ=255)
                                 tATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGc                                                                                                                 <  2:431022/149‑1 (MQ=255)
                                 tATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGc                                                                                                                 <  1:646394/149‑1 (MQ=255)
                                  aTTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGct                                                                                                                >  1:3804210/1‑149 (MQ=255)
                                  aTTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGct                                                                                                                >  1:2508566/1‑149 (MQ=255)
                                   ttttAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGctc                                                                                                               <  2:2306167/149‑1 (MQ=255)
                                              tCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGtt                                                                                                    <  2:1429558/149‑1 (MQ=255)
                                                 aGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACt                                                                                                 <  1:2623512/149‑1 (MQ=255)
                                                  gCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCACGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCTCGTGCCTGTGCATCCGCCCAGGTTGTGTTGCTCACTGCCGCTCTGAGCGCAGTTACTg                                                                                                >  1:4714438/1‑149 (MQ=255)
                                                   cGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTgg                                                                                               <  2:3556113/149‑1 (MQ=255)
                                                   cGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTgg                                                                                               <  2:476414/149‑1 (MQ=255)
                                                     gTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGcc                                                                                             <  1:2431863/149‑1 (MQ=255)
                                                     gTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATACGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGcc                                                                                             >  2:3435573/1‑149 (MQ=255)
                                                        ggCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                          <  1:4300853/149‑1 (MQ=255)
                                                        ggCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                          >  1:619346/1‑149 (MQ=255)
                                                        ggCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                          <  1:4576769/149‑1 (MQ=255)
                                                           aCCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTa                                                                                       >  1:2299615/1‑149 (MQ=255)
                                                            ccGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAg                                                                                      <  1:3914995/149‑1 (MQ=255)
                                                             cGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAgc                                                                                     <  1:3828648/149‑1 (MQ=255)
                                                             cGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAgc                                                                                     <  1:3633579/149‑1 (MQ=255)
                                                              ggTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAgcg                                                                                    >  2:1386367/1‑149 (MQ=255)
                                                                  aGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaa                                                                                >  2:4312135/1‑149 (MQ=255)
                                                                     cgcgGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAAt                                                                             <  1:3319292/149‑1 (MQ=255)
                                                                     cgcgGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAAt                                                                             <  2:4743663/149‑1 (MQ=255)
                                                                       cgGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTa                                                                           <  2:635526/149‑1 (MQ=255)
                                                                        ggATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTat                                                                          >  1:1050534/1‑149 (MQ=255)
                                                                           tCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatcat                                                                       <  1:4120183/149‑1 (MQ=255)
                                                                             cAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCg                                                                     >  2:2288798/1‑149 (MQ=255)
                                                                               gcgcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGtt                                                                   <  1:3779444/149‑1 (MQ=255)
                                                                                  ccccGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGttatt                                                                <  2:1813665/149‑1 (MQ=255)
                                                                                    ccGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATtag                                                              <  1:4078729/149‑1 (MQ=255)
                                                                                        cAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAgg                                                          >  2:1000114/1‑149 (MQ=255)
                                                                                         aGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGt                                                         >  2:1483598/1‑149 (MQ=255)
                                                                                           gAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTa                                                       <  2:2216571/149‑1 (MQ=255)
                                                                                             gTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAAc                                                     >  1:3938645/1‑149 (MQ=255)
                                                                                             gTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAAc                                                     >  1:4661874/1‑149 (MQ=255)
                                                                                               gACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGAGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActc                                                   >  2:895613/1‑149 (MQ=255)
                                                                                                aCGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActct                                                  >  2:2037542/1‑149 (MQ=255)
                                                                                                 cGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACtctc                                                 <  2:2797527/149‑1 (MQ=255)
                                                                                                  gAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCa                                                <  2:2700382/149‑1 (MQ=255)
                                                                                                    aaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAtt                                              >  1:2412252/1‑149 (MQ=255)
                                                                                                      aaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttct                                            >  1:1557127/1‑149 (MQ=255)
                                                                                                       aGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctc                                           <  1:1113759/149‑1 (MQ=255)
                                                                                                        gTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctct                                          <  2:399475/149‑1 (MQ=255)
                                                                                                           gtcgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttct                                       <  2:2505339/149‑1 (MQ=255)
                                                                                                            tcgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                      <  2:3877347/149‑1 (MQ=255)
                                                                                                            tcgtcTGAGTCGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                      <  1:1976648/149‑1 (MQ=255)
                                                                                                             cgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTcc                                     <  1:4185603/149‑1 (MQ=255)
                                                                                                              gtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCt                                    >  1:2750340/1‑149 (MQ=255)
                                                                                                                   aGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTg                               <  1:421906/149‑1 (MQ=255)
                                                                                                                    gttgttTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGc                              >  1:2874034/1‑149 (MQ=255)
                                                                                                                     ttgttTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCg                             <  2:1312039/149‑1 (MQ=255)
                                                                                                                      tgttTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGt                            <  2:648335/149‑1 (MQ=255)
                                                                                                                      tgttTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGt                            >  1:147246/1‑149 (MQ=255)
                                                                                                                       gttTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTc                           >  1:2559039/1‑149 (MQ=255)
                                                                                                                          ttCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAcc                        >  2:4625169/1‑149 (MQ=255)
                                                                                                                          ttCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGGTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAcc                        >  1:1227091/1‑149 (MQ=255)
                                                                                                                              acaaacacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAccgccg                    >  2:2080041/1‑149 (MQ=255)
                                                                                                                                 aacacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGccc                 >  2:813003/1‑149 (MQ=255)
                                                                                                                                  acacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCa                <  2:3077589/149‑1 (MQ=255)
                                                                                                                                   cacaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCAc               >  1:4497740/1‑149 (MQ=255)
                                                                                                                                     caaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGa             >  1:4538492/1‑149 (MQ=255)
                                                                                                                                        acTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATga          >  1:2346005/1‑149 (MQ=255)
                                                                                                                                        acTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATga          <  1:4023410/149‑1 (MQ=255)
                                                                                                                                        acTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATga          >  1:3612286/1‑149 (MQ=255)
                                                                                                                                          tCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATgaga        >  2:3476893/1‑149 (MQ=255)
                                                                                                                                               tGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCg   >  1:2308010/1‑149 (MQ=255)
                                                                                                                                               tGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCg   >  1:2492465/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGC  >  CP009273/302281‑302573

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: