Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 301,998 T→A 15.3% T177T (ACA→ACT ecpD ← polymerized tip adhesin of ECP fibers

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273301,9980TA15.3% 540.8 / 80.5 196T177T (ACA→ACTecpDpolymerized tip adhesin of ECP fibers
Reads supporting (aligned to +/- strand):  ref base T (90/76);  new base A (19/11);  total (109/87)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 198 total aligned reads displayed.

GCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTG  >  CP009273/301853‑302145
                                                                                                                                                 |                                                                                                                                                   
gCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGt                                                                                                                                                   >  1:389991/1‑148 (MQ=255)
   aGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTaaa                                                                                                                                                <  1:2299746/148‑1 (MQ=255)
    gTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAg                                                                                                                                               <  1:1468493/148‑1 (MQ=255)
     tGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgc                                                                                                                                              <  2:50513/148‑1 (MQ=255)
       ggTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgctg                                                                                                                                            <  1:1038230/148‑1 (MQ=255)
       ggTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgctg                                                                                                                                            <  1:1197626/148‑1 (MQ=255)
        gTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgctgc                                                                                                                                           <  1:2850864/148‑1 (MQ=255)
        gTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgctgc                                                                                                                                           <  1:642380/148‑1 (MQ=255)
         tATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAgctgct                                                                                                                                          >  1:331145/1‑148 (MQ=255)
          aTTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGctgctg                                                                                                                                         >  2:2470968/1‑148 (MQ=255)
          aTTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGctgctg                                                                                                                                         >  1:263367/1‑148 (MQ=255)
              aTCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGa                                                                                                                                     >  2:2886231/1‑148 (MQ=255)
                  aCCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACggg                                                                                                                                 <  2:2999405/148‑1 (MQ=255)
                   ccGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGc                                                                                                                                >  1:727339/1‑148 (MQ=255)
                    cGTAGATTTGCTGCTTTCGTAGGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCa                                                                                                                               <  2:1243351/148‑1 (MQ=255)
                     gTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCAt                                                                                                                              <  1:1075019/148‑1 (MQ=255)
                      tAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCAtt                                                                                                                             >  2:820878/1‑148 (MQ=255)
                        gATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATttg                                                                                                                           <  2:2924711/148‑1 (MQ=255)
                         aTTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATttgt                                                                                                                          <  2:2306206/148‑1 (MQ=255)
                         aTTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATttgt                                                                                                                          <  2:2086867/148‑1 (MQ=255)
                         aTTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATttgt                                                                                                                          >  2:1619934/1‑148 (MQ=255)
                              ctgctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCat                                                                                                                     <  1:2185769/148‑1 (MQ=255)
                              ctgctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCat                                                                                                                     >  2:1886658/1‑148 (MQ=255)
                               tgctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCata                                                                                                                    >  1:1065824/1‑148 (MQ=255)
                               tgctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCata                                                                                                                    <  2:3211281/148‑1 (MQ=255)
                                gctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCatat                                                                                                                   >  2:3115446/1‑148 (MQ=255)
                                 ctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGTTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACTGGCATTTGTTGCAaata                                                                                                                  >  1:853489/1‑148 (MQ=255)
                                 ctTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCAtata                                                                                                                  <  2:2044887/148‑1 (MQ=255)
                                    tCGTATGGGTGACGGTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTg                                                                                                               <  1:161177/148‑1 (MQ=255)
                                    tCGTAGGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTg                                                                                                               <  2:2884048/148‑1 (MQ=255)
                                       tATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTa                                                                                                            >  2:506284/1‑148 (MQ=255)
                                        aTGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTaa                                                                                                           >  1:3129022/1‑148 (MQ=255)
                                         tGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTaaa                                                                                                          <  1:1146174/148‑1 (MQ=255)
                                          gggTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTaaaa                                                                                                         <  2:3195652/148‑1 (MQ=255)
                                          gggTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTaaaa                                                                                                         <  1:1996509/148‑1 (MQ=255)
                                           ggTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAAc                                                                                                        >  2:2466899/1‑148 (MQ=255)
                                               aCGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGt                                                                                                    >  1:1030954/1‑148 (MQ=255)
                                                  ttGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTg                                                                                                 <  1:3088993/148‑1 (MQ=255)
                                                   tGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGa                                                                                                <  1:1379634/148‑1 (MQ=255)
                                                        aCATACGAGTTGCCGGAGGACTGATGCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTTGCAGGCATTGATGGTCATTGTTAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGacaaca                                                                                           <  2:1920232/148‑1 (MQ=255)
                                                        aCATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGacaaca                                                                                           <  2:2425036/148‑1 (MQ=255)
                                                        aCATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGacaaca                                                                                           <  1:3067679/148‑1 (MQ=255)
                                                        aCATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGacaaca                                                                                           >  1:3198716/1‑148 (MQ=255)
                                                         cATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACAt                                                                                          >  1:1015344/1‑148 (MQ=255)
                                                           tACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGc                                                                                        >  1:244329/1‑148 (MQ=255)
                                                           tACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGc                                                                                        <  1:590378/148‑1 (MQ=255)
                                                            aCCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGcc                                                                                       >  2:526299/1‑148 (MQ=255)
                                                             ccAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCg                                                                                      >  2:1294822/1‑148 (MQ=255)
                                                             ccAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCAGTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCg                                                                                      >  2:2947587/1‑148 (MQ=255)
                                                               aGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTg                                                                                    >  2:326539/1‑148 (MQ=255)
                                                                gTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGc                                                                                   >  2:1187655/1‑148 (MQ=255)
                                                                 ttGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCa                                                                                  >  2:2860552/1‑148 (MQ=255)
                                                                  tGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCAt                                                                                 <  1:1214380/148‑1 (MQ=255)
                                                                  tGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCAt                                                                                 <  1:2662060/148‑1 (MQ=255)
                                                                      ggAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCAc                                                                             <  1:407265/148‑1 (MQ=255)
                                                                        aGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACtt                                                                           >  2:855777/1‑148 (MQ=255)
                                                                        aGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACtt                                                                           >  2:854531/1‑148 (MQ=255)
                                                                         ggCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTACAACGCGTCTGACAACATGCCGTGCATCCACTTc                                                                          >  1:487989/1‑148 (MQ=255)
                                                                         ggCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTc                                                                          >  2:1842587/1‑148 (MQ=255)
                                                                         ggCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTc                                                                          >  1:2150168/1‑148 (MQ=255)
                                                                          gCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTcg                                                                         <  2:3081423/148‑1 (MQ=255)
                                                                           ccTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTcgc                                                                        >  2:1982612/1‑148 (MQ=255)
                                                                           ccTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTcgc                                                                        >  2:2732545/1‑148 (MQ=255)
                                                                             tGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCa                                                                              <  2:1358551/140‑1 (MQ=255)
                                                                                  cTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGcgccgccg                                                                 <  2:1065824/148‑1 (MQ=255)
                                                                                         cGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCa                                                          <  2:2131469/148‑1 (MQ=255)
                                                                                           tGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCAcc                                                        <  2:2620947/148‑1 (MQ=255)
                                                                                               ccgcTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTcg                                                    >  1:913888/1‑148 (MQ=255)
                                                                                               ccgcTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTcg                                                    >  1:1766419/1‑148 (MQ=255)
                                                                                                  ctgctgGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCgcgc                                                 >  1:1412459/1‑148 (MQ=255)
                                                                                                  ctgctgGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCgcgc                                                 >  1:3057211/1‑148 (MQ=255)
                                                                                                        gCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATaaa                                           >  2:2920867/1‑148 (MQ=255)
                                                                                                        gCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATaaa                                           >  2:2921048/1‑148 (MQ=255)
                                                                                                           tCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATc                                        >  1:2578196/1‑148 (MQ=255)
                                                                                                             aTAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGAGGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCa                                      <  2:356174/148‑1 (MQ=255)
                                                                                                             aTAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCa                                      <  2:1814477/148‑1 (MQ=255)
                                                                                                             aTAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACGTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCa                                      <  2:356171/148‑1 (MQ=255)
                                                                                                              tAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCAc                                     <  1:2256610/148‑1 (MQ=255)
                                                                                                              tAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCAc                                     >  1:2942816/1‑148 (MQ=255)
                                                                                                                gTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTg                                   >  2:2416933/1‑148 (MQ=255)
                                                                                                                 ttCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTgg                                  >  1:2780110/1‑148 (MQ=255)
                                                                                                                     caGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTc                              >  1:612356/1‑148 (MQ=255)
                                                                                                                      agagGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCg                             >  2:551099/1‑148 (MQ=255)
                                                                                                                       gagGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCgg                            <  1:2698626/148‑1 (MQ=255)
                                                                                                                        agGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCggt                           >  1:1403039/1‑148 (MQ=255)
                                                                                                                         ggTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCggtg                          <  1:1911880/148‑1 (MQ=255)
                                                                                                                         ggTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCggtg                          <  1:820878/148‑1 (MQ=255)
                                                                                                                         ggTCTGGCAGGCATTGATGGTCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCggtg                          >  2:1083967/1‑148 (MQ=255)
                                                                                                                            cTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGtt                       <  2:1434514/148‑1 (MQ=255)
                                                                                                                              ggcaggcaTTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTg                     <  2:3077492/148‑1 (MQ=255)
                                                                                                                              ggcaggcaTTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTg                     <  2:920444/148‑1 (MQ=255)
                                                                                                                               gcaggcaTTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGc                    >  1:2395422/1‑148 (MQ=255)
                                                                                                                                caggcaTTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCg                   >  1:2026203/1‑148 (MQ=255)
                                                                                                                                     aTTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAgg              >  1:1685377/1‑148 (MQ=255)
                                                                                                                                      ttGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGa             >  1:2924928/1‑148 (MQ=255)
                                                                                                                                          tGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAg         <  1:854535/148‑1 (MQ=255)
                                                                                                                                          tGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAg         <  1:855777/148‑1 (MQ=255)
                                                                                                                                             tCATAGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACt      >  1:15315/1‑148 (MQ=255)
                                                                                                                                                ttGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACtggt   <  2:623167/148‑1 (MQ=255)
                                                                                                                                                 tGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACtggtg  <  2:2498731/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGGCATTGATGGTCATTGTAAAGCTGCTGCCGACGGGCATTTGTTGCATATACTGGTAAAACGCGTCTGACAACATGCCGTGCATCCACTTCGCGCCGCCGCTGGTCACCGTCGCGCCATAAAATCCACTGGCGTCGGTGGTTTGCGGCAGGATAAGACTGGTG  >  CP009273/301853‑302145

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: