Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 992,767 C→A 100% D68Y (GAT→TAT)  ssuE ← NAD(P)H‑dependent FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273992,7670CA98.6% 950.9 / 2.2 300D68Y (GAT→TAT) ssuENAD(P)H‑dependent FMN reductase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (157/138);  total (159/141)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 300 total aligned reads displayed.

CCGCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCA  >  CP009273/992620‑992913
                                                                                                                                                   |                                                                                                                                                  
ccGCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTGTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAgg                                                                                                                                                  >  2:4934972/1‑150 (MQ=255)
  gCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGcc                                                                                                                                                <  1:892789/150‑1 (MQ=255)
    cagcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                              >  1:5028978/1‑150 (MQ=255)
     agcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCGGt                                                                                                                                             >  1:1732626/1‑150 (MQ=255)
          aTGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAg                                                                                                                                        >  2:3353673/1‑150 (MQ=255)
           tGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGc                                                                                                                                       <  1:1996065/150‑1 (MQ=255)
                cACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTc                                                                                                                                  >  2:3698881/1‑150 (MQ=255)
                 aCGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCg                                                                                                                                 <  2:512676/150‑1 (MQ=255)
                        cGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAgg                                                                                                                          >  2:468621/1‑150 (MQ=255)
                          cccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTc                                                                                                                        <  1:521762/150‑1 (MQ=255)
                           ccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCt                                                                                                                       >  1:2571821/1‑150 (MQ=255)
                           ccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCt                                                                                                                       >  1:799536/1‑150 (MQ=255)
                            cGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCtt                                                                                                                      >  2:1077940/1‑150 (MQ=255)
                            cGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCtt                                                                                                                      <  2:3056530/150‑1 (MQ=255)
                              tCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGACGGTCTTGa                                                                                                                    >  2:4510554/1‑150 (MQ=255)
                                gCCAGCGGTAGCACCACTTTGCCTTGCACAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGt                                                                                                                  <  2:3741717/150‑1 (MQ=255)
                                  cAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGc                                                                                                                >  2:1278793/1‑150 (MQ=255)
                                  cAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGc                                                                                                                >  2:2467530/1‑150 (MQ=255)
                                  cAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGc                                                                                                                >  2:3136828/1‑150 (MQ=255)
                                   aGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGcc                                                                                                               <  1:1282383/150‑1 (MQ=255)
                                   aGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGcc                                                                                                               <  1:4407507/150‑1 (MQ=255)
                                     cGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCgg                                                                                                             >  2:2779604/1‑150 (MQ=255)
                                      ggTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                            >  2:1147141/1‑150 (MQ=255)
                                       gTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGAc                                                                                                           >  2:3543091/1‑150 (MQ=255)
                                        tAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACt                                                                                                          <  2:1974435/150‑1 (MQ=255)
                                         aGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTa                                                                                                         >  1:2144095/1‑150 (MQ=255)
                                               acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCTGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGaaa                                                                                                   <  1:1391242/150‑1 (MQ=255)
                                               acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGaaa                                                                                                   >  2:2955078/1‑150 (MQ=255)
                                                cTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTCTCGAAAc                                                                                                  >  2:157766/1‑150 (MQ=255)
                                                cTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAAc                                                                                                  >  1:2376258/1‑150 (MQ=255)
                                                   tGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAg                                                                                               >  2:583013/1‑150 (MQ=255)
                                                     ccTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCa                                                                                             <  2:124915/150‑1 (MQ=255)
                                                     ccTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCa                                                                                             <  2:4366664/150‑1 (MQ=255)
                                                       ttGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATa                                                                                           >  2:3586145/1‑150 (MQ=255)
                                                        tGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaa                                                                                          <  1:4957005/150‑1 (MQ=255)
                                                         gCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaaa                                                                                         <  2:870781/150‑1 (MQ=255)
                                                          cAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAg                                                                                        >  1:1830883/1‑150 (MQ=255)
                                                           aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                       >  1:3518057/1‑150 (MQ=255)
                                                            aaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTa                                                                                      <  1:3209006/150‑1 (MQ=255)
                                                             aGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGATTAGGCTGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAg                                                                                     <  1:2932974/150‑1 (MQ=255)
                                                               cgcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGAt                                                                                   >  1:4759298/1‑150 (MQ=255)
                                                               cgcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGAt                                                                                   >  1:4325308/1‑150 (MQ=255)
                                                                      tGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                            >  2:1499175/1‑150 (MQ=255)
                                                                      tGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                            <  2:47967/150‑1 (MQ=255)
                                                                      tGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                            >  1:1540520/1‑150 (MQ=255)
                                                                        gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgggg                                                                          <  2:2104451/150‑1 (MQ=255)
                                                                         cagcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                         >  2:902892/1‑150 (MQ=255)
                                                                         cagcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                         <  1:1459781/150‑1 (MQ=255)
                                                                          agcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCg                                                                        >  1:2859500/1‑150 (MQ=255)
                                                                              ggTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGt                                                                    >  1:4878239/1‑150 (MQ=255)
                                                                                tCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGttt                                                                  >  1:1511939/1‑150 (MQ=255)
                                                                                tCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGttt                                                                  >  1:4075746/1‑150 (MQ=255)
                                                                                 cGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtttt                                                                 <  2:2580345/150‑1 (MQ=255)
                                                                                  gAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTg                                                                <  1:678692/150‑1 (MQ=255)
                                                                                  gAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTg                                                                >  2:4452636/1‑150 (MQ=255)
                                                                                    gCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCa                                                              <  1:3877121/150‑1 (MQ=255)
                                                                                     cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                             >  1:3444249/1‑150 (MQ=255)
                                                                                     cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                             <  1:4552151/150‑1 (MQ=255)
                                                                                     cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                             >  2:3824490/1‑150 (MQ=255)
                                                                                         gTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTc                                                         <  1:4975323/150‑1 (MQ=255)
                                                                                             tCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGt                                                     >  1:4013507/1‑150 (MQ=255)
                                                                                                  gCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaa                                                <  1:2313386/150‑1 (MQ=255)
                                                                                                   cACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaaa                                               >  1:728069/1‑150 (MQ=255)
                                                                                                   cACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaaa                                               >  1:3216642/1‑150 (MQ=255)
                                                                                                    aCCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAc                                              >  2:1816628/1‑150 (MQ=255)
                                                                                                     ccGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAcc                                             <  2:3059738/150‑1 (MQ=255)
                                                                                                      cGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                            >  1:919113/1‑150 (MQ=255)
                                                                                                       ggAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctc                                           >  2:5014511/1‑150 (MQ=255)
                                                                                                       ggAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctc                                           >  2:3954932/1‑150 (MQ=255)
                                                                                                        gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                          <  2:1969015/150‑1 (MQ=255)
                                                                                                        gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                          <  1:1196600/150‑1 (MQ=255)
                                                                                                        gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                          <  2:1442745/150‑1 (MQ=255)
                                                                                                         aaTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTa                                         >  2:3245252/1‑150 (MQ=255)
                                                                                                          aTAGGCGGCTGTATACACAGGCGTGGCGACAAGCAGCCCCTAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTAc                                        <  2:5104186/150‑1 (MQ=255)
                                                                                                             ggcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATc                                     <  2:4057316/150‑1 (MQ=255)
                                                                                                             ggcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATc                                     >  1:2501740/1‑150 (MQ=255)
                                                                                                              gcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATcc                                    >  2:3414237/1‑150 (MQ=255)
                                                                                                               nggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATccc                                   >  1:1039879/2‑149 (MQ=255)
                                                                                                               cggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCa                                   <  2:3566400/150‑1 (MQ=255)
                                                                                                               cggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCa                                   <  2:2376258/150‑1 (MQ=255)
                                                                                                                ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                                  <  1:2897000/150‑1 (MQ=255)
                                                                                                                ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                                  <  2:1795402/150‑1 (MQ=255)
                                                                                                                      ataCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAtt                            >  1:4616758/1‑150 (MQ=255)
                                                                                                                       taCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAttt                           >  1:1765518/1‑150 (MQ=255)
                                                                                                                       taCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAttt                           <  2:3485566/150‑1 (MQ=255)
                                                                                                                            aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                      >  1:3623584/1‑150 (MQ=255)
                                                                                                                            aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                      <  1:3781366/150‑1 (MQ=255)
                                                                                                                            aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                      >  2:771638/1‑150 (MQ=255)
                                                                                                                             ggCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGTTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGtttt                     >  1:4253346/1‑150 (MQ=255)
                                                                                                                              gCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttttt                    <  1:2340362/150‑1 (MQ=255)
                                                                                                                                 tGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTccc                 >  1:1601246/1‑150 (MQ=255)
                                                                                                                                  ggCGACAATCAGCCCATAGGCCTGGTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCcg                >  2:4577967/1‑150 (MQ=255)
                                                                                                                                    cGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCcgcg              <  2:4737076/150‑1 (MQ=255)
                                                                                                                                     gACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCgcgc             >  2:4884694/1‑150 (MQ=255)
                                                                                                                                      aCAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCa            >  2:2241045/1‑150 (MQ=255)
                                                                                                                                      aCAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCa            <  2:5038634/150‑1 (MQ=255)
                                                                                                                                        aaTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata          >  2:162813/1‑150 (MQ=255)
                                                                                                                                            aGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc      <  1:2779604/150‑1 (MQ=255)
                                                                                                                                                cATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCa  >  1:1251011/1‑150 (MQ=255)
                                                                                                                                                cATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCa  <  2:4153021/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCA  >  CP009273/992620‑992913

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