Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,316 G→T 9.6% D233Y (GAT→TAT)  gpmM → phosphoglycero mutase III, cofactor‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,779,3160GT 9.6% 572.2 / 4.3 222D233Y (GAT→TAT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base G (111/75);  new base T (4/32);  total (115/107)
Fisher's exact test for biased strand distribution p-value = 3.85e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

ATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTT  >  CP009273/3779169‑3779464
                                                                                                                                                   |                                                                                                                                                    
atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGGTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                                                                     >  1:1023283/1‑149 (MQ=255)
       gTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTc                                                                                                                                              <  2:1199864/149‑1 (MQ=255)
       gTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTc                                                                                                                                              >  2:193079/1‑149 (MQ=255)
         cGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGt                                                                                                                                            <  1:805184/149‑1 (MQ=255)
         cGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGt                                                                                                                                            <  1:1145774/149‑1 (MQ=255)
         cGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGt                                                                                                                                            <  1:1970881/149‑1 (MQ=255)
              ctacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                                                       >  1:938197/1‑149 (MQ=255)
                  gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGAc                                                                                                                                   >  1:548790/1‑149 (MQ=255)
                           cgcgATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCCGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCg                                                                                                                          <  2:2568881/149‑1 (MQ=255)
                            gcgATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGt                                                                                                                         >  1:1779791/1‑149 (MQ=255)
                              gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGCCCGTTATCCGTGc                                                                                                                       >  2:2614391/1‑149 (MQ=255)
                              gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGc                                                                                                                       >  2:1176954/1‑149 (MQ=255)
                               aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAACTGATGAATTCGTGACAGCGACCGTTATCCGTGCt                                                                                                                      >  2:2439084/1‑149 (MQ=255)
                               aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCt                                                                                                                      >  2:4319858/1‑149 (MQ=255)
                               aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCt                                                                                                                      >  2:4319683/1‑149 (MQ=255)
                                tAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTg                                                                                                                     >  2:2739895/1‑149 (MQ=255)
                                 aaCCGTTGGGAGCGCGTAGAAAACGCGTTGGACCTGCTGACTCTGGCGCCGGGCGTGTTCCAGGCCGATCCCGCCGTGGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAGGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGa                                                                                                                    <  2:4266402/149‑1 (MQ=255)
                                  aCCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGaa                                                                                                                   >  2:1309888/1‑149 (MQ=255)
                                   ccGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAg                                                                                                                  >  2:3432047/1‑149 (MQ=255)
                                      ttGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTc                                                                                                               <  2:1763489/149‑1 (MQ=255)
                                       tGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCTGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCa                                                                                                              <  2:1593148/149‑1 (MQ=255)
                                         gtaTCGCTTATAAAAAGCTTATGACCTGCTGACTCTGGCGCGGGGCGAGTGCCTGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTTTGCTCGCGACTAAAATTTTTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                            <  2:3772575/147‑1 (MQ=255)
                                         ggATCGCGTATAAAATGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCCGGCCGATACCGCCGTGGCTGGTTGGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                            <  2:2336918/149‑1 (MQ=255)
                                         ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                            >  2:679771/1‑149 (MQ=255)
                                          gATCGCGTAGAAAAATCTTATGACCTGCTGACTCTGGCGCTGGGCGTGTTCCTGGCCTATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGcc                                                                                                           <  1:1743155/149‑1 (MQ=255)
                                           aTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCCGGCTGCGTATGCTCGCGACGAAAATGATTAATTAGTGAAATCGACCGGTATACGTGCAGAAGGTCAGCCa                                                                                                          >  1:3020056/1‑149 (MQ=255)
                                           aTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                                                          >  2:3615236/1‑149 (MQ=255)
                                              gcgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCCTTATACGTGCTGAAGGTCAGCCCGAc                                                                                                       >  1:1948509/1‑149 (MQ=255)
                                                gTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgc                                                                                                     >  2:1328777/1‑149 (MQ=255)
                                                  aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                                                                   >  2:1638215/1‑149 (MQ=255)
                                                  aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                                                                   >  2:371525/1‑149 (MQ=255)
                                                    aaaaaGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCGGTGGCTGGTTTGCAGGCTGCTTTTGCTCGCGACGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCa                                                                                                 <  2:1869665/149‑1 (MQ=255)
                                                    aaaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCa                                                                                                 >  1:2541562/1‑149 (MQ=255)
                                                     aaaaGTTTTTGACGTGCTGACTCGGGCGCAGGGCGAGTTCCAGGCCTATACCGCCGTTGCTGTTTTGCTGGCTGTTTTTGCTCGCGACTAAATTTATTAATTGGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                                                                <  2:3347611/149‑1 (MQ=255)
                                                     aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAAATCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCCGACGCGGCaa                                                                                                >  1:2677929/1‑149 (MQ=255)
                                                     aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCaa                                                                                                >  1:2677783/1‑149 (MQ=255)
                                                       aaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                                                              >  2:4156515/1‑149 (MQ=255)
                                                        aGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGccctgg                                                                                             >  1:3142830/1‑144 (MQ=255)
                                                           ttATGACCTGCTGACTCTGGCGCAGGGCGAGTGCCCGGCCGATACCGCGGTTGGTGGTTTGCAGGCTGCTTAGGCGCGCGACTAAAATGTTTAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                                                          <  1:2297379/149‑1 (MQ=255)
                                                            tATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                                                         >  1:995713/1‑149 (MQ=255)
                                                            tATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                                                         >  2:4123512/1‑149 (MQ=255)
                                                             ttcaCCGGCTGACTGGGGCGCTGGGCGTGTTCCAGGCCTATACCGCGGTTGGTGTTTTGCTGGCTTATTATGCTCGCCAAGAAAATTATTATTTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGAc                                                                                        <  1:1460370/146‑1 (MQ=255)
                                                              tGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGATTTTGCTCGCGACGAAAATTATGATTTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                                                       <  2:3044749/149‑1 (MQ=255)
                                                              tGACCTGCTGACTCGGGCGCTGGGCGAGTTCCCGGCCGATACCGCGTTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAATTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                                                       <  2:3121263/149‑1 (MQ=255)
                                                                 ccTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg                                                                                    <  2:2858486/149‑1 (MQ=255)
                                                                   tgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                                                                  <  2:3960773/149‑1 (MQ=255)
                                                                      tgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCa                                                                               >  2:1037661/1‑149 (MQ=255)
                                                                       gACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                                                              >  1:1058092/1‑149 (MQ=255)
                                                                       gACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                                                              >  1:3509833/1‑149 (MQ=255)
                                                                             ggCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAttt                                                                        >  2:3511318/1‑149 (MQ=255)
                                                                                gctgggCGTGTGCCCGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATTATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                                                     <  1:2473002/146‑1 (MQ=255)
                                                                                gcAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                                                     <  2:1993072/149‑1 (MQ=255)
                                                                                gcAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                                                     <  2:1992116/149‑1 (MQ=255)
                                                                                   gggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTTCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTAATGa                                                                  >  1:3142185/1‑149 (MQ=255)
                                                                                   gggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                                                                  >  2:1253332/1‑149 (MQ=255)
                                                                                      cgtgTTCCAGGCCTATCCCGCGTTTGTTTGTTTGCCGGCTGCTTATGCTCCCTACTAAATTGTTGATTCCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt                                                               <  1:840655/146‑1 (MQ=255)
                                                                                      cGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGATTATGCTCGCGACGAAAATGTTGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt                                                               <  2:1044047/149‑1 (MQ=255)
                                                                                      cGAGTGCCAGGCCGATACGGCGGTTGGTGGTTTGCTGGCGGTTTTTGCTCGCTACTAAATTTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt                                                               <  1:3779998/149‑1 (MQ=255)
                                                                                         gtgcctggCCGTTACCGCGTTTGTTGGTTGGCTGGTTGCTTTTGCTGGCGACTAAATTTATTATTTGGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                                                            <  1:3349088/143‑1 (MQ=255)
                                                                                           tCCAGGCCTATACCGCCGTGGGTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATGATGATTTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGt                                                          <  1:4393391/149‑1 (MQ=255)
                                                                                            ccAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTg                                                         >  2:919778/1‑149 (MQ=255)
                                                                                               ggCCGATACCGGGGTTGGTGTTGTGCTGGGGTCTTAGGCTCGCGACGAAAATTTTTAATAGGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                      <  1:3389630/149‑1 (MQ=255)
                                                                                               ggCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCAAGAACTTCCGTGCTg                                                      >  1:734216/1‑149 (MQ=255)
                                                                                                gCCGATCCCGCGGTTGCTGGTTTGCAGGCTGTTTTTGCGCGCTACTAAATTGATTATTTCGAGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGa                                                     <  2:4059714/149‑1 (MQ=255)
                                                                                                 ccGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAc                                                    <  2:736513/149‑1 (MQ=255)
                                                                                                 cagATCCCGCGTTTGGTGTTTTGCGGGCTGGTTATGCTCGCTACTAAAATTATGATTTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAc                                                    <  1:165380/147‑1 (MQ=255)
                                                                                                       ccgccgTTGCTGGTTTTCAGGCTGCTTATTCTCGCGACGAAAAAGATGAAATAGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTAATTAACTTACGTGCTGACcgcgcg                                              >  1:4321971/1‑149 (MQ=255)
                                                                                                               gCTGTTTTGCTGGCTTCTTATGCTCGCGACTAAATTGTTGATTTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAt                                      <  1:1572042/149‑1 (MQ=255)
                                                                                                               gCTGGTTGGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAt                                      <  2:1546698/149‑1 (MQ=255)
                                                                                                                cTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATc                                     <  2:4045157/149‑1 (MQ=255)
                                                                                                                 tGGTTTGCTGGGTGATTTGGCGCGCGACTAAATTTATGAATTGGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCa                                    <  1:3812204/149‑1 (MQ=255)
                                                                                                                 gggTTTGCTGGCTGATTATGCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCa                                    <  1:3628163/148‑1 (MQ=255)
                                                                                                                  ggTTTGCTGGCTGCTTATGCTCGCGACTAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCAc                                   <  1:3653921/149‑1 (MQ=255)
                                                                                                                       gCAGGCTGGTTATGCTCGCGACGAAATTTTTGAATTCTTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTg                              <  2:1146987/149‑1 (MQ=255)
                                                                                                                        cAGGCTGCTTTTGCTCGCGACTAAAATGTTTAATTCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGc                             <  1:555237/149‑1 (MQ=255)
                                                                                                                           gctgctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCttt                          >  1:4281190/1‑149 (MQ=255)
                                                                                                                              gctTATGCTCGCGGCGAAAATGATGAATTCGTGAAAGCGACCGTGATCCGTGCTGTAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACTGATTTTCTGGAACTTCCGTGCTGACCGCGCGCGTGACATCACTCGGGCTTTCGt                       >  1:2796832/1‑149 (MQ=255)
                                                                                                                                ttATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGa                     >  1:36805/1‑149 (MQ=255)
                                                                                                                                ttATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGa                     <  2:2342215/149‑1 (MQ=255)
                                                                                                                                   tGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACg                  >  1:1274485/1‑149 (MQ=255)
                                                                                                                                   tGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACg                  >  1:2106427/1‑149 (MQ=255)
                                                                                                                                    gCTCGCGACGAAAATGATGAATTCTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGc                 <  1:1811585/149‑1 (MQ=255)
                                                                                                                                     cTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCt                >  2:3786847/1‑149 (MQ=255)
                                                                                                                                      tCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTg               >  2:2926287/1‑149 (MQ=255)
                                                                                                                                        gcgACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAt             >  1:3829914/1‑149 (MQ=255)
                                                                                                                                          gactaAATTTATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAttt           <  1:3238779/149‑1 (MQ=255)
                                                                                                                                          gacgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAttt           <  1:3595929/149‑1 (MQ=255)
                                                                                                                                           acgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTc          >  2:19044/1‑149 (MQ=255)
                                                                                                                                           acgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTc          >  2:4033665/1‑149 (MQ=255)
                                                                                                                                             gaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa        <  1:3293169/149‑1 (MQ=255)
                                                                                                                                             gaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa        >  2:1271992/1‑149 (MQ=255)
                                                                                                                                              aaaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGAt       <  2:2648738/149‑1 (MQ=255)
                                                                                                                                              aaaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGAt       <  1:1007923/149‑1 (MQ=255)
                                                                                                                                              aaaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGAt       >  1:1129291/1‑149 (MQ=255)
                                                                                                                                               aaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATg      >  1:1576151/1‑149 (MQ=255)
                                                                                                                                                aaTGATTAATTCGTGAAAGCGACCGTTATCAGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGATGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATgg     >  2:1454225/1‑149 (MQ=255)
                                                                                                                                                 atgatgAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGc    >  1:921360/1‑149 (MQ=255)
                                                                                                                                                   gattaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCtt  <  1:948691/149‑1 (MQ=255)
                                                                                                                                                   gatgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCtt  >  1:4075993/1‑149 (MQ=255)
                                                                                                                                                   gatgaATTCGGGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCtt  <  2:210613/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTT  >  CP009273/3779169‑3779464

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