Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,319 G→T 7.9% E234* (GAA→TAA)  gpmM → phosphoglycero mutase III, cofactor‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,779,3190GT 7.9% 786.2 / 4.2 281E234* (GAA→TAA) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base G (142/98);  new base T (3/38);  total (145/136)
Fisher's exact test for biased strand distribution p-value = 2.07e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

CATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGC  >  CP009273/3779171‑3779467
                                                                                                                                                    |                                                                                                                                                    
catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgatg                                                                                                                                                      <  2:3273325/149‑1 (MQ=255)
 atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                                                                     >  1:1266345/1‑149 (MQ=255)
 atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                                                                     >  2:2575273/1‑149 (MQ=255)
    ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAtt                                                                                                                                                  <  1:442064/149‑1 (MQ=255)
    ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAtt                                                                                                                                                  <  2:2312525/149‑1 (MQ=255)
     gTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTc                                                                                                                                                 <  1:3266098/149‑1 (MQ=255)
      tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAttag                                                                                                                                                >  1:637833/1‑147 (MQ=255)
        gTTACTTCGCGATGGACCGCGATTACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTg                                                                                                                                              <  1:181367/149‑1 (MQ=255)
        gTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTg                                                                                                                                              <  2:4857156/149‑1 (MQ=255)
        gTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTg                                                                                                                                              <  2:153087/149‑1 (MQ=255)
                  gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCg                                                                                                                                    >  2:2889707/1‑149 (MQ=255)
                  gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCg                                                                                                                                    >  2:2889006/1‑149 (MQ=255)
                   aTGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGt                                                                                                                                   >  2:2372963/1‑149 (MQ=255)
                       aCCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCGGCTTATTCTCGCGACGAAAAAGATGTATTCGTGAAAGCGACCGTTATc                                                                                                                               >  1:2609357/1‑149 (MQ=255)
                        ccGCTATAACCTTTGGGATCGCTTATAAAAAGCTAATGACCTGCTGACTCTGGGGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCTGGCTGTTTTTGCTCGCGACGAAATTTATTATTTCTTTAAAGCGACCGTTATcc                                                                                                                              <  2:2259647/149‑1 (MQ=255)
                            gATAACCGTTGGGATCGCGTAGAAAAAGCTTATAACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGc                                                                                                                          >  1:2118989/1‑149 (MQ=255)
                              tAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTg                                                                                                                        >  1:4888852/1‑149 (MQ=255)
                              tAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTg                                                                                                                        >  1:2280522/1‑149 (MQ=255)
                                  cGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAgg                                                                                                                    >  1:4208832/1‑149 (MQ=255)
                                     tGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCTGTATCCGTTCTGAAGGTCa                                                                                                                 >  2:5063534/1‑149 (MQ=255)
                                      gggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAg                                                                                                                >  2:952801/1‑149 (MQ=255)
                                       ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                               <  2:286144/149‑1 (MQ=255)
                                       ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGc                                                                                                               >  2:2801050/1‑149 (MQ=255)
                                       ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAATCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                               >  2:2801096/1‑149 (MQ=255)
                                        gATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGcc                                                                                                              >  2:1633616/1‑149 (MQ=255)
                                            gcgTAGAAAAAGCTTATGACCTGCTGACTCGGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGATTATACTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                                                                                          <  2:3609229/149‑1 (MQ=255)
                                               tAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                                                       >  2:1784401/1‑149 (MQ=255)
                                               tAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                                                       >  2:921858/1‑149 (MQ=255)
                                               tAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                                                       >  2:770218/1‑149 (MQ=255)
                                               tAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                                                       >  1:274044/1‑149 (MQ=255)
                                                  aaaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGTATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGCa                                                                                                    >  2:3471189/1‑149 (MQ=255)
                                                     aaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                                                                 >  2:2948047/1‑149 (MQ=255)
                                                          tATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTTTGCTCGCGACGAAAATGATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCAGCAATGGAAGa                                                                                            <  2:1470244/149‑1 (MQ=255)
                                                          tATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                                                            >  2:3330130/1‑149 (MQ=255)
                                                          tATGACCTGCTGACTCGGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTGGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATGCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCAGCAATGGAAGa                                                                                            <  2:4084377/149‑1 (MQ=255)
                                                           ttGACCTGCTTACTGTGGCGCTGGGCTTTTTCCCGGCCGATACGGCGGTGGCTGGTTTGCTGGCTGATTTTGCTCGCTACGAAATTTATTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGAc                                                                                           <  1:1559160/148‑1 (MQ=255)
                                                              aCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAAGGAAGACGGc                                                                                        >  1:3182287/1‑149 (MQ=255)
                                                              aCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATCCGTGCGGAAGGTCAGCCAGACGCGGCCATGGAAGACGGc                                                                                        >  1:592355/1‑149 (MQ=255)
                                                                 tgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                                                                     >  2:4246844/1‑149 (MQ=255)
                                                                 tgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                                                                     >  2:380096/1‑149 (MQ=255)
                                                                 ggctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGCCGGCGAt                                                                                     <  2:1432795/148‑1 (MQ=255)
                                                                  gctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTTCAGGCTGGGTATTCTCGCGACGAAAAAGATGAAATCGTGGAATCGACCGGTATACGTGCTGTATGTCAGCCACACCCGGCCATGGAAGAAGGCGaag                                                                                    >  1:3468986/1‑147 (MQ=255)
                                                                   ctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGc                                                                                   <  2:1173721/149‑1 (MQ=255)
                                                                   ctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAGGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGc                                                                                   <  2:4682983/149‑1 (MQ=255)
                                                                   ctgACTCGGGCGCTGGGCGAGTTCCAGGCCGATCCCGCCGTTGCGGGTTTGCCGGCTGCTTATGCTCGCGACTAAATTGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGc                                                                                   <  1:4712870/149‑1 (MQ=255)
                                                                    tgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCa                                                                                  >  2:1643130/1‑149 (MQ=255)
                                                                     gACTGTGGCGCTGGGCGTGTTCCTGGCCGATCCCGCGGTTGGGGTTTTGCAGGCTGTTTTTGCTCGCTACTAAATTTTTGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                                                                 <  1:3714518/149‑1 (MQ=255)
                                                                     gACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                                                                 >  2:3168329/1‑149 (MQ=255)
                                                                        tctGGCGCTGGGCGAGGTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGa                                                                              <  2:2754919/149‑1 (MQ=255)
                                                                              gctgggCGTGTTCCAGGCCGATACCGCGGTTGGTGGTTTGAAGGCTGCTTATGCTCGCGACTAAATTGATTAATTCTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                                                        <  1:1944075/146‑1 (MQ=255)
                                                                              gcAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTTCAGGCTGCGTATTCTCGCGACGAAAAATATGAAATAGTGAAATCGACCGGTATACCTGCTGGAGGTCAGCCACACGCGGCCAATGAAGACGGCGATGCACAGAttttaa                                                                        >  1:1719754/1‑147 (MQ=255)
                                                                               ctgggCGAGTACCAGGCCGATCCCGCCGTTGCTGGTTTGCTGGCTGCTTAGGCTCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCAt                                                                       <  2:1386727/147‑1 (MQ=255)
                                                                               cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCAt                                                                       <  2:3707297/149‑1 (MQ=255)
                                                                                aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTAGTGTAAGCCCACGATAACCGTGCGGTAGGTCAGCCAGACGCGGCACTGGAAGACGGCGATGCACTGATTTTCATg                                                                      >  2:4805339/1‑149 (MQ=255)
                                                                                  ggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                                                                    >  2:2504833/1‑149 (MQ=255)
                                                                                  ggCGAGTGCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGATGCTTATGCTCGCGACGAAAATGAGGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                                                                    <  2:4496229/149‑1 (MQ=255)
                                                                                      aGTTCCAGGCCGATACCGCCGTGGCTGGTTTGCAGGCTGCGTATGCGCGCGACGAAAAAGATGAATTAGTGAAAGCGACCGTTATACGTGCGGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACAGATTTTCATGAACATc                                                                >  1:1004860/1‑149 (MQ=255)
                                                                                       tttCCTGGCCTATACCGCGTTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACTAAAATTATGATTTCGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                                                               <  1:3221462/148‑1 (MQ=255)
                                                                                         tCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGt                                                             <  2:2996512/149‑1 (MQ=255)
                                                                                            tggCCTATCCCGCGGTTGGTGTTTTGCTGGCTGCTTTTGCGCGCGACTAAATTTTTTATTTGGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCt                                                          <  1:2331001/148‑1 (MQ=255)
                                                                                            aGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCt                                                          >  2:1200718/1‑149 (MQ=255)
                                                                                             ggCCTATACCGCCGTTGCTGGTTTGCCGGCTGATTTTGCTCGCGACTAAAATTATTAATGCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                         <  1:2410766/149‑1 (MQ=255)
                                                                                             ggCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                         <  2:2326367/149‑1 (MQ=255)
                                                                                             ggCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                         <  2:2576710/149‑1 (MQ=255)
                                                                                             ggCAGATACCGCGGTTGCTGTTTTGCAGGATGATTATGCTCGCTACGAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                         <  2:2576803/149‑1 (MQ=255)
                                                                                                   tACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCAAAGGAAGACGGCGATGCACTGATTTTCATGAACATCCGTGCTGACcgcg                                                   >  2:2200988/1‑149 (MQ=255)
                                                                                                       gccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcgcgcg                                               >  1:462154/1‑149 (MQ=255)
                                                                                                       gccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcgcgcg                                               >  1:2772318/1‑149 (MQ=255)
                                                                                                        gtgttgttgttTTGCTGGCTGTTTTTGATCGCTACTAAATTTTTGATTTGGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt                                              <  1:4776850/147‑1 (MQ=255)
                                                                                                        ccgTTGCTGGTTTGCTGTCTGCTTATGCTCGCGACTAAAATGATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt                                              <  1:2011261/149‑1 (MQ=255)
                                                                                                          gTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGa                                            >  2:3347231/1‑149 (MQ=255)
                                                                                                             gCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAt                                         >  1:932084/1‑149 (MQ=255)
                                                                                                                ggTTGGCTGGCTGCTTTTGCTCGCGCCGAAACTGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCAc                                      <  2:3920181/149‑1 (MQ=255)
                                                                                                                      ccggCTGCTTATGCTCGCGACTAAAATGATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGc                                <  1:5162611/147‑1 (MQ=255)
                                                                                                                      cAGGCTTCTATTGCTCGCTACTAAATTGATTATTTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGc                                <  1:1045367/149‑1 (MQ=255)
                                                                                                                         ggtgctTTTGCTCGCGACTAAATTGATGAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCttt                             <  1:4352134/147‑1 (MQ=255)
                                                                                                                         gctgctTAGGCTCGCGACGAAAATGATTATTTGGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCttt                             <  1:4116467/149‑1 (MQ=255)
                                                                                                                         gctgatTTTGCGCGCGACCAAAATTATTAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCttt                             <  1:126852/149‑1 (MQ=255)
                                                                                                                          ctgctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTc                            <  2:5002957/149‑1 (MQ=255)
                                                                                                                            gattttGCTCGCTACGAAATTTTTGAATTCTTTAAAGCGACCTTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGt                          <  2:2738166/144‑1 (MQ=255)
                                                                                                                              ttATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGa                        >  1:2501333/1‑149 (MQ=255)
                                                                                                                                aTGCTCGCGACTAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAAc                      <  2:3043076/149‑1 (MQ=255)
                                                                                                                                aTGCTCGCGACGAAAATGATGATTTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAAc                      <  2:451204/149‑1 (MQ=255)
                                                                                                                                    tCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTg                  >  1:896207/1‑149 (MQ=255)
                                                                                                                                      gcgACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAt                <  2:417538/149‑1 (MQ=255)
                                                                                                                                       ctacgaAAATGTTTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACACTGAtt               <  1:744603/147‑1 (MQ=255)
                                                                                                                                         acgaAAATGATTAATGCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTc             <  1:4083550/149‑1 (MQ=255)
                                                                                                                                         acgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTc             <  2:3237233/149‑1 (MQ=255)
                                                                                                                                           gaAATTGATTATTGCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa           <  1:3572313/149‑1 (MQ=255)
                                                                                                                                           gaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa           >  1:97275/1‑149 (MQ=255)
                                                                                                                                           gaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa           >  1:4376900/1‑149 (MQ=255)
                                                                                                                                              aaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATgg        >  2:4776103/1‑149 (MQ=255)
                                                                                                                                                ttatgaTTTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCt      <  2:4374497/147‑1 (MQ=255)
                                                                                                                                                tgatgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCt      >  1:3020075/1‑149 (MQ=255)
                                                                                                                                                 tttgaTTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCtt     <  1:2821204/147‑1 (MQ=255)
                                                                                                                                                 gatgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCtt     >  1:256542/1‑149 (MQ=255)
                                                                                                                                                  atgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTc    >  1:4528017/1‑149 (MQ=255)
                                                                                                                                                   tgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTcg   >  2:4970969/1‑149 (MQ=255)
                                                                                                                                                   tgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTcg   >  2:549631/1‑149 (MQ=255)
                                                                                                                                                   tgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTcg   >  1:3884503/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGC  >  CP009273/3779171‑3779467

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: