Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 877,980 G→T 5.3% Q71K (CAA→AAA)  deoR ← deoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273877,9800GT 5.3% 713.9 / 3.0 213Q71K (CAA→AAA) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (90/103);  new base T (18/2);  total (108/105)
Fisher's exact test for biased strand distribution p-value = 2.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

ACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAAC  >  CP009273/877832‑878128
                                                                                                                                                    |                                                                                                                                                    
aCGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTg                                                                                                                                                      >  2:3507594/1‑149 (MQ=255)
 cggcggTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGa                                                                                                                                                     <  2:3432399/149‑1 (MQ=255)
  ggcggTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGAt                                                                                                                                                    <  1:812490/149‑1 (MQ=255)
  ggcggTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGAt                                                                                                                                                    <  1:812517/149‑1 (MQ=255)
     ggTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTCTTCCACCAGGCGGGATTTTTGATCGc                                                                                                                                                 >  2:1351462/1‑149 (MQ=255)
       tAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTATTCCACCAGGCGGGATTTTTGATCGCtt                                                                                                                                               >  1:3663577/1‑149 (MQ=255)
        aaaaGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTa                                                                                                                                              <  2:4171232/149‑1 (MQ=255)
        aaaaGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTa                                                                                                                                              <  2:4172337/149‑1 (MQ=255)
         aaaGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTaa                                                                                                                                             >  2:110358/1‑149 (MQ=255)
             ggATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAgg                                                                                                                                         >  2:4327975/1‑149 (MQ=255)
             ggATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAgg                                                                                                                                         >  1:2437698/1‑149 (MQ=255)
              gATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCTGGAGGCGGGATTTTTGATCGCTTAACAGGt                                                                                                                                        <  1:2026342/149‑1 (MQ=255)
              gATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTATTACACCACGCGGGATTTTTTATAGCTTAACAGGt                                                                                                                                        >  1:2643915/1‑149 (MQ=255)
               aTTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTa                                                                                                                                       <  1:755215/149‑1 (MQ=255)
                tttCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTCTTAAACCAGGCGGGATTTTTTATAGCTTAACAGGTaa                                                                                                                                      >  2:2581591/1‑149 (MQ=255)
                 ttCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAt                                                                                                                                     <  1:2909445/149‑1 (MQ=255)
                     ttATCAATCGCTTCCATACTCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGGAAtggct                                                                                                                                 <  2:704229/149‑1 (MQ=255)
                     ttATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAtggct                                                                                                                                 <  1:1363498/149‑1 (MQ=255)
                        tCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTGTTACAACAGGCGGGATGTTTTATATCTTAAAACGTAAtttctgtc                                                                                                                              >  1:3001132/1‑142 (MQ=255)
                          aaTCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCAc                                                                                                                            <  1:1602657/149‑1 (MQ=255)
                          aaTCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAAAATGTAAAGGCTGGCAc                                                                                                                            >  2:2959672/1‑149 (MQ=255)
                           aTCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACt                                                                                                                           >  2:3980606/1‑149 (MQ=255)
                            tCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTg                                                                                                                          >  2:2737475/1‑149 (MQ=255)
                               cTTCAATCATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTgcgc                                                                                                                       <  2:2901614/149‑1 (MQ=255)
                               cTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTgcgc                                                                                                                       >  2:3582906/1‑149 (MQ=255)
                               cTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTgcgc                                                                                                                       >  2:3905657/1‑149 (MQ=255)
                                 tCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCgg                                                                                                                     >  2:1967376/1‑149 (MQ=255)
                                 tCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCgg                                                                                                                     <  2:1617015/149‑1 (MQ=255)
                                 tCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCgg                                                                                                                     >  2:2376250/1‑149 (MQ=255)
                                 tCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAAAAGGTAATGGCTGGCACTGCGCgg                                                                                                                     >  1:338690/1‑149 (MQ=255)
                                   aataatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATAGCTTAACAGGTAATGGCTGGCACTGCGCGGtt                                                                                                                   >  1:1692191/1‑149 (MQ=255)
                                    ataatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTc                                                                                                                  >  2:3170945/1‑149 (MQ=255)
                                    ataatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTc                                                                                                                  <  1:501911/149‑1 (MQ=255)
                                    ataatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTc                                                                                                                  >  2:2130281/1‑149 (MQ=255)
                                      aatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTGTTCCACCAGGCGGGATTTTTTTTCGCTTAACAGGGAATTGGGGGCACAGCGCCGTTGCa                                                                                                                >  1:122260/1‑149 (MQ=255)
                                      aatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCa                                                                                                                >  2:2310379/1‑149 (MQ=255)
                                      aatCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTATTCAAACAGGCGGGATTTTTTTTAGCATAAAAGGTAATTGCTGGCACATCCCGGTTcaa                                                                                                                >  1:122420/1‑147 (MQ=255)
                                        tCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGa                                                                                                              <  1:306313/149‑1 (MQ=255)
                                         ccACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTCTTCCCCCCGGCGGGGTTTTTTTTAGCTTAACAGGGAATGTGTGGCACTGCGCGGTTCCAGaa                                                                                                             >  1:908559/1‑149 (MQ=255)
                                            cGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACaa                                                                                                          <  1:297211/149‑1 (MQ=255)
                                             ggCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAat                                                                                                         <  1:1246049/149‑1 (MQ=255)
                                                gtggtggtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTTATCGCTTAACAGGTAATGGCTGGCACAGCGCGGTTCCAGAACAAtata                                                                                                      >  2:1948159/1‑149 (MQ=255)
                                                  ggtggtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAgc                                                                                                    <  1:3163224/149‑1 (MQ=255)
                                                  ggtggtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAgc                                                                                                    >  2:4342739/1‑149 (MQ=255)
                                                   gtggtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAgcc                                                                                                   <  2:352445/149‑1 (MQ=255)
                                                     ggtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAgccgc                                                                                                 <  1:3720589/149‑1 (MQ=255)
                                                     ggtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTCTTCCCCCAGGCGGGATTTTTGATCGCGTAACAGGGAAAGGCGGGCACAGCGCCGTTCCAGAACAAAAAAgcccc                                                                                                 >  2:1194/1‑147 (MQ=255)
                                                      gtgCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAgccgcc                                                                                                <  1:597367/149‑1 (MQ=255)
                                                              gTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCaaaa                                                                                        <  1:3833306/149‑1 (MQ=255)
                                                               tCAAAAAAGATGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCATCCCGGCGTTTTTGTTCCAACCGGCGGGATTTTTGATCGCTTAAAAGGTAATGGCTGGCACTGGGCGGGTACCGAACAATAGTGCAGACGAGCGAAac                                                                                       >  1:2060623/1‑149 (MQ=255)
                                                                 aaaaaaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCTGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCTCTGCGCGGCTCCTGATGAATATTGCCGCCGAGCAAAtcga                                                                                     <  2:3600012/149‑1 (MQ=255)
                                                                 aaaaaaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacga                                                                                     <  2:1024122/149‑1 (MQ=255)
                                                                 aaaaaaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacga                                                                                     <  1:1902009/149‑1 (MQ=255)
                                                                 aaaaaaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacga                                                                                     <  1:1701258/149‑1 (MQ=255)
                                                                   aaaaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTTATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAAAATAGCCGCCGAGCAAAacgacg                                                                                   >  1:4368077/1‑149 (MQ=255)
                                                                   aaaaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacgacg                                                                                   <  2:2486770/149‑1 (MQ=255)
                                                                     aaGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCATAAAAGGTAATGGCTGGCACAGCGCGGTTCCAGAACAATAAAGACGCCGAGCAAAACCACggg                                                                                 >  2:3051196/1‑149 (MQ=255)
                                                                       gagGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTGTTCCCCCCCGCGGGATTTTTGTTAGCTTAAAATGTAAAGGCGGGCACAGCGCGGTTCCACAACAAAATAGCCCCCGCGCCAAaccacccgcg                                                                               >  1:3458897/1‑144 (MQ=255)
                                                                        agGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTGTTACACCAGGCGGGATTTTTGATCGCGTAAAAGGTAAGGGCTGGCACTGCGCGGGTACAGAAAAATATATCCGCCGAGCAAAACCACCGgcgc                                                                              >  1:2741835/1‑149 (MQ=255)
                                                                         gggTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTGTTCCACCAGGCGGGATTTTTGATAGCTTAACAGGGAATTGCTGGCACAGCGCGGTTCCAGAACCATATAGCCGCCGAGCAAAAAGACGGGCGCa                                                                             >  1:2018254/1‑149 (MQ=255)
                                                                             cTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACtgt                                                                         <  1:3617066/149‑1 (MQ=255)
                                                                               gATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACCACGGGCGCACTGTgg                                                                       >  2:3880590/1‑149 (MQ=255)
                                                                                aTCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTGTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACACGTAATGGCTGGCACTGCGCGGTTACAGAAAAATATATCCGCCCAGCAAAACGACGGGCGCACTGTggg                                                                      >  2:819063/1‑148 (MQ=255)
                                                                                   gggTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGgttgt                                                                   <  2:1752619/149‑1 (MQ=255)
                                                                                       tCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGTATTTTTGATCGATTAACAGGTAATGGCTGGCACTGCGCGGTTCCATAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAg                                                               <  1:1918499/149‑1 (MQ=255)
                                                                                         tACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCCGGCGGGATTTTTGATCGCTTAACAGGTAATGGCGGGCACTGCGCGGTTCCAGAACAATATCGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGAt                                                             <  1:4109218/149‑1 (MQ=255)
                                                                                             aGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAACACGACGGGCGCACTGTGGTTGTTCAGATcgcg                                                         <  1:2169327/149‑1 (MQ=255)
                                                                                              gCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGcga                                                        <  2:3725086/149‑1 (MQ=255)
                                                                                                  cGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGTTTTTTGATCGATTAACTGGTAATGGCTGGCACTGCGCGGTACCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGaa                                                    <  1:130609/149‑1 (MQ=255)
                                                                                                        cAGTTTTGCCGCCCGGCGTTTTTGTTCCACCAGGCGGGATTTTTGATAGCTTAACAGGTAATGGCGGGCCCGGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAAAGACGGGCGCACAGTGGTTGTTCAGATCGCGACGAAACGTCa                                              >  1:688040/1‑149 (MQ=255)
                                                                                                        cAGTTTTGCCGCCCGGCGTTTTTCTTACACCAGGCGGGATTTTTTATAGCTTAACAGGTAATGGCTGGCACTGCGCGGTTACAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCa                                              >  2:1548539/1‑149 (MQ=255)
                                                                                                              tGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCg                                        <  2:1452560/149‑1 (MQ=255)
                                                                                                                  gccCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAc                                    <  2:859769/149‑1 (MQ=255)
                                                                                                                    ccGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAccc                                  <  2:2665670/149‑1 (MQ=255)
                                                                                                                    ccGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAccc                                  <  1:256059/149‑1 (MQ=255)
                                                                                                                      ggCGTTTTTCTTCCACCTGGCGGGATTTTTGATCGCTTAACAGGTTATGGCTGGCACTGCGCGGTTCCTGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCa                                <  2:4068611/149‑1 (MQ=255)
                                                                                                                      ggCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCa                                >  2:4304659/1‑149 (MQ=255)
                                                                                                                        cGTTTTTCTTCCACCAGGCGGGATTTTTTATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAg                              <  1:795895/149‑1 (MQ=255)
                                                                                                                         gTTTTTCTTCCACCAGGCGGTATTTTTGATCGCTTAACAGGTAAGGGCTGGCACTGCGCGGTTCCAGAACAATATCGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGc                             <  1:3406378/149‑1 (MQ=255)
                                                                                                                         gTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGc                             <  1:4353725/149‑1 (MQ=255)
                                                                                                                          tttttCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCa                            >  2:1070138/1‑149 (MQ=255)
                                                                                                                            tttCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAGGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAACACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAgg                          <  1:3400867/149‑1 (MQ=255)
                                                                                                                                ttcCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGggcgg                      <  1:1069795/149‑1 (MQ=255)
                                                                                                                                 tcCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGggcggc                     >  2:1169942/1‑149 (MQ=255)
                                                                                                                                       aGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGACGTc               <  2:764741/149‑1 (MQ=255)
                                                                                                                                        ggCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCt              >  1:474406/1‑149 (MQ=255)
                                                                                                                                        ggCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCt              >  2:1331207/1‑149 (MQ=255)
                                                                                                                                        ggCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGCTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCt              >  2:1331980/1‑149 (MQ=255)
                                                                                                                                         gCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCtt             >  1:1534128/1‑149 (MQ=255)
                                                                                                                                          cGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCttt            <  1:108794/149‑1 (MQ=255)
                                                                                                                                          cGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCttt            <  1:3985324/149‑1 (MQ=255)
                                                                                                                                          cGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCttt            <  1:2310787/149‑1 (MQ=255)
                                                                                                                                           gggATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTa           >  2:3104698/1‑149 (MQ=255)
                                                                                                                                           gggATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTa           >  2:4244568/1‑149 (MQ=255)
                                                                                                                                           gggATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTa           >  1:4411090/1‑149 (MQ=255)
                                                                                                                                           gggATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTa           >  1:1658008/1‑149 (MQ=255)
                                                                                                                                            ggATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTaa          >  1:34821/1‑149 (MQ=255)
                                                                                                                                                ttttGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATg      <  2:706913/149‑1 (MQ=255)
                                                                                                                                                 tttGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGt     <  1:907563/149‑1 (MQ=255)
                                                                                                                                                    gATCGCTTAACAGGTTATGGATGGCACTGCGCGTTTCAAGAAAAATATATCCTCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAAc  <  2:1651255/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGGCGGTAAAAGGGATTTCATTATCAATCGCTTCAATAATCCACGGCGTGGTGGTGCCACAGTCAAAAAAGAGGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAAC  >  CP009273/877832‑878128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: