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,7670CA99.6% 1023.0 / ‑4.3 318D68Y (GAT→TAT) ssuENAD(P)H‑dependent FMN reductase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (161/154);  total (162/156)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 318 total aligned reads displayed.

GCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCA  >  CP009273/992622‑992913
                                                                                                                                                 |                                                                                                                                                  
gCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGcc                                                                                                                                                >  2:4282360/1‑150 (MQ=255)
  cagcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                              >  2:5307381/1‑150 (MQ=255)
  cagcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                              <  1:3213607/150‑1 (MQ=255)
   agcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGt                                                                                                                                             <  2:1002751/150‑1 (MQ=255)
    gcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGtt                                                                                                                                            <  2:2590040/150‑1 (MQ=255)
     cagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTg                                                                                                                                           >  1:1525729/1‑150 (MQ=255)
      agATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACACTCAGCCCATAGGCCTGTTGc                                                                                                                                          >  2:3497099/1‑150 (MQ=255)
      agATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGc                                                                                                                                          <  1:1382252/150‑1 (MQ=255)
      agATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGc                                                                                                                                          >  2:3813022/1‑150 (MQ=255)
       gATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCa                                                                                                                                         <  2:721850/150‑1 (MQ=255)
       gATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCa                                                                                                                                         <  2:4380299/150‑1 (MQ=255)
        aTGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAg                                                                                                                                        >  2:2816399/1‑150 (MQ=255)
               aCGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCg                                                                                                                                 <  1:3778416/150‑1 (MQ=255)
                cGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCgg                                                                                                                                <  2:3749672/150‑1 (MQ=255)
                cGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCgg                                                                                                                                <  2:1227824/150‑1 (MQ=255)
                 ggTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGt                                                                                                                               >  1:2788770/1‑150 (MQ=255)
                 ggTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGt                                                                                                                               <  2:741273/150‑1 (MQ=255)
                  gTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCATAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTg                                                                                                                              <  1:4182636/150‑1 (MQ=255)
                    aCCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGaa                                                                                                                            <  1:3084436/150‑1 (MQ=255)
                       gCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCATCTGTTCGGTGAAGGt                                                                                                                         <  2:3329777/150‑1 (MQ=255)
                         ccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCt                                                                                                                       <  2:3539681/150‑1 (MQ=255)
                          cGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCtt                                                                                                                      <  2:106051/150‑1 (MQ=255)
                          cGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCtt                                                                                                                      <  2:5201962/150‑1 (MQ=255)
                               ccAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTg                                                                                                                 <  1:70805/150‑1 (MQ=255)
                                cAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGc                                                                                                                >  1:3383168/1‑150 (MQ=255)
                                 aGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGcc                                                                                                               <  2:3960674/150‑1 (MQ=255)
                                  gCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCg                                                                                                              <  2:1554198/150‑1 (MQ=255)
                                   cGGTAGCACCACTTTGCCTTGCACAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCgg                                                                                                             <  1:1572756/150‑1 (MQ=255)
                                    ggTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                            <  2:2517006/150‑1 (MQ=255)
                                    ggTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                            >  2:4765855/1‑150 (MQ=255)
                                     gTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGAc                                                                                                           >  2:3604509/1‑150 (MQ=255)
                                        gCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTAt                                                                                                        >  2:1683042/1‑150 (MQ=255)
                                         caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATc                                                                                                       <  1:564344/150‑1 (MQ=255)
                                         caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATc                                                                                                       <  2:2158465/150‑1 (MQ=255)
                                         caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATc                                                                                                       >  2:3215690/1‑150 (MQ=255)
                                         caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATc                                                                                                       >  2:1766383/1‑150 (MQ=255)
                                          accacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTCTCg                                                                                                      >  2:4051274/1‑150 (MQ=255)
                                          accacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCg                                                                                                      <  2:940299/150‑1 (MQ=255)
                                             acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGaaa                                                                                                   >  1:4052792/1‑150 (MQ=255)
                                              cTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAAc                                                                                                  <  2:1306243/150‑1 (MQ=255)
                                                    cTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCAt                                                                                            >  1:3950836/1‑150 (MQ=255)
                                                     ttGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATa                                                                                           <  1:4609499/150‑1 (MQ=255)
                                                      tGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaa                                                                                          >  1:2497031/1‑150 (MQ=255)
                                                        cAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGTCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAg                                                                                        >  2:4564715/1‑150 (MQ=255)
                                                         aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                       >  2:1949267/1‑150 (MQ=255)
                                                         aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAAAAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGATCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                       >  2:2850823/1‑150 (MQ=255)
                                                            gcgcGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCTACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGTAACGAGCATAAAGTAGa                                                                                    <  2:2371311/150‑1 (MQ=255)
                                                            gcgcGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGa                                                                                    <  2:3245191/150‑1 (MQ=255)
                                                             cgcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGAt                                                                                   <  1:4351755/150‑1 (MQ=255)
                                                                    tGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                            >  1:4520656/1‑150 (MQ=255)
                                                                    tGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                            <  2:2193227/150‑1 (MQ=255)
                                                                    tGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                            <  1:91759/150‑1 (MQ=255)
                                                                      gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgggg                                                                          <  1:2394939/150‑1 (MQ=255)
                                                                      gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgggg                                                                          <  1:4003121/150‑1 (MQ=255)
                                                                       cagcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                         <  1:5157126/150‑1 (MQ=255)
                                                                        agcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGACGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCg                                                                        >  1:4488516/1‑150 (MQ=255)
                                                                         gcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGa                                                                       <  2:436202/150‑1 (MQ=255)
                                                                             gTCGAGCAGGGTTTTCAACGCACCGGTATAGGCGGCTTTATACACCGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtt                                                                   <  1:1298240/150‑1 (MQ=255)
                                                                              tCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGttt                                                                  <  1:3412284/150‑1 (MQ=255)
                                                                               cGAGCAGGGTTTTCAACGCACCGGAATGGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtttt                                                                 >  2:1347107/1‑150 (MQ=255)
                                                                               cGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGACGtttt                                                                 >  2:847629/1‑150 (MQ=255)
                                                                                gAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTg                                                                >  1:4921046/1‑150 (MQ=255)
                                                                                 aGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGc                                                               >  1:2718144/1‑150 (MQ=255)
                                                                                  gCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCa                                                              <  1:3550378/150‑1 (MQ=255)
                                                                                   cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACGCAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                             <  2:489010/150‑1 (MQ=255)
                                                                                   cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                             >  2:2610565/1‑150 (MQ=255)
                                                                                    aGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGa                                                            >  2:968472/1‑150 (MQ=255)
                                                                                    aGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGa                                                            >  2:2931527/1‑150 (MQ=255)
                                                                                         tttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCa                                                       >  2:4510624/1‑150 (MQ=255)
                                                                                         tttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCa                                                       >  1:4748823/1‑150 (MQ=255)
                                                                                           tCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGt                                                     <  2:3699275/150‑1 (MQ=255)
                                                                                            cAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTg                                                    <  1:1472604/150‑1 (MQ=255)
                                                                                            cAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTg                                                    >  2:3842183/1‑150 (MQ=255)
                                                                                             aaCGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGa                                                   <  1:123486/150‑1 (MQ=255)
                                                                                              aCGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGAt                                                  <  2:2022079/150‑1 (MQ=255)
                                                                                               cGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATa                                                 <  2:5066849/150‑1 (MQ=255)
                                                                                                      gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                          >  1:2893763/1‑150 (MQ=255)
                                                                                                       aaTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTa                                         <  1:1029416/150‑1 (MQ=255)
                                                                                                         tAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACa                                       >  2:492787/1‑150 (MQ=255)
                                                                                                             cggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCa                                   <  1:2972661/150‑1 (MQ=255)
                                                                                                              ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                                  >  2:608203/1‑150 (MQ=255)
                                                                                                              ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                                  <  1:1841147/150‑1 (MQ=255)
                                                                                                                    ataCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAtt                            >  1:2416329/1‑150 (MQ=255)
                                                                                                                    ataCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAtt                            >  1:2074554/1‑150 (MQ=255)
                                                                                                                      acacAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTa                          >  1:111330/1‑150 (MQ=255)
                                                                                                                      acacAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTa                          >  1:1439590/1‑150 (MQ=255)
                                                                                                                        acaGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGt                        >  1:5190078/1‑150 (MQ=255)
                                                                                                                          aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                      >  1:3911358/1‑150 (MQ=255)
                                                                                                                            gCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttttt                    <  2:4207573/150‑1 (MQ=255)
                                                                                                                            gCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttttt                    <  2:1525729/150‑1 (MQ=255)
                                                                                                                              gTGGCGACAATCAGCCCATAGGCCTTTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTcc                  <  1:5300128/150‑1 (MQ=255)
                                                                                                                              gTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTcc                  <  2:3180493/150‑1 (MQ=255)
                                                                                                                              gTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTcc                  >  1:4098073/1‑150 (MQ=255)
                                                                                                                                ggCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCcg                <  1:2136784/150‑1 (MQ=255)
                                                                                                                                   gACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCgcgc             >  1:2341379/1‑150 (MQ=255)
                                                                                                                                   gACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCgcgc             <  1:4282360/150‑1 (MQ=255)
                                                                                                                                      aaTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata          <  2:3373496/150‑1 (MQ=255)
                                                                                                                                      aaTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata          >  1:2298326/1‑150 (MQ=255)
                                                                                                                                       aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCatat         >  2:1954859/1‑150 (MQ=255)
                                                                                                                                        tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATAtt        >  1:275651/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCA  >  CP009273/992622‑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: