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,9310GA96.3% 776.3 / 13.1 215S7F (TCC→TTC) fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (101/106);  total (105/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 232 total aligned reads displayed.

TCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGT  >  CP009273/607784‑608077
                                                                                                                                                   |                                                                                                                                                  
tCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaa                                                                                                                                                  >  1:3061714/1‑150 (MQ=255)
  tGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatg                                                                                                                                                <  1:538955/150‑1 (MQ=255)
  tGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatg                                                                                                                                                <  2:3895026/150‑1 (MQ=255)
   gCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatga                                                                                                                                               >  1:1567260/1‑150 (MQ=255)
        gATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtctt                                                                                                                                          <  2:759220/150‑1 (MQ=255)
           ggTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATcttctt                                                                                                                                       <  1:1166597/150‑1 (MQ=255)
            gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTg                                                                                                                                      >  1:3541629/1‑150 (MQ=255)
            gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTg                                                                                                                                      <  2:2134880/150‑1 (MQ=255)
               gAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTc                                                                                                                                   >  2:1442427/1‑150 (MQ=255)
                aaaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCa                                                                                                                                  <  2:3188163/150‑1 (MQ=255)
                  aCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAtt                                                                                                                                <  2:3152573/150‑1 (MQ=255)
                   cGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTg                                                                                                                               <  1:3316716/150‑1 (MQ=255)
                     ccAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGtt                                                                                                                             <  2:1455046/150‑1 (MQ=255)
                      cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGACTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                            >  1:1319836/1‑150 (MQ=255)
                         gcgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                         <  1:2103695/150‑1 (MQ=255)
                          cgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAACTGAATCTTCTTGTTCATTGTTTTAtt                                                                                                                        >  2:1217139/1‑150 (MQ=255)
                          cgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAtt                                                                                                                        >  1:3655975/1‑150 (MQ=255)
                          cgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAtt                                                                                                                        >  1:2908941/1‑150 (MQ=255)
                             cTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCt                                                                                                                     >  1:2170217/1‑150 (MQ=255)
                               gTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGc                                                                                                                   <  1:590626/150‑1 (MQ=255)
                               gTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGc                                                                                                                   <  2:2871623/150‑1 (MQ=255)
                                tAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCa                                                                                                                  >  2:1819941/1‑150 (MQ=255)
                                  aGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAtt                                                                                                                <  1:1303891/150‑1 (MQ=255)
                                   gTTCTGCTCGGCGGCGGTCACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttt                                                                                                               <  1:2449583/150‑1 (MQ=255)
                                   gTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttt                                                                                                               <  2:1610445/150‑1 (MQ=255)
                                   gTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttt                                                                                                               <  2:2027912/150‑1 (MQ=255)
                                   gTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttt                                                                                                               <  2:1860034/150‑1 (MQ=255)
                                    ttCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAtttt                                                                                                              <  2:115612/150‑1 (MQ=255)
                                      cTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTg                                                                                                            >  2:2149991/1‑150 (MQ=255)
                                       tGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGc                                                                                                           >  2:1987089/1‑150 (MQ=255)
                                        gCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGcc                                                                                                          <  1:3610721/150‑1 (MQ=255)
                                         cTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAACTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCa                                                                                                         <  1:3167250/150‑1 (MQ=255)
                                         cTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCa                                                                                                         <  1:3169367/150‑1 (MQ=255)
                                         cTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCa                                                                                                         <  1:933550/150‑1 (MQ=255)
                                          tCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCAc                                                                                                        >  2:3770380/1‑150 (MQ=255)
                                          tCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCAc                                                                                                        >  2:3553513/1‑150 (MQ=255)
                                            ggcggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTGTTATTCCTGCATTTTTGCCACGa                                                                                                      >  2:3484501/1‑150 (MQ=255)
                                             gcggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGaa                                                                                                     >  1:3593634/1‑150 (MQ=255)
                                               ggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAtt                                                                                                   >  1:3199014/1‑150 (MQ=255)
                                               ggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAtt                                                                                                   >  1:2303618/1‑150 (MQ=255)
                                                  ggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCa                                                                                                >  2:365215/1‑150 (MQ=255)
                                                  ggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCa                                                                                                >  2:906635/1‑150 (MQ=255)
                                                      aCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTg                                                                                            >  1:4072548/1‑150 (MQ=255)
                                                      aCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTg                                                                                            >  2:412486/1‑150 (MQ=255)
                                                       cGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGt                                                                                           <  2:1024677/150‑1 (MQ=255)
                                                        gACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTc                                                                                          <  1:2545960/150‑1 (MQ=255)
                                                        gACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTc                                                                                          <  2:1011663/150‑1 (MQ=255)
                                                            aTAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGc                                                                                      >  2:3348117/1‑150 (MQ=255)
                                                            aTAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGc                                                                                      >  2:3346621/1‑150 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    >  1:4139702/1‑150 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    >  1:1949159/1‑150 (MQ=255)
                                                              aGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCAt                                                                                    <  1:1173573/150‑1 (MQ=255)
                                                                tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                  >  2:533843/1‑150 (MQ=255)
                                                                 aTCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgt                                                                                 >  2:1407651/1‑150 (MQ=255)
                                                                    gTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcgt                                                                              <  1:3825647/150‑1 (MQ=255)
                                                                    gTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcgt                                                                              <  2:4291861/150‑1 (MQ=255)
                                                                       aTGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtcat                                                                           >  2:177519/1‑150 (MQ=255)
                                                                            aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                      <  2:1509787/150‑1 (MQ=255)
                                                                            aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAGCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                      >  2:3942403/1‑150 (MQ=255)
                                                                                aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                  >  2:1049701/1‑150 (MQ=255)
                                                                                aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                  <  1:2720517/150‑1 (MQ=255)
                                                                                 ggAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgacg                                                                 <  2:199528/150‑1 (MQ=255)
                                                                                      aTCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCGACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATcc                                                            >  1:473963/1‑150 (MQ=255)
                                                                                          gtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCGTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCt                                                        <  1:1259398/150‑1 (MQ=255)
                                                                                             ggCTCTTGTGCCTGCGCTACCCCATAAATCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcg                                                    >  1:1022625/1‑150 (MQ=255)
                                                                                                tctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                  >  2:3258541/1‑150 (MQ=255)
                                                                                                   tgtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaa                                               <  1:2149991/150‑1 (MQ=255)
                                                                                                     tgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAAc                                             >  1:4093147/1‑150 (MQ=255)
                                                                                                         tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                         <  2:1063243/150‑1 (MQ=255)
                                                                                                         tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                         <  2:2419772/150‑1 (MQ=255)
                                                                                                           cgcTACCCCATACATCCCCAGATTGACCAACACGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGa                                       <  1:3698581/150‑1 (MQ=255)
                                                                                                           cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGa                                       <  1:3697848/150‑1 (MQ=255)
                                                                                                             cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     <  2:3284359/150‑1 (MQ=255)
                                                                                                             cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     >  1:983273/1‑150 (MQ=255)
                                                                                                             cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     >  2:1244699/1‑150 (MQ=255)
                                                                                                             cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                     >  1:1019867/1‑150 (MQ=255)
                                                                                                              tACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATcc                                    <  1:2350748/150‑1 (MQ=255)
                                                                                                              tACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATcc                                    <  1:608784/150‑1 (MQ=255)
                                                                                                                  ccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaa                                <  2:701773/150‑1 (MQ=255)
                                                                                                                   cATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaaa                               >  2:2434817/1‑150 (MQ=255)
                                                                                                                    aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                              >  1:519621/1‑150 (MQ=255)
                                                                                                                     tAAATCCCCGGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGa                             >  2:2652326/1‑150 (MQ=255)
                                                                                                                     tAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGa                             <  1:3687218/150‑1 (MQ=255)
                                                                                                                          ccccAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgc                        >  2:4252661/1‑150 (MQ=255)
                                                                                                                           cccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgca                       <  1:1987089/150‑1 (MQ=255)
                                                                                                                              aGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg                    <  2:98007/150‑1 (MQ=255)
                                                                                                                              aGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg                    <  2:2274202/150‑1 (MQ=255)
                                                                                                                              aGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg                    >  2:2364981/1‑150 (MQ=255)
                                                                                                                              aGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg                    <  1:2824745/150‑1 (MQ=255)
                                                                                                                                      caacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGt            <  1:73150/150‑1 (MQ=255)
                                                                                                                                       aacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt           <  2:3010104/150‑1 (MQ=255)
                                                                                                                                       aacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt           <  2:3011161/150‑1 (MQ=255)
                                                                                                                                         caaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaa         >  1:1461651/1‑150 (MQ=255)
                                                                                                                                          aaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaat        >  1:4244001/1‑150 (MQ=255)
                                                                                                                                             gCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaattaa     <  1:2174339/150‑1 (MQ=255)
                                                                                                                                              ccAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat    <  2:633064/150‑1 (MQ=255)
                                                                                                                                              ccAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat    >  1:871814/1‑150 (MQ=255)
                                                                                                                                               cAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg   >  1:3339322/1‑150 (MQ=255)
                                                                                                                                                aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt  >  1:2574603/1‑150 (MQ=255)
                                                                                                                                                aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt  >  1:1616984/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: