Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 607,931 G→A 100% 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,9310GA98.0% 719.2 / 2.4 204S7F (TCC→TTC) fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (95/105);  total (98/106)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 213 total aligned reads displayed.

TCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGT  >  CP009273/607784‑608077
                                                                                                                                                   |                                                                                                                                                  
 cTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaat                                                                                                                                                 <  1:664438/150‑1 (MQ=255)
   gCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatga                                                                                                                                               <  1:1144633/150‑1 (MQ=255)
     ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                             <  1:573450/150‑1 (MQ=255)
     ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                             >  1:529387/1‑150 (MQ=255)
       tGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtct                                                                                                                                           <  1:451643/150‑1 (MQ=255)
          tGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtcttct                                                                                                                                        <  2:3275806/150‑1 (MQ=255)
           ggTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATcttctt                                                                                                                                       <  1:727494/150‑1 (MQ=255)
           ggTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATcttctt                                                                                                                                       <  1:1793142/150‑1 (MQ=255)
           ggTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATcttctt                                                                                                                                       >  1:1944722/1‑150 (MQ=255)
            gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTg                                                                                                                                      <  1:1051519/150‑1 (MQ=255)
              cGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGtt                                                                                                                                    <  2:4373609/150‑1 (MQ=255)
               gAAACGCCAGGCGCCTGTTAGTTCTGCTCGGCGGCTGTAACGACAATAGTATCGTCATGTGAAACAGGATTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTc                                                                                                                                   <  1:1291350/150‑1 (MQ=255)
               gAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTc                                                                                                                                   <  1:4458189/150‑1 (MQ=255)
                aaaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCa                                                                                                                                  >  1:2733305/1‑150 (MQ=255)
                aaaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCa                                                                                                                                  <  2:532157/150‑1 (MQ=255)
                 aaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAt                                                                                                                                 >  1:3455100/1‑150 (MQ=255)
                 aaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAt                                                                                                                                 <  1:1558611/150‑1 (MQ=255)
                   cGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTg                                                                                                                               <  1:3264915/150‑1 (MQ=255)
                      cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                            <  2:1806049/150‑1 (MQ=255)
                      cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                            <  1:153763/150‑1 (MQ=255)
                      cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                            <  2:4411976/150‑1 (MQ=255)
                       aGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGtttt                                                                                                                           <  2:1568362/150‑1 (MQ=255)
                           gcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTc                                                                                                                       >  1:2259840/1‑150 (MQ=255)
                              tGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTg                                                                                                                    <  2:62635/150‑1 (MQ=255)
                                       tGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGc                                                                                                           >  2:4263261/1‑150 (MQ=255)
                                              cggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAt                                                                                                    >  1:2206760/1‑150 (MQ=255)
                                                 cggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGc                                                                                                 <  1:2629885/150‑1 (MQ=255)
                                                  ggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCa                                                                                                >  1:1096391/1‑150 (MQ=255)
                                                  ggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCa                                                                                                >  2:2459414/1‑150 (MQ=255)
                                                    tAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAAc                                                                                              >  1:3768080/1‑150 (MQ=255)
                                                     aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTGTTATTCCTGCATTTTTGCCACGAATTGCAACt                                                                                             >  2:2294878/1‑150 (MQ=255)
                                                       cGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATATTTGCCACGAATTGCAACTGt                                                                                           <  2:2982029/150‑1 (MQ=255)
                                                         aCAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCg                                                                                         >  1:4230021/1‑150 (MQ=255)
                                                         aCAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCg                                                                                         >  2:1417598/1‑150 (MQ=255)
                                                         aCAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCACCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCg                                                                                         >  1:1902119/1‑150 (MQ=255)
                                                           aaTAGTATCGTCATGTGAAACAGGATTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGACAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTTCCACGAATTGCAACTGTCggg                                                                                       >  1:4146643/1‑150 (MQ=255)
                                                             tAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCa                                                                                     >  2:2726476/1‑150 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTTCGCTACCCCATAAATCCCCAGATTGACCATCAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    <  1:4100724/150‑1 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    <  1:4349486/150‑1 (MQ=255)
                                                                tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                  >  2:830405/1‑150 (MQ=255)
                                                                   cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcg                                                                               <  1:180838/150‑1 (MQ=255)
                                                                    gTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcgt                                                                              >  2:325582/1‑150 (MQ=255)
                                                                      cATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtca                                                                            <  1:3067218/150‑1 (MQ=255)
                                                                      cATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtca                                                                            >  2:108367/1‑150 (MQ=255)
                                                                          tgAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCaa                                                                        <  2:3953430/150‑1 (MQ=255)
                                                                            aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                      <  2:1691150/150‑1 (MQ=255)
                                                                               cAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacga                                                                   >  1:2185744/1‑150 (MQ=255)
                                                                                aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                  >  1:2948762/1‑150 (MQ=255)
                                                                                aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                  <  1:2835787/150‑1 (MQ=255)
                                                                                aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                  <  2:2841790/150‑1 (MQ=255)
                                                                                 ggAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgacg                                                                 >  1:1686612/1‑150 (MQ=255)
                                                                                   aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                               >  1:712117/1‑150 (MQ=255)
                                                                                        cggtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATAGCAACTGACGGGCATGGTCGTCATCAACACGACGCATCCCg                                                          >  1:2016897/1‑150 (MQ=255)
                                                                                          gtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCt                                                        <  1:435164/150‑1 (MQ=255)
                                                                                           tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTa                                                       <  2:1275483/150‑1 (MQ=255)
                                                                                           tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTa                                                       <  2:4054752/150‑1 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     >  2:579669/1‑150 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     >  1:715116/1‑150 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     <  1:2357055/150‑1 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     <  1:2082407/150‑1 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     >  2:4142863/1‑150 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     >  2:1475100/1‑150 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                     >  1:4173195/1‑150 (MQ=255)
                                                                                              gCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcg                                                    <  2:1587476/150‑1 (MQ=255)
                                                                                              gCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTACTTGTTCATTGTTTTATTACTGCATATTTGCCACGAATTGCAAATGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcct                                                    >  2:4195345/1‑149 (MQ=255)
                                                                                                    gtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaaa                                              >  2:1606698/1‑150 (MQ=255)
                                                                                                    gtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaaa                                              >  1:2141715/1‑150 (MQ=255)
                                                                                                    gtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaaa                                              >  1:1204630/1‑150 (MQ=255)
                                                                                                      gCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAAcc                                            >  2:3262124/1‑150 (MQ=255)
                                                                                                         tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                         <  2:4346477/150‑1 (MQ=255)
                                                                                                         tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                         <  2:118574/150‑1 (MQ=255)
                                                                                                         tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                         <  1:596886/150‑1 (MQ=255)
                                                                                                          gcgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGGTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTg                                        <  2:3728278/150‑1 (MQ=255)
                                                                                                             cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     >  2:2098188/1‑150 (MQ=255)
                                                                                                             cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     >  1:2507713/1‑150 (MQ=255)
                                                                                                                ccccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                                  <  2:3564431/150‑1 (MQ=255)
                                                                                                                 cccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                                 <  1:40053/150‑1 (MQ=255)
                                                                                                                   cATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaaa                               >  1:1135167/1‑150 (MQ=255)
                                                                                                                    aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                              <  1:83061/150‑1 (MQ=255)
                                                                                                                    aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                              <  2:1097301/150‑1 (MQ=255)
                                                                                                                     tAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGa                             <  1:3473230/150‑1 (MQ=255)
                                                                                                                       aaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGaca                           <  1:1094107/150‑1 (MQ=255)
                                                                                                                       aaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGaca                           <  2:3798113/150‑1 (MQ=255)
                                                                                                                         tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                         <  1:2990228/150‑1 (MQ=255)
                                                                                                                         tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                         <  2:643283/150‑1 (MQ=255)
                                                                                                                         tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                         >  2:4544613/1‑150 (MQ=255)
                                                                                                                         tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                         >  2:2655701/1‑150 (MQ=255)
                                                                                                                            ccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcag                      >  2:528461/1‑150 (MQ=255)
                                                                                                                             cAGATTGACCACCAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt                     <  2:1359703/150‑1 (MQ=255)
                                                                                                                              aGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg                    >  1:4517336/1‑150 (MQ=255)
                                                                                                                                aTTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtgca                  >  2:4293743/1‑150 (MQ=255)
                                                                                                                                  tGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtgcagt                <  1:994330/150‑1 (MQ=255)
                                                                                                                                       aacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt           >  1:1802421/1‑150 (MQ=255)
                                                                                                                                         caaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaa         >  2:4188857/1‑150 (MQ=255)
                                                                                                                                         caaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaa         >  2:2350662/1‑150 (MQ=255)
                                                                                                                                            ggCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaatta      >  1:229886/1‑150 (MQ=255)
                                                                                                                                              ccAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat    >  2:271207/1‑150 (MQ=255)
                                                                                                                                               cAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg   <  1:464959/150‑1 (MQ=255)
                                                                                                                                                aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt  >  2:1509549/1‑150 (MQ=255)
                                                                                                                                                aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt  >  1:1397508/1‑150 (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: