Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 302,428 G→T 100% 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,4280GT97.6% 656.8 / 6.4 212P34Q (CCG→CAG) ecpDpolymerized tip adhesin of ECP fibers
Reads supporting (aligned to +/- strand):  ref base G (0/5);  new base T (98/109);  total (98/114)
Fisher's exact test for biased strand distribution p-value = 6.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 212 total aligned reads displayed.

GCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGG  >  CP009273/302283‑302569
                                                                                                                                                 |                                                                                                                                             
gCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGt                                                                                                                                             >  1:2097596/1‑148 (MQ=255)
  tttGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGc                                                                                                                                           <  1:2452110/148‑1 (MQ=255)
  tttGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGc                                                                                                                                           <  2:3670635/148‑1 (MQ=255)
    tGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCAt                                                                                                                                         <  1:2348892/148‑1 (MQ=255)
           aGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCa                                                                                                                                  >  2:2707270/1‑148 (MQ=255)
             aTGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                >  1:3768614/1‑148 (MQ=255)
             aTGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                >  1:2581860/1‑148 (MQ=255)
             aTGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                >  1:2581511/1‑148 (MQ=255)
             aTGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAgg                                                                                                                                >  1:2465693/1‑148 (MQ=255)
               ggTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgtt                                                                                                                              <  1:3302703/148‑1 (MQ=255)
               ggTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgtt                                                                                                                              >  1:747407/1‑148 (MQ=255)
                gTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttg                                                                                                                             <  2:1640079/148‑1 (MQ=255)
                gTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCAACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttg                                                                                                                             <  1:989108/148‑1 (MQ=255)
                  tCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGgttgtt                                                                                                                           <  2:53517/148‑1 (MQ=255)
                   ccTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGttt                                                                                                                          <  1:1635966/148‑1 (MQ=255)
                         ccAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCa                                                                                                                    <  2:2643797/148‑1 (MQ=255)
                         ccAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCa                                                                                                                    >  2:3036718/1‑148 (MQ=255)
                                 ttttAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGct                                                                                                            >  2:3731653/1‑148 (MQ=255)
                                 ttttAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGct                                                                                                            >  1:2823097/1‑148 (MQ=255)
                                  tttACACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGctc                                                                                                           <  2:466657/148‑1 (MQ=255)
                                    tAAACCGGTCCCGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTg                                                                                                         <  2:1738571/148‑1 (MQ=255)
                                     aaaCCGGGCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGa                                                                                                        <  2:2394470/148‑1 (MQ=255)
                                       aCCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAgc                                                                                                      >  2:2712128/1‑148 (MQ=255)
                                       aCCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAgc                                                                                                      >  2:1067311/1‑148 (MQ=255)
                                       aCCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAgc                                                                                                      >  2:3675322/1‑148 (MQ=255)
                                            tCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGt                                                                                                 >  2:3107631/1‑148 (MQ=255)
                                              cAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTa                                                                                               <  1:492895/148‑1 (MQ=255)
                                               aGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTAc                                                                                              >  1:765960/1‑148 (MQ=255)
                                               aGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTAc                                                                                              >  1:765856/1‑148 (MQ=255)
                                                gCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACt                                                                                             <  1:3391272/148‑1 (MQ=255)
                                                   gTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCCGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGc                                                                                          >  1:1304329/1‑148 (MQ=255)
                                                     tGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCa                                                                                        >  1:691584/1‑148 (MQ=255)
                                                     tGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCa                                                                                        >  1:3408956/1‑148 (MQ=255)
                                                     tGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCa                                                                                        >  2:624138/1‑148 (MQ=255)
                                                       gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                      <  2:714500/148‑1 (MQ=255)
                                                       gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                      <  1:3666327/148‑1 (MQ=255)
                                                       gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                      >  1:3784547/1‑148 (MQ=255)
                                                       gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                      >  1:1538494/1‑148 (MQ=255)
                                                       gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                      <  2:2806707/148‑1 (MQ=255)
                                                       gCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGa                                                                                      <  1:2707270/148‑1 (MQ=255)
                                                         aCCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAACAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGAtt                                                                                    <  2:532651/148‑1 (MQ=255)
                                                         aCCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGAtt                                                                                    <  2:1097272/148‑1 (MQ=255)
                                                              tCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAgcgc                                                                               >  1:1901771/1‑148 (MQ=255)
                                                                  gcgcGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaaa                                                                           <  2:720605/148‑1 (MQ=255)
                                                                  gcgcGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaaa                                                                           <  1:735147/148‑1 (MQ=255)
                                                                  gcgcGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaaa                                                                           <  2:2726055/148‑1 (MQ=255)
                                                                  gcgcGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGaaa                                                                           <  1:2239159/148‑1 (MQ=255)
                                                                     cgGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAAtt                                                                        <  1:2900112/148‑1 (MQ=255)
                                                                         tCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatca                                                                    >  2:1684756/1‑148 (MQ=255)
                                                                         tCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatca                                                                    <  2:3552523/148‑1 (MQ=255)
                                                                         tCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatca                                                                    <  2:1125731/148‑1 (MQ=255)
                                                                          ccAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatcat                                                                   <  2:2333974/148‑1 (MQ=255)
                                                                          ccAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTatcat                                                                   <  1:1603280/148‑1 (MQ=255)
                                                                            aGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCg                                                                 <  2:3194600/148‑1 (MQ=255)
                                                                            aGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCg                                                                 <  2:1293235/148‑1 (MQ=255)
                                                                              cgcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGtt                                                               <  1:2790300/148‑1 (MQ=255)
                                                                              cgcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGtt                                                               >  1:28220/1‑148 (MQ=255)
                                                                               gcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGtta                                                              <  1:2754457/148‑1 (MQ=255)
                                                                               gcCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGtta                                                              <  1:2164029/148‑1 (MQ=255)
                                                                                ccccGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGttat                                                             <  2:533819/148‑1 (MQ=255)
                                                                                      cAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATtagtag                                                       <  2:2603296/148‑1 (MQ=255)
                                                                                       aGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAgg                                                      <  2:1715124/148‑1 (MQ=255)
                                                                                       aGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAgg                                                      <  1:2173817/148‑1 (MQ=255)
                                                                                        ggAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGt                                                     >  2:139181/1‑148 (MQ=255)
                                                                                         gAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGtt                                                    <  1:1405382/148‑1 (MQ=255)
                                                                                         gAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGtt                                                    <  2:2182830/148‑1 (MQ=255)
                                                                                           gTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTaa                                                  <  2:2569082/148‑1 (MQ=255)
                                                                                           gTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTaa                                                  <  1:3163686/148‑1 (MQ=255)
                                                                                            tGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAAc                                                 >  2:1373865/1‑148 (MQ=255)
                                                                                             gACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAAct                                                <  1:3133463/148‑1 (MQ=255)
                                                                                               cGAAAAAGTTGTCCTCTGAGTTGATTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAActct                                              <  1:3557811/148‑1 (MQ=255)
                                                                                                 aaaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCa                                            <  1:2406817/148‑1 (MQ=255)
                                                                                                 aaaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCa                                            <  1:12661/148‑1 (MQ=255)
                                                                                                 aaaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCa                                            >  1:524135/1‑148 (MQ=255)
                                                                                                 aaaaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGGTATTAGTAGGTTAACTCTCa                                            >  1:455863/1‑148 (MQ=255)
                                                                                                   aaaGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAtt                                          <  1:1960288/148‑1 (MQ=255)
                                                                                                     aGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttct                                        <  1:911857/148‑1 (MQ=255)
                                                                                                      gTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctc                                       <  1:2712128/148‑1 (MQ=255)
                                                                                                       ttGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctct                                      <  1:3036718/148‑1 (MQ=255)
                                                                                                       ttGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctct                                      >  1:553498/1‑148 (MQ=255)
                                                                                                        tGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctctt                                     >  1:2415565/1‑148 (MQ=255)
                                                                                                           cgtcTGAGTTGTTTTCGACCAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                  >  1:1774592/1‑148 (MQ=255)
                                                                                                           cgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                  >  1:425491/1‑148 (MQ=255)
                                                                                                           cgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                  <  2:2581816/148‑1 (MQ=255)
                                                                                                           cgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                  <  2:2581511/148‑1 (MQ=255)
                                                                                                           cgtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCAttctcttctc                                  <  1:2916148/148‑1 (MQ=255)
                                                                                                            gtcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTcc                                 <  2:142610/148‑1 (MQ=255)
                                                                                                             tcTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCACATTAGCGCGAAAATTATCAACGTTATTAGTAGGTTAACTCTCATTCTCTTATCCt                                >  1:895448/1‑148 (MQ=255)
                                                                                                              cTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTg                               >  1:2477545/1‑148 (MQ=255)
                                                                                                               tGAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGt                              <  1:767145/148‑1 (MQ=255)
                                                                                                                gAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTc                             >  1:1065203/1‑148 (MQ=255)
                                                                                                                gAGTTGTTTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTc                             >  1:2142503/1‑148 (MQ=255)
                                                                                                                   ttgttTTCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGc                          >  1:2455528/1‑148 (MQ=255)
                                                                                                                      ttttCGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTc                       >  1:2287967/1‑148 (MQ=255)
                                                                                                                          cGACAAACACAAACTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCAccg                   >  2:1597034/1‑148 (MQ=255)
                                                                                                                                  acaaacTCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCAc           >  1:1703636/1‑148 (MQ=255)
                                                                                                                                        tCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATgag     >  2:3427383/1‑148 (MQ=255)
                                                                                                                                        tCGCGTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGGTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATgag     >  2:1324221/1‑148 (MQ=255)
                                                                                                                                           cgTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAgg  <  2:509655/148‑1 (MQ=255)
                                                                                                                                           cgTGCCTGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAgg  <  1:3288355/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GCTTTGCTGATAGATGGTTCCTGAACCAGTGTATTTTAAACCGGTCCAGCGGTTGGCACCGGTCAGGCGCGGATCCAGCGCCCCGCCAGGAGTGACGAAAAAGTTGTCGTCTGAGTTGTTTTCGACAAACACAAACTCGCGTGCCGGTGCATCCGCCCAGGTTGTTTTGCTCACTGCCGCTCTGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGG  >  CP009273/302283‑302569

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: ‑