Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,647,598 G→T 7.6% R93L (CGA→CTA)  slp → outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,647,5980GT 7.6% 1239.4 / 53.1 328R93L (CGA→CTA) slpouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (152/151);  new base T (11/14);  total (163/165)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 328 total aligned reads displayed.

CAGCCGGGGTTATATGTTGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATCCAG  >  CP009273/3647450‑3647746
                                                                                                                                                    |                                                                                                                                                    
 aGCCGGGGTTATATGTTGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGAc                                                                                                                                                    >  2:4428275/1‑150 (MQ=255)
  gCCGGGGTTATATGTTGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACt                                                                                                                                                   <  1:392924/150‑1 (MQ=255)
       ggTTATATGTTGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCg                                                                                                                                              >  2:3697887/1‑150 (MQ=255)
             atGTTGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGAc                                                                                                                                        >  2:2371858/1‑150 (MQ=255)
                 tGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAg                                                                                                                                    <  2:2526932/150‑1 (MQ=255)
                   gTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCg                                                                                                                                  >  1:3755253/1‑150 (MQ=255)
                    tCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCTACTGCTCGCCAGACAAAGCgg                                                                                                                                 <  2:816134/150‑1 (MQ=255)
                       acaaGCGCGCGTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCtt                                                                                                                              <  2:1590941/150‑1 (MQ=255)
                          aGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCt                                                                                                                           <  1:259721/150‑1 (MQ=255)
                            cgcgcTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAACTCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTg                                                                                                                         <  2:2445312/150‑1 (MQ=255)
                            cgcgcTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCTACTGCTCGCCAGACAAAGCGGCTTCCTTg                                                                                                                         <  1:414351/150‑1 (MQ=255)
                             gcgcGTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGa                                                                                                                        <  2:2006197/150‑1 (MQ=255)
                               gcTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCTAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATc                                                                                                                      >  1:1217423/1‑150 (MQ=255)
                                 tttGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCa                                                                                                                    <  2:2784971/150‑1 (MQ=255)
                                       gggAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAAc                                                                                                              <  2:3307853/150‑1 (MQ=255)
                                        ggAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACt                                                                                                             <  2:3027778/150‑1 (MQ=255)
                                        ggAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACt                                                                                                             >  2:513592/1‑150 (MQ=255)
                                         gAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGACCCAGTGAACTa                                                                                                            <  1:2864194/150‑1 (MQ=255)
                                           aGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATc                                                                                                          <  2:1895491/150‑1 (MQ=255)
                                           aGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATc                                                                                                          >  1:3266077/1‑150 (MQ=255)
                                            ggTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTTTCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCg                                                                                                         <  2:49326/150‑1 (MQ=255)
                                            ggTTATCAACGTGATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCg                                                                                                         <  1:162647/150‑1 (MQ=255)
                                              ttATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTa                                                                                                       <  2:1764439/150‑1 (MQ=255)
                                               tATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTaa                                                                                                      <  1:1758227/150‑1 (MQ=255)
                                                aTCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAAt                                                                                                     >  2:3059101/1‑150 (MQ=255)
                                                aTCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAAt                                                                                                     >  2:123636/1‑150 (MQ=255)
                                                aTCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAAt                                                                                                     >  1:538215/1‑150 (MQ=255)
                                                 tCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATc                                                                                                    <  1:4448500/150‑1 (MQ=255)
                                                  cAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCTACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCa                                                                                                   >  1:497126/1‑150 (MQ=255)
                                                  cAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCa                                                                                                   >  2:2935939/1‑150 (MQ=255)
                                                   aaCGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCAc                                                                                                  >  2:731114/1‑150 (MQ=255)
                                                      gTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACttt                                                                                               >  2:3912139/1‑150 (MQ=255)
                                                      gTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACttt                                                                                               <  1:443526/150‑1 (MQ=255)
                                                      gTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTAGTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACttt                                                                                               <  2:3984960/150‑1 (MQ=255)
                                                       ttATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTg                                                                                              >  2:1663623/1‑150 (MQ=255)
                                                       ttATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTg                                                                                              >  2:1223023/1‑150 (MQ=255)
                                                        tATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGt                                                                                             <  1:680991/150‑1 (MQ=255)
                                                         aTCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGtt                                                                                            >  2:4234128/1‑150 (MQ=255)
                                                          tCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTa                                                                                           >  1:3620058/1‑150 (MQ=255)
                                                          tCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTa                                                                                           <  2:826943/150‑1 (MQ=255)
                                                           cAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTAc                                                                                          <  1:3829639/150‑1 (MQ=255)
                                                           cAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTAc                                                                                          >  1:2249338/1‑150 (MQ=255)
                                                               ggCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATc                                                                                      <  1:4099759/150‑1 (MQ=255)
                                                               ggCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATc                                                                                      >  1:299118/1‑150 (MQ=255)
                                                                gCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATcc                                                                                     >  2:4409726/1‑150 (MQ=255)
                                                                 cAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCt                                                                                    <  2:4220624/150‑1 (MQ=255)
                                                                  aaaaCGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCTACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTc                                                                                   <  1:1810425/150‑1 (MQ=255)
                                                                  aaaaCGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTc                                                                                   >  1:3882411/1‑150 (MQ=255)
                                                                  aaaaCGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTc                                                                                   <  1:2661207/150‑1 (MQ=255)
                                                                   aaaCGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGCTCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCg                                                                                  >  2:607900/1‑150 (MQ=255)
                                                                   aaaCGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCg                                                                                  >  1:2118638/1‑150 (MQ=255)
                                                                     aCGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGc                                                                                >  1:4065112/1‑150 (MQ=255)
                                                                               tgttAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAgg                                                                      <  2:4419674/150‑1 (MQ=255)
                                                                                gttAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAggg                                                                     <  1:1979430/150‑1 (MQ=255)
                                                                                 ttAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGt                                                                    <  2:1981064/150‑1 (MQ=255)
                                                                                  tAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTg                                                                   >  1:1319470/1‑150 (MQ=255)
                                                                                  tAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTg                                                                   >  1:325434/1‑150 (MQ=255)
                                                                                   aGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGa                                                                  <  2:246386/150‑1 (MQ=255)
                                                                                   aGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGa                                                                  >  1:3695067/1‑150 (MQ=255)
                                                                                      aaTCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGaaca                                                               <  2:1321282/150‑1 (MQ=255)
                                                                                        tctcTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGaacaac                                                             >  1:4318102/1‑150 (MQ=255)
                                                                                         ctctGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAAcc                                                            >  1:3179858/1‑150 (MQ=255)
                                                                                         ctctGAATTACCGTTGGATAGCTATTCGAAGCCTGATATTGAAGCCAACTATCAGGGCCTACTGCTCGCCAGACAAAGCGGATTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGaacaaac                                                            >  2:751848/1‑148 (MQ=255)
                                                                                         cgctGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAAcc                                                            >  1:1082728/3‑150 (MQ=255)
                                                                                           ctGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTg                                                          >  1:3210758/1‑150 (MQ=255)
                                                                                           ctGTATGACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTg                                                          <  2:235/150‑1 (MQ=255)
                                                                                             gTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGc                                                        >  1:923822/1‑150 (MQ=255)
                                                                                               aTTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCtt                                                      <  1:2170544/150‑1 (MQ=255)
                                                                                                ttACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCttt                                                     <  2:2201892/150‑1 (MQ=255)
                                                                                                 tACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTa                                                    >  2:1420393/1‑150 (MQ=255)
                                                                                                      ttGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAAt                                               <  2:3225566/150‑1 (MQ=255)
                                                                                                            aGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTc                                         >  1:1399554/1‑150 (MQ=255)
                                                                                                             gCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTcc                                        <  2:3430308/150‑1 (MQ=255)
                                                                                                               tATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCg                                      <  1:4428275/150‑1 (MQ=255)
                                                                                                               tATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCg                                      >  1:628051/1‑150 (MQ=255)
                                                                                                                aTGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGt                                     >  2:1555753/1‑150 (MQ=255)
                                                                                                                 tGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGta                                    >  2:3567541/1‑150 (MQ=255)
                                                                                                                   cGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGtata                                  >  2:205334/1‑150 (MQ=255)
                                                                                                                     aaGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAAc                                >  1:4522900/1‑150 (MQ=255)
                                                                                                                     aaGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAAc                                <  2:3193664/150‑1 (MQ=255)
                                                                                                                      aGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACt                               <  1:871497/150‑1 (MQ=255)
                                                                                                                        ccTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTc                             <  2:1058731/150‑1 (MQ=255)
                                                                                                                          tGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCt                           <  1:3513426/150‑1 (MQ=255)
                                                                                                                              atTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGaa                       >  2:3823294/1‑150 (MQ=255)
                                                                                                                                tGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGt                     <  2:3389112/150‑1 (MQ=255)
                                                                                                                                tGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGt                     <  1:2023369/150‑1 (MQ=255)
                                                                                                                                 gAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTg                    >  1:1049276/1‑150 (MQ=255)
                                                                                                                                 gAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTg                    <  1:2806830/150‑1 (MQ=255)
                                                                                                                                    gCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAat                 <  2:4525069/150‑1 (MQ=255)
                                                                                                                                    gCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAat                 <  1:184140/150‑1 (MQ=255)
                                                                                                                                      cAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAatat               <  1:232977/150‑1 (MQ=255)
                                                                                                                                      cAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAatat               <  2:1108022/150‑1 (MQ=255)
                                                                                                                                       aaCTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATg              <  1:2490136/150‑1 (MQ=255)
                                                                                                                                       aaCTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATg              >  1:2287605/1‑150 (MQ=255)
                                                                                                                                          tATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAg           >  2:2896032/1‑150 (MQ=255)
                                                                                                                                                ggCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATc     <  1:2396834/150‑1 (MQ=255)
                                                                                                                                                 gCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATcc    <  2:1646819/150‑1 (MQ=255)
                                                                                                                                                  ccGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATCCa   <  1:15118/150‑1 (MQ=255)
                                                                                                                                                   cGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATCCAg  <  1:1554889/150‑1 (MQ=255)
                                                                                                                                                   cGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATCCAg  >  2:2276977/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCCGGGGTTATATGTTGGTCAACAAGCGCGCTTTGGTGGGAAGGTTATCAACGTTATCAATGGCAAAACGGATACGTTGTTAGAAATCTCTGTATTACCGTTGGATAGCTATGCGAAGCCTGATATTGAAGCCAACTATCAGGGCCGACTGCTCGCCAGACAAAGCGGCTTCCTTGATCCAGTGAACTATCGTAATCACTTTGTTACCATCCTCGGCACCATTCAGGGTGAACAACCTGGCTTTATCAATAAAGTCCCGTATAACTTCCTGGAAGTGAATATGCAGGGCATCCAG  >  CP009273/3647450‑3647746

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