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.1% 568.0 / 6.6 188D68Y (GAT→TAT) ssuENAD(P)H‑dependent FMN reductase
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base A (106/76);  total (108/80)
Fisher's exact test for biased strand distribution p-value = 4.04e-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 188 total aligned reads displayed.

GCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAG  >  CP009273/992626‑992911
                                                                                                                                             |                                                                                                                                                
gcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTg                                                                                                                                            >  2:3269012/1‑148 (MQ=255)
 cagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGt                                                                                                                                           >  2:6807/1‑148 (MQ=255)
         ccACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTg                                                                                                                                   <  2:1487894/148‑1 (MQ=255)
           aCGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGtt                                                                                                                                 <  1:133875/148‑1 (MQ=255)
            cGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTc                                                                                                                                <  2:907529/148‑1 (MQ=255)
             ggTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTNCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTCGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCg                                                                                                                               <  1:1801597/148‑1 (MQ=255)
                aCCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTg                                                                                                                            >  1:2370129/1‑148 (MQ=255)
                 ccGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGa                                                                                                                           >  2:384443/1‑148 (MQ=255)
                    cccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAgg                                                                                                                        >  2:1815576/1‑148 (MQ=255)
                    cccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAgg                                                                                                                        >  2:2880142/1‑148 (MQ=255)
                      cGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTc                                                                                                                      <  2:1467790/148‑1 (MQ=255)
                         cGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTg                                                                                                                   <  2:570841/148‑1 (MQ=255)
                            cAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGt                                                                                                                >  2:61542/1‑148 (MQ=255)
                                ggTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCg                                                                                                            <  1:1343992/148‑1 (MQ=255)
                                  tAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                          <  2:254394/148‑1 (MQ=255)
                                  tAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                          <  2:790210/148‑1 (MQ=255)
                                  tAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGa                                                                                                          <  2:1953794/148‑1 (MQ=255)
                                    gctccacTTTGCCGTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTGTCAACGCACCGGTATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACt                                                                                                        <  2:2997573/145‑1 (MQ=255)
                                     caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTa                                                                                                       >  1:295397/1‑148 (MQ=255)
                                      accacGTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTAt                                                                                                      <  2:2342739/148‑1 (MQ=255)
                                         acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGa                                                                                                   >  2:2231494/1‑148 (MQ=255)
                                         acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGa                                                                                                   >  2:1313925/1‑148 (MQ=255)
                                          cTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGaa                                                                                                  <  1:3003737/148‑1 (MQ=255)
                                              gCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGa                                                                                              <  2:1666955/148‑1 (MQ=255)
                                                cTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGc                                                                                            >  2:2468464/1‑148 (MQ=255)
                                                 ttGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCa                                                                                           >  2:2427049/1‑148 (MQ=255)
                                                  tGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCAt                                                                                          <  1:2826927/148‑1 (MQ=255)
                                                  tGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCAt                                                                                          <  1:1051909/148‑1 (MQ=255)
                                                      aaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAg                                                                                      >  1:3539532/1‑148 (MQ=255)
                                                      aaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAg                                                                                      <  1:1912849/148‑1 (MQ=255)
                                                       aGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                     >  2:2379887/1‑148 (MQ=255)
                                                         cgcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAg                                                                                   <  1:747957/148‑1 (MQ=255)
                                                         cgcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAg                                                                                   <  1:270060/148‑1 (MQ=255)
                                                          gcgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGa                                                                                  <  1:1636159/148‑1 (MQ=255)
                                                             ttCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCt                                                                               >  1:3406773/1‑148 (MQ=255)
                                                                  gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                          >  2:1730279/1‑148 (MQ=255)
                                                                  gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                          >  2:1230717/1‑148 (MQ=255)
                                                                  gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                          >  1:1181050/1‑148 (MQ=255)
                                                                   cagaagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGACGGACTATCGAAACGAGCATAAAGTAGATCTTCCggg                                                                         <  1:1992709/148‑1 (MQ=255)
                                                                     gcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                       >  1:3285852/1‑148 (MQ=255)
                                                                         gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                   >  1:2991315/1‑148 (MQ=255)
                                                                         gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                   >  2:2136161/1‑148 (MQ=255)
                                                                           cGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtt                                                                 <  1:80279/148‑1 (MQ=255)
                                                                             aGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtttt                                                               >  1:419841/1‑148 (MQ=255)
                                                                              gCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTg                                                              >  1:3118364/1‑148 (MQ=255)
                                                                                    ttttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAtt                                                        >  1:486582/1‑148 (MQ=255)
                                                                                    ttttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAtt                                                        >  2:660521/1‑148 (MQ=255)
                                                                                     tttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTc                                                       >  1:2909341/1‑148 (MQ=255)
                                                                                      ttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTcc                                                      <  2:2633083/148‑1 (MQ=255)
                                                                                      ttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTcc                                                      <  2:2633169/148‑1 (MQ=255)
                                                                                      ttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTcc                                                      <  2:1355678/148‑1 (MQ=255)
                                                                                       tCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCa                                                     <  1:974895/148‑1 (MQ=255)
                                                                                       tCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCa                                                     <  1:501490/148‑1 (MQ=255)
                                                                                          aCGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTg                                                  <  2:2691178/148‑1 (MQ=255)
                                                                                           cGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGa                                                 >  1:2399873/1‑148 (MQ=255)
                                                                                           cGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGa                                                 >  2:2862126/1‑148 (MQ=255)
                                                                                            gCACCGGAATAGGCGGCTTTAGACACAGGCGTGGCGACAATCAGCCCATAGGCCTGGTGCAGCGGATCGCACAACGACTAGACTCCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGTAGTTTTGCAGATTCCAGTGAt                                                <  2:3388662/148‑1 (MQ=255)
                                                                                              aCCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaa                                              >  2:1338686/1‑148 (MQ=255)
                                                                                              aCCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaa                                              >  2:2203805/1‑148 (MQ=255)
                                                                                               ccGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaaa                                             >  2:3036758/1‑148 (MQ=255)
                                                                                                cGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAc                                            >  1:2624285/1‑148 (MQ=255)
                                                                                                 ggAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAAcc                                           >  1:2659531/1‑148 (MQ=255)
                                                                                                   aaTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctc                                         >  1:816506/1‑148 (MQ=255)
                                                                                                    aTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                        >  2:2626893/1‑148 (MQ=255)
                                                                                                    aTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAACATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctct                                        >  1:2368097/1‑148 (MQ=255)
                                                                                                     tAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTa                                       >  1:652689/1‑148 (MQ=255)
                                                                                                       ggcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACa                                     >  2:3222360/1‑148 (MQ=255)
                                                                                                        gcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACAt                                    >  1:3453406/1‑148 (MQ=255)
                                                                                                         cggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACAta                                   >  1:1976945/1‑147 (MQ=255)
                                                                                                         cggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATc                                   >  1:1116573/1‑148 (MQ=255)
                                                                                                         cggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATc                                   <  2:2370121/148‑1 (MQ=255)
                                                                                                          ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATcc                                  >  2:1552507/1‑148 (MQ=255)
                                                                                                           gcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCa                                 <  1:3219521/148‑1 (MQ=255)
                                                                                                            cTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                                >  2:2118680/1‑148 (MQ=255)
                                                                                                            cTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAg                                <  1:1815576/148‑1 (MQ=255)
                                                                                                              ttATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGc                              >  1:1687355/1‑148 (MQ=255)
                                                                                                                 taCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAt                           >  2:432034/1‑148 (MQ=255)
                                                                                                                 taCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAt                           >  2:431551/1‑148 (MQ=255)
                                                                                                                   cacaGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAttt                         <  1:3217507/148‑1 (MQ=255)
                                                                                                                    acaGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTa                        <  1:574343/148‑1 (MQ=255)
                                                                                                                     caGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAg                       >  1:2495752/1‑148 (MQ=255)
                                                                                                                      aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGt                      <  2:998693/148‑1 (MQ=255)
                                                                                                                      aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGt                      >  2:343354/1‑148 (MQ=255)
                                                                                                                      aGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGt                      <  1:3108191/148‑1 (MQ=255)
                                                                                                                        gCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                    >  2:3493001/1‑148 (MQ=255)
                                                                                                                        gCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                    >  1:94284/1‑148 (MQ=255)
                                                                                                                         cGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATATAGtttt                   >  2:2085963/1‑148 (MQ=255)
                                                                                                                                 cAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATGTAGTTTTTCCCgcgc           >  2:1092833/1‑148 (MQ=255)
                                                                                                                                  aaTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCa          >  1:2103308/1‑148 (MQ=255)
                                                                                                                                   aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat         >  1:1217019/1‑148 (MQ=255)
                                                                                                                                   aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat         >  1:1685182/1‑148 (MQ=255)
                                                                                                                                   aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat         <  1:2184254/148‑1 (MQ=255)
                                                                                                                                   aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat         <  1:2184016/148‑1 (MQ=255)
                                                                                                                                    tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGTCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAGTTAGTTTTTCCCGCGCata        >  1:2578063/1‑148 (MQ=255)
                                                                                                                                    tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata        <  1:2127402/148‑1 (MQ=255)
                                                                                                                                      aGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATAtt      >  2:1104889/1‑148 (MQ=255)
                                                                                                                                        cccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc    >  1:404370/1‑148 (MQ=255)
                                                                                                                                        cccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc    >  2:283755/1‑148 (MQ=255)
                                                                                                                                          cATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAg  >  2:307512/1‑148 (MQ=255)
                                                                                                                                          cATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATAGTCCAg  <  2:2772474/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAG  >  CP009273/992626‑992911

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