Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 747,973 C→A 86.0% P93P (CCG→CCT)  ybgD ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273747,9730CA86.0% 728.7 / 131.7 307P93P (CCG→CCT) ybgDputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (9/34);  new base A (115/149);  total (124/183)
Fisher's exact test for biased strand distribution p-value = 6.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 307 total aligned reads displayed.

TGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAT  >  CP009273/747829‑748118
                                                                                                                                                |                                                                                                                                                 
   aCCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTc                                                                                                                                            <  2:1821392/149‑1 (MQ=255)
   aCCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTc                                                                                                                                            <  2:2705381/149‑1 (MQ=255)
     cTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCg                                                                                                                                          >  1:3925561/1‑149 (MQ=255)
      tAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                         <  1:2527265/149‑1 (MQ=255)
         tACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTAc                                                                                                                                      >  2:4673581/1‑149 (MQ=255)
          aCGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTAcc                                                                                                                                     <  1:1348369/149‑1 (MQ=255)
           cGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCg                                                                                                                                    >  1:848142/1‑149 (MQ=255)
              tGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCATCAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTg                                                                                                                                 <  1:1659078/149‑1 (MQ=255)
              tGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTg                                                                                                                                 >  2:4173077/1‑149 (MQ=255)
               gTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGt                                                                                                                                >  1:1166106/1‑149 (MQ=255)
                  aCCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGGTGCTTCGCCTGCACTGGTTTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAg                                                                                                                             <  2:2500724/149‑1 (MQ=255)
                      tCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                         <  2:3122567/149‑1 (MQ=255)
                       cATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgca                                                                                                                        >  1:304198/1‑149 (MQ=255)
                        aTTTTTGTCCATCAGTCGTACACCGACCCCAGTAGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcag                                                                                                                       >  1:3557021/1‑149 (MQ=255)
                            ttGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcaggcag                                                                                                                   >  2:482997/1‑149 (MQ=255)
                              gTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGAt                                                                                                                 >  1:6213/1‑149 (MQ=255)
                              gTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGAt                                                                                                                 <  2:367915/149‑1 (MQ=255)
                                 cATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATcac                                                                                                              >  1:1017036/1‑149 (MQ=255)
                                   tCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                            <  1:2226789/149‑1 (MQ=255)
                                      gTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTg                                                                                                         >  2:3746545/1‑149 (MQ=255)
                                      gTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTg                                                                                                         >  1:2318344/1‑149 (MQ=255)
                                      gTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTg                                                                                                         >  1:3944120/1‑149 (MQ=255)
                                           acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCaa                                                                                                    >  2:198899/1‑149 (MQ=255)
                                           acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCaa                                                                                                    <  1:3156538/149‑1 (MQ=255)
                                           acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCaa                                                                                                    >  1:961598/1‑149 (MQ=255)
                                           acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCaa                                                                                                    <  1:1211862/149‑1 (MQ=255)
                                            cacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAg                                                                                                   <  2:2334186/149‑1 (MQ=255)
                                             acCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGc                                                                                                  >  1:4549220/1‑149 (MQ=255)
                                                   cccAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGt                                                                                            <  2:190047/149‑1 (MQ=255)
                                                         tGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGGTACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACgg                                                                                      >  2:1513042/1‑149 (MQ=255)
                                                           cTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGcc                                                                                    <  1:1818950/149‑1 (MQ=255)
                                                            ttCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                                                                   >  2:2426364/1‑149 (MQ=255)
                                                                  tGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                                                                             >  2:4518253/1‑149 (MQ=255)
                                                                  tGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                                                                             <  1:2478442/149‑1 (MQ=255)
                                                                   gCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAtt                                                                            <  1:141774/149‑1 (MQ=255)
                                                                     aCTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGc                                                                          <  2:4235139/149‑1 (MQ=255)
                                                                      cTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCt                                                                         >  2:694309/1‑149 (MQ=255)
                                                                       tGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTa                                                                        <  1:4056155/149‑1 (MQ=255)
                                                                          tgtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatga                                                                     <  1:3101139/149‑1 (MQ=255)
                                                                           gtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgat                                                                    <  1:2360425/149‑1 (MQ=255)
                                                                           gtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgat                                                                    >  1:4249649/1‑149 (MQ=255)
                                                                            ttGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgatg                                                                   >  2:4459209/1‑149 (MQ=255)
                                                                                tCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc                                                               <  2:1647190/149‑1 (MQ=255)
                                                                                tCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc                                                               <  2:1118867/149‑1 (MQ=255)
                                                                                     aaaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGtt                                                          <  2:4636750/149‑1 (MQ=255)
                                                                                     aaaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGtt                                                          <  2:1268336/149‑1 (MQ=255)
                                                                                      aaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt                                                         <  1:3443609/149‑1 (MQ=255)
                                                                                       aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTta                                                        >  2:3472628/1‑149 (MQ=255)
                                                                                        ggCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTtat                                                       <  2:2956390/149‑1 (MQ=255)
                                                                                         gCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTtata                                                      <  1:1976382/149‑1 (MQ=255)
                                                                                           gTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATg                                                    <  2:4312180/149‑1 (MQ=255)
                                                                                              gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 <  2:322742/149‑1 (MQ=255)
                                                                                              gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 >  1:1200576/1‑149 (MQ=255)
                                                                                              gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 <  2:2715185/149‑1 (MQ=255)
                                                                                                  cAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTc                                             <  2:2224863/149‑1 (MQ=255)
                                                                                                    gTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGt                                           <  2:3615927/149‑1 (MQ=255)
                                                                                                       gTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGAc                                        <  1:740574/149‑1 (MQ=255)
                                                                                                       gTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGAc                                        >  1:2134115/1‑149 (MQ=255)
                                                                                                       gTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTGTATATGGGTTGTCGTGAc                                        >  2:818333/1‑149 (MQ=255)
                                                                                                         cTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACtt                                      <  2:3769240/149‑1 (MQ=255)
                                                                                                          ttAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTg                                     >  1:1098156/1‑149 (MQ=255)
                                                                                                          ttAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTg                                     <  1:2088324/149‑1 (MQ=255)
                                                                                                          ttAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTg                                     >  1:2790858/1‑149 (MQ=255)
                                                                                                           tAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGt                                    <  2:50906/149‑1 (MQ=255)
                                                                                                                gTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAg                               >  2:2271201/1‑149 (MQ=255)
                                                                                                                 tGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCATAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGa                              <  2:3211757/149‑1 (MQ=255)
                                                                                                                 tGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGa                              >  2:1277942/1‑149 (MQ=255)
                                                                                                                  gCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGAt                             <  2:1555491/149‑1 (MQ=255)
                                                                                                                  gCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGAt                             <  1:4397229/149‑1 (MQ=255)
                                                                                                                  gCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCCAGAt                             >  1:4678607/1‑149 (MQ=255)
                                                                                                                   cTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATc                            <  2:1324829/149‑1 (MQ=255)
                                                                                                                   cTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATc                            <  2:3809678/149‑1 (MQ=255)
                                                                                                                   cTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATc                            <  2:3428671/149‑1 (MQ=255)
                                                                                                                     aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                          <  1:694309/149‑1 (MQ=255)
                                                                                                                     aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                          <  1:3679641/149‑1 (MQ=255)
                                                                                                                      tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTAt                         <  2:1209039/149‑1 (MQ=255)
                                                                                                                      tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTAt                         <  2:4137494/149‑1 (MQ=255)
                                                                                                                      tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTAt                         >  2:2603149/1‑149 (MQ=255)
                                                                                                                         aaGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                      <  1:3086777/149‑1 (MQ=255)
                                                                                                                         aaGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                      >  1:4715050/1‑149 (MQ=255)
                                                                                                                         aaGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                      >  1:1287964/1‑149 (MQ=255)
                                                                                                                         aaGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                      >  1:2540974/1‑149 (MQ=255)
                                                                                                                          aGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttt                     >  1:3743754/1‑149 (MQ=255)
                                                                                                                           ggTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtttt                    <  1:4021182/149‑1 (MQ=255)
                                                                                                                           ggTTACGCCAACTTTGGAGACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtttt                    <  2:4397776/149‑1 (MQ=255)
                                                                                                                            tttACGCAAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGTCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttttt                   <  1:749212/148‑1 (MQ=255)
                                                                                                                            gTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttttt                   <  2:4178417/149‑1 (MQ=255)
                                                                                                                             ttACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                  <  1:2940360/149‑1 (MQ=255)
                                                                                                                               aCGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATc                <  2:2096966/149‑1 (MQ=255)
                                                                                                                                 gCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa              >  2:2528840/1‑149 (MQ=255)
                                                                                                                                  ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat             <  2:2642018/149‑1 (MQ=255)
                                                                                                                                   cAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAata            >  1:1578019/1‑149 (MQ=255)
                                                                                                                                    aaCTGTGGATACAGGCAGTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat           <  2:3342161/149‑1 (MQ=255)
                                                                                                                                     aCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtc          >  1:3467354/1‑149 (MQ=255)
                                                                                                                                      cTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtct         <  2:1056299/149‑1 (MQ=255)
                                                                                                                                       tttGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtctt        >  1:3761629/1‑149 (MQ=255)
                                                                                                                                       tttGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtctt        >  2:2942546/1‑149 (MQ=255)
                                                                                                                                             tACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAt  >  2:3889125/1‑149 (MQ=255)
                                                                                                                                             tACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAt  >  1:3593833/1‑149 (MQ=255)
                                                                                                                                             tACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAt  >  1:1560867/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAT  >  CP009273/747829‑748118

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