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,7670CA97.2% 651.8 / 7.1 213D68Y (GAT→TAT) ssuENAD(P)H‑dependent FMN reductase
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (85/122);  total (88/125)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 213 total aligned reads displayed.

GCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCA  >  CP009273/992622‑992910
                                                                                                                                                 |                                                                                                                                               
gCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAgg                                                                                                                                               >  2:2346000/1‑148 (MQ=255)
 ccAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGc                                                                                                                                              <  2:1283538/148‑1 (MQ=255)
    gcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                           >  2:1997317/1‑148 (MQ=255)
    gcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                           >  2:2542891/1‑148 (MQ=255)
    gcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                           <  2:1231870/148‑1 (MQ=255)
     cagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGt                                                                                                                                          <  1:2004593/148‑1 (MQ=255)
     cagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGt                                                                                                                                          <  2:509334/148‑1 (MQ=255)
             ccACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTg                                                                                                                                  >  1:2654511/1‑148 (MQ=255)
             ccACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTg                                                                                                                                  <  2:2591566/148‑1 (MQ=255)
                    aCCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTg                                                                                                                           <  1:2264646/148‑1 (MQ=255)
                     ccGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGa                                                                                                                          <  2:1222356/148‑1 (MQ=255)
                      cGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGACTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGaa                                                                                                                         >  2:645648/1‑148 (MQ=255)
                       gCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAg                                                                                                                        <  2:1626681/148‑1 (MQ=255)
                         ccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCANCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGt                                                                                                                      <  1:1780287/148‑1 (MQ=255)
                         ccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGt                                                                                                                      <  1:1964909/148‑1 (MQ=255)
                          cGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTc                                                                                                                     <  1:2186132/148‑1 (MQ=255)
                           gTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCAGAGGCCTGTTGCAGCTGTTCGGTGAAGGTCt                                                                                                                    >  1:3425477/1‑148 (MQ=255)
                             cGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTg                                                                                                                  <  1:1691865/148‑1 (MQ=255)
                               ccAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTACAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAg                                                                                                                <  1:3403725/148‑1 (MQ=255)
                                  gCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGc                                                                                                             >  1:598993/1‑148 (MQ=255)
                                      tAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                         <  1:892713/148‑1 (MQ=255)
                                       aGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGAc                                                                                                        <  2:2284691/148‑1 (MQ=255)
                                       aGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGAc                                                                                                        >  1:767786/1‑148 (MQ=255)
                                        gCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACt                                                                                                       >  1:3236568/1‑148 (MQ=255)
                                         caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTa                                                                                                      <  2:767255/148‑1 (MQ=255)
                                         caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTa                                                                                                      >  2:847649/1‑148 (MQ=255)
                                             acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGa                                                                                                  <  2:3397900/148‑1 (MQ=255)
                                                 tGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACg                                                                                              <  2:2235601/148‑1 (MQ=255)
                                                 tGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACg                                                                                              <  1:1225511/148‑1 (MQ=255)
                                                   ccTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAg                                                                                            <  2:1366449/148‑1 (MQ=255)
                                                    cTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGc                                                                                           <  1:2412286/148‑1 (MQ=255)
                                                     ttGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCa                                                                                          <  1:2630037/148‑1 (MQ=255)
                                                        cAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaa                                                                                       <  1:622727/148‑1 (MQ=255)
                                                         aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaaa                                                                                      <  2:2858496/148‑1 (MQ=255)
                                                          aaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAg                                                                                     >  2:188790/1‑148 (MQ=255)
                                                            gcgcGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTa                                                                                   <  2:2799297/148‑1 (MQ=255)
                                                              gcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGa                                                                                 <  2:1517533/148‑1 (MQ=255)
                                                               cgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGAt                                                                                >  1:165735/1‑148 (MQ=255)
                                                               cgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGAt                                                                                >  2:1590958/1‑148 (MQ=255)
                                                                  tCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCtt                                                                             >  2:1553349/1‑148 (MQ=255)
                                                                  tCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCtt                                                                             <  2:3116038/148‑1 (MQ=255)
                                                                   cTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTc                                                                            <  1:1033705/148‑1 (MQ=255)
                                                                     ggCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCg                                                                          <  1:3502320/148‑1 (MQ=255)
                                                                         gcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                      >  2:597470/1‑148 (MQ=255)
                                                                            ggTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGaa                                                                   >  2:783367/1‑148 (MQ=255)
                                                                             gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                  <  1:1852833/148‑1 (MQ=255)
                                                                             gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                  <  1:2321890/148‑1 (MQ=255)
                                                                              tCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGt                                                                 <  2:147342/148‑1 (MQ=255)
                                                                                 aGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtttt                                                              >  2:1541101/1‑148 (MQ=255)
                                                                                  gCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTg                                                             <  2:3242466/148‑1 (MQ=255)
                                                                                   cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGc                                                            >  2:398381/1‑148 (MQ=255)
                                                                                     gggTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                          >  2:1857628/1‑148 (MQ=255)
                                                                                       gTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAt                                                        >  2:3085033/1‑148 (MQ=255)
                                                                                        ttttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGACACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAtt                                                       >  1:3302292/1‑148 (MQ=255)
                                                                                        ttttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAtt                                                       <  1:3151946/148‑1 (MQ=255)
                                                                                             aaCGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGt                                                  <  1:479784/148‑1 (MQ=255)
                                                                                              aCGCACCGGAATAGGCGGCTTTATACACCGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTg                                                 <  2:2044379/148‑1 (MQ=255)
                                                                                                gCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGAt                                               >  1:2700847/1‑148 (MQ=255)
                                                                                                   ccGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaaa                                            >  1:745618/1‑148 (MQ=255)
                                                                                                    cGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAc                                           >  2:1672619/1‑148 (MQ=255)
                                                                                                    cGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGACGTTTTGCAGATTCCAGTGATAAAc                                           >  1:710269/1‑148 (MQ=255)
                                                                                                    cGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAc                                           >  1:2084583/1‑148 (MQ=255)
                                                                                                    cGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTCCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAc                                           >  1:2443458/1‑148 (MQ=255)
                                                                                                     ggAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAcc                                          >  2:1844884/1‑148 (MQ=255)
                                                                                                     ggAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAcc                                          <  2:3445878/148‑1 (MQ=255)
                                                                                                      taaTAGGCGGCTTTATACTCAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTATTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                         <  2:1894367/147‑1 (MQ=255)
                                                                                                      gAATAGGCGGCTTTATACCCAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                         <  2:1655124/148‑1 (MQ=255)
                                                                                                      gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                         <  2:1895519/148‑1 (MQ=255)
                                                                                                      gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                         <  1:438223/148‑1 (MQ=255)
                                                                                                       aaTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctc                                        <  2:2386227/148‑1 (MQ=255)
                                                                                                        aTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                       <  2:2654506/148‑1 (MQ=255)
                                                                                                            gcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACAt                                   >  1:3012825/1‑148 (MQ=255)
                                                                                                             cggcTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATc                                  <  2:3350688/148‑1 (MQ=255)
                                                                                                              ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATcc                                 <  2:1798049/148‑1 (MQ=255)
                                                                                                               gcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCa                                >  2:3421342/1‑148 (MQ=255)
                                                                                                                cTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                               >  1:351473/1‑148 (MQ=255)
                                                                                                                cTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                               >  1:2120898/1‑148 (MQ=255)
                                                                                                                  ttATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGc                             >  1:1648506/1‑148 (MQ=255)
                                                                                                                  ttATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGc                             >  2:3015065/1‑148 (MQ=255)
                                                                                                                   tataCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCTAGTGATAAACCTCTACATCCAGGcc                            <  1:647060/148‑1 (MQ=255)
                                                                                                                   tataCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGcc                            <  1:2533062/148‑1 (MQ=255)
                                                                                                                    ataCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCa                           <  1:509412/148‑1 (MQ=255)
                                                                                                                     taCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAt                          >  2:66766/1‑148 (MQ=255)
                                                                                                                      acacAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAtt                         >  2:420625/1‑148 (MQ=255)
                                                                                                                       cacaGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAttt                        <  1:1837788/148‑1 (MQ=255)
                                                                                                                       cacaGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAttt                        <  2:2426465/148‑1 (MQ=255)
                                                                                                                         caGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAg                      <  1:188850/148‑1 (MQ=255)
                                                                                                                         caGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAg                      <  1:2918221/148‑1 (MQ=255)
                                                                                                                          aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGt                     <  1:2348873/148‑1 (MQ=255)
                                                                                                                                ggCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTcc               <  1:1541101/148‑1 (MQ=255)
                                                                                                                                 gCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTccc              >  1:1461686/1‑148 (MQ=255)
                                                                                                                                  cGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCcg             >  1:1543422/1‑148 (MQ=255)
                                                                                                                                  cGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCcg             >  2:2666870/1‑148 (MQ=255)
                                                                                                                                     cAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCgcgc          >  1:2596408/1‑148 (MQ=255)
                                                                                                                                       aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat        >  2:1324882/1‑148 (MQ=255)
                                                                                                                                        tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata       >  1:1496473/1‑148 (MQ=255)
                                                                                                                                          aGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATAtt     <  1:1553349/148‑1 (MQ=255)
                                                                                                                                            cccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc   <  2:1130145/148‑1 (MQ=255)
                                                                                                                                            cccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc   >  2:886604/1‑148 (MQ=255)
                                                                                                                                             ccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCa  >  2:1742025/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCA  >  CP009273/992622‑992910

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