Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 607,931 G→A 84.4% S7F (TCC→TTC)  fepA ← iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273607,9310GA84.4% 641.9 / 101.8 239S7F (TCC→TTC) fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (15/22);  new base A (96/105);  total (111/128)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 256 total aligned reads displayed.

TCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGT  >  CP009273/607784‑608077
                                                                                                                                                   |                                                                                                                                                  
 cTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaat                                                                                                                                                 <  2:4218263/150‑1 (MQ=255)
   gCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatga                                                                                                                                               <  2:2113123/150‑1 (MQ=255)
    cGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaa                                                                                                                                              >  1:4301572/1‑150 (MQ=255)
    cGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaa                                                                                                                                              <  2:2927427/150‑1 (MQ=255)
     ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                             >  1:449511/1‑150 (MQ=255)
     ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                             >  1:2714340/1‑150 (MQ=255)
      gTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtc                                                                                                                                            >  1:909070/1‑150 (MQ=255)
         aTGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAAtcttc                                                                                                                                         <  2:1724902/150‑1 (MQ=255)
           ggTCGCAACGCCAGGCGCCTGTAAGGTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCCACACGGCCAGGAAATGAATcttctt                                                                                                                                       <  2:3395156/150‑1 (MQ=255)
              cGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGtt                                                                                                                                    >  1:3788594/1‑150 (MQ=255)
               gAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTc                                                                                                                                   >  1:4292071/1‑150 (MQ=255)
               gAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTc                                                                                                                                   >  1:601250/1‑150 (MQ=255)
                 aaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAt                                                                                                                                 >  1:3480906/1‑150 (MQ=255)
                  aCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAtt                                                                                                                                <  2:1466203/150‑1 (MQ=255)
                  aCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAtt                                                                                                                                <  1:468727/150‑1 (MQ=255)
                  aCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAtt                                                                                                                                <  2:3853591/150‑1 (MQ=255)
                   cGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTg                                                                                                                               >  2:1517876/1‑150 (MQ=255)
                    gCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGt                                                                                                                              <  2:947384/150‑1 (MQ=255)
                    gCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGt                                                                                                                              <  2:4176697/150‑1 (MQ=255)
                    gCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGt                                                                                                                              <  1:556247/150‑1 (MQ=255)
                      cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                            <  1:21282/150‑1 (MQ=255)
                       aGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGtttt                                                                                                                           <  1:1722534/150‑1 (MQ=255)
                       aGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGtttt                                                                                                                           <  2:3102146/150‑1 (MQ=255)
                       aGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGtttt                                                                                                                           >  2:1348164/1‑150 (MQ=255)
                       aGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGtttt                                                                                                                           <  2:843013/150‑1 (MQ=255)
                           gcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTc                                                                                                                       >  2:2019533/1‑150 (MQ=255)
                           gcCAGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTc                                                                                                                       <  2:4503847/150‑1 (MQ=255)
                             cTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCt                                                                                                                     <  1:190371/150‑1 (MQ=255)
                             cTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCt                                                                                                                     >  1:4213717/1‑150 (MQ=255)
                             cTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCt                                                                                                                     >  1:1188750/1‑150 (MQ=255)
                              tGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTg                                                                                                                    <  1:541813/150‑1 (MQ=255)
                                tAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCa                                                                                                                  <  1:3015978/150‑1 (MQ=255)
                                  aGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAtt                                                                                                                <  1:2317568/150‑1 (MQ=255)
                                  aGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAtt                                                                                                                >  1:1772194/1‑150 (MQ=255)
                                     tCTGCTCGGCGGCGGTAACGACAATAGTATCGTCCTGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATACCTCCCCCGATTGACCACCAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttttt                                                                                                             <  1:504608/150‑1 (MQ=255)
                                     tCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttttt                                                                                                             >  2:26924/1‑150 (MQ=255)
                                          tCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATACATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCAc                                                                                                        <  2:392447/150‑1 (MQ=255)
                                          tCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCAc                                                                                                        <  1:1496927/150‑1 (MQ=255)
                                          tCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCAc                                                                                                        <  2:4348717/150‑1 (MQ=255)
                                              cggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAt                                                                                                    <  2:1267366/150‑1 (MQ=255)
                                                gcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTg                                                                                                  <  1:2190392/150‑1 (MQ=255)
                                                gcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTg                                                                                                  <  1:679503/150‑1 (MQ=255)
                                                gcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTACCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAACCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTg                                                                                                  <  1:3317563/150‑1 (MQ=255)
                                                    tAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAAc                                                                                              >  2:1138977/1‑150 (MQ=255)
                                                     aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACt                                                                                             <  1:4468630/150‑1 (MQ=255)
                                                        gACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTc                                                                                          <  1:2354249/150‑1 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    <  2:1878531/150‑1 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    <  2:3346229/150‑1 (MQ=255)
                                                                tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                  >  2:3483499/1‑150 (MQ=255)
                                                                      cATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtca                                                                            >  2:3429386/1‑150 (MQ=255)
                                                                           gAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAac                                                                       <  2:465289/150‑1 (MQ=255)
                                                                           gAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAac                                                                       >  2:1908749/1‑150 (MQ=255)
                                                                            aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                      >  2:752349/1‑150 (MQ=255)
                                                                            aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                      >  2:3545927/1‑150 (MQ=255)
                                                                                aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCACCAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                  <  1:319421/150‑1 (MQ=255)
                                                                                   aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCCTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                               <  1:1818720/150‑1 (MQ=255)
                                                                                   aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                               >  1:1994477/1‑150 (MQ=255)
                                                                                   aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                               >  1:4126708/1‑150 (MQ=255)
                                                                                   aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                               >  1:1851227/1‑150 (MQ=255)
                                                                                    gTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCAt                                                              >  2:2685963/1‑150 (MQ=255)
                                                                                        cggtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCg                                                          >  2:158377/1‑150 (MQ=255)
                                                                                          gtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCt                                                        <  2:4041995/150‑1 (MQ=255)
                                                                                          gtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCGTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCt                                                        <  1:3196987/150‑1 (MQ=255)
                                                                                           tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTa                                                       <  1:3709958/150‑1 (MQ=255)
                                                                                            cggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAc                                                      >  2:372770/1‑150 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAACTCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAGTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     <  1:2110314/150‑1 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     >  1:4075860/1‑150 (MQ=255)
                                                                                                  ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaa                                                <  2:1993368/150‑1 (MQ=255)
                                                                                                   tgtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaa                                               >  2:4122397/1‑150 (MQ=255)
                                                                                                       ccTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCt                                           >  1:2845544/1‑150 (MQ=255)
                                                                                                        cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                          <  1:866709/150‑1 (MQ=255)
                                                                                                            gcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                                      <  1:4188903/150‑1 (MQ=255)
                                                                                                            gcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                                      <  1:668088/150‑1 (MQ=255)
                                                                                                             cgaCCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     <  1:2941140/148‑1 (MQ=255)
                                                                                                                ccccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                                  <  1:2801826/150‑1 (MQ=255)
                                                                                                                 cccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                                 <  2:2086202/150‑1 (MQ=255)
                                                                                                                   cATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaaa                               >  1:1070460/1‑150 (MQ=255)
                                                                                                                    aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                              >  1:2356464/1‑150 (MQ=255)
                                                                                                                      aaaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGac                            >  1:1325975/1‑150 (MQ=255)
                                                                                                                       aaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCCTGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGaca                           >  2:4474560/1‑150 (MQ=255)
                                                                                                                        aTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGacac                          <  1:932018/150‑1 (MQ=255)
                                                                                                                        aTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGacac                          >  1:37648/1‑150 (MQ=255)
                                                                                                                         tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                         <  1:4365492/150‑1 (MQ=255)
                                                                                                                         tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                         >  2:4211556/1‑150 (MQ=255)
                                                                                                                          ccccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgc                        >  2:2148143/1‑150 (MQ=255)
                                                                                                                          ccccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgc                        <  2:1061808/150‑1 (MQ=255)
                                                                                                                           cccAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgca                       >  1:3946549/1‑150 (MQ=255)
                                                                                                                            ccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcag                      >  1:1496038/1‑150 (MQ=255)
                                                                                                                             cAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt                     >  1:1611065/1‑150 (MQ=255)
                                                                                                                               gATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtgc                   >  1:1422461/1‑150 (MQ=255)
                                                                                                                                   gACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGtt               >  2:2428266/1‑150 (MQ=255)
                                                                                                                                       aacaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt           >  1:2587734/1‑150 (MQ=255)
                                                                                                                                       aacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt           >  2:4445399/1‑150 (MQ=255)
                                                                                                                                         caaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaa         >  1:2823637/1‑150 (MQ=255)
                                                                                                                                          aaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaat        <  1:3698152/150‑1 (MQ=255)
                                                                                                                                              ccAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat    <  1:220933/150‑1 (MQ=255)
                                                                                                                                               cAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg   <  2:2762599/150‑1 (MQ=255)
                                                                                                                                               cAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg   <  2:3533911/150‑1 (MQ=255)
                                                                                                                                               cAGGAAATGAATCTTCTTGGTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg   <  2:4168744/150‑1 (MQ=255)
                                                                                                                                                aGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt  <  2:3219207/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGT  >  CP009273/607784‑608077

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: