Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,908 T→G 25.5% T252P (ACC→CCC)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,782,9080TG25.5% 83.3 / 6.5 67T252P (ACC→CCC) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base T (18/27);  new base G (16/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 5.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCC  >  NC_000913/3782670‑3783142
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           
gCACGACCTAGTAAAGTCAATTCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                             <  2:23438/239‑1 (MQ=255)
gCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                             <  1:983328/239‑1 (MQ=255)
gCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                             <  1:387423/239‑1 (MQ=255)
         aGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTTCGCGCACGGTTGTCGGGACAGGTGAGCAccca                                                                                                                                                                                                                                    >  2:387423/1‑237 (MQ=255)
         aGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAg                                                                                                                                                                                                                                    <  2:667302/239‑1 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTTACCGACTGGGTGTCGGTGCAGGGAAGCACCaaaaagcca                                                                                                                                                                                                                             >  2:530998/1‑231 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGCCTGGTTGTTGGTACAGGGAAGCACCAAATCCCCa                                                                                                                                                                                                                             >  1:433923/1‑239 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGCTAAGCGACCGGGTGTTGGGACAGGTAAGCACCAAAACGCCa                                                                                                                                                                                                                             >  2:1206744/1‑239 (MQ=255)
                    tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTTGGGACCGGCAAGCACCAGAACGCCaaac                                                                                                                                                                                                                          >  1:82946/1‑236 (MQ=255)
                    tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTCCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGcc                                                                                                                                                                                                                         <  1:30353/239‑1 (MQ=255)
                         ccTCGCGCGCTTTTTTTCCGAAATATAATTCTTGATAAATTTCCTCGTTTATTGGCTTTTCACCGCCG‑AAGCGATGCGCCATTTCGCGTACTTCTTTCGTATTGCGGTAGCCTTCCCCCACCTGCCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGcc                                                                                                                                                                                                                    <  2:966700/239‑1 (MQ=255)
                         ccTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTTCGCGACTGGTTGTCCGTACAGGGAAGAACCAGAACGCCCAAcccccc                                                                                                                                                                                                                     >  1:501866/1‑236 (MQ=255)
                         ccTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTGTCACGCCTCGCTGTCGGCTCACGGCAGCACCTAACAGCCAAGcccccc                                                                                                                                                                                                                     >  2:273988/1‑236 (MQ=255)
                           ccgggcgcgTTTTTTCCGCAATATAATACTTGAGAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                  <  1:767203/235‑1 (MQ=255)
                              cgcgcgTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCCGAACGCCAAGCCCCGCCATga                                                                                                                                                                                                                >  1:97778/1‑237 (MQ=255)
                                                        ttGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCTACTGGTTGTTGGTACAGGTAAGCCCCAAAACGCCCAGCCCCGCCCTGCCCATAAAGGTGGCACGGTCGGCAcca                                                                                                                                                                                      >  1:557134/1‑237 (MQ=255)
                                                        ttGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCCGCTGGTTGTCGGTACAGGGAAGAACCAGATCGCCCAGCCCCGCCATGCCCCTAAAGGTGGCAGGGGTGGCAccc                                                                                                                                                                                      >  1:985111/1‑238 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTTCGCGGCTGGTTGTCGGGAAAGGGAAGGACCAAATCGCCAAACCCCGCCATGCCCACAAAGGGGGGAGGGTCGGCCCCCCgc                                                                                                                                                                                   >  1:186201/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGCCTGGTTGTCGGGACAGGGAAGAACCAGATCGCCAAGCCCCGCCCTGCCCATAAAAGTGGCAGGGTGGGCACCCCgc                                                                                                                                                                                   >  1:1161013/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCACAACGCCCAGCCCCCCCATGCCCCCAAAGGTGGCAGGGTCGGGACCCac                                                                                                                                                                                    >  2:637378/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAGCACCAGATCCCCCAGCCCCCCCCTGCCCATAAAGGGGGCCGGGGCGGGCCCCAgc                                                                                                                                                                                   >  1:882197/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAACACCAAAACCCCAAGCCCCCCCCTGCCCATAAAAGAGGGAGGGGCGGGACCCAg                                                                                                                                                                                    >  2:682614/1‑238 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTAAAGGTAAGCACCAAAACGCCAAACCCCGCCATGCCCCTAAAAGTGGCCGGGGCGGCACCCac                                                                                                                                                                                    >  1:743100/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACCGGTTGTCGGTACAAGTAAGCCCCAAATCCCCCAGCCCCCCCCCGCCCCTAAAAGTGGCAGGGGCGGCccccct                                                                                                                                                                                    >  1:281495/1‑236 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCCACTGGTTGTCGGTTCAGGTAAGGACCAAATCGCCAAACCCCCCCCTGCCCCTAAAGGGGGCAGGGGCGGCACCCAg                                                                                                                                                                                    >  2:732277/1‑238 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGGTACAGGTAAGCACCAGAACGCCCAACCCCGCCATTCCCCTAAAGGGGGGAGGGTCGGCATCCac                                                                                                                                                                                    >  1:711061/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACCCCGCTGGTTGTCGGGAAAGGTCAACACCAAAACCGCCAGCCCCCCCCAGCCCCTAAAGGTGGGCGGGGCGGCACCCac                                                                                                                                                                                    >  1:526344/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACCCCACTGGTTGTGTGTAAAAGGCAGAGCCAGAAAGCTAACCCCCCCCCAGCCCCTAAACGTCGTATGGCAGGCACCCaa                                                                                                                                                                                    >  2:1125273/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACACGACTTGTTTTCGGGACAGGGAAGCACCAAAACGCCAAACCCCCCCATGGCCCTAAAGGGGGCAGGGGCGGGACCCAgc                                                                                                                                                                                   >  1:427033/1‑239 (MQ=255)
                                                            tAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCCGCTGGTTGTCGGTACAGGTAAACACCAGAACGCCCAACCCCGCCCTGCCCATAAAGGGGGCAGGGGCGGCACCCAccg                                                                                                                                                                                  >  1:588266/1‑238 (MQ=255)
                                                            gaaaTTCCCGCGGTTAGT‑GCCTTTCAAACCCGAAAGTGAGGCCCCTGTTCGCTGACTTTTTTCTTGTTGCGGTCGCCCCCCCCCTCCTGCCCAGTCTGCTCCGGCGGGTTTTGTTC‑TCCCGTCCCTGGCCCAGGCCTCACGACAAACCGGCGTTTCCGCGCCGGGGTGCGGGTACAGGTAACCCCAAAATCACCACGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAgcgc                                                                                                                                                                                 <  2:248860/238‑1 (MQ=255)
                                                                         ttATTGGCATTTCACCGCCG‑AAGCGATGCGCCAGTCCGCGGACTTCTTTCGTATGGCGGTAGCCTTCCACCACCGGACCAATCTCCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                       <  2:214218/220‑1 (MQ=255)
                                                                                    tCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                         <  1:1103245/239‑1 (MQ=255)
                                                                                           cG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGGCCGG‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                                                                  <  1:1125273/239‑1 (MQ=255)
                                                                                           cG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGCCCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                                                                  <  1:530998/239‑1 (MQ=255)
                                                                                                      cgcCAGTTCGCGGACTTCTTTCGTATGGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCAGGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCa                                                                                                                                        <  1:1224288/239‑1 (MQ=255)
                                                                                                            ttCGCGGACTTTTTTCTTATTGCGGTAGCCTCCCCCCCCCTCACCATTCTTCTCCTGCGCTCTTTGTCCATCCTTGCCCTGCCCGA‑GCATCATGCCAAACCGGCGGTTACGCGCCTGGTTGTGGGTACAGGTAAGCCCCAGAGCGCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCg                                                                                                                                  <  2:108712/239‑1 (MQ=255)
                                                                                                              cgcgGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGCAAGCACCAAAACGCCAAACCCCCCCCCGCCCCTAAAGGGGGCAGGGGCGGGCCCCCACCCCGCCCCCAGACGCCGCATTTTCGCCAGCCCCCGGGGGGTCAGCCCCgg                                                                                                                                 >  2:419411/1‑237 (MQ=255)
                                                                                                                      tcttTCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                        <  1:1206744/239‑1 (MQ=255)
                                                                                                                      tcttTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                        <  1:732277/239‑1 (MQ=255)
                                                                                                                      tcttTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTAACCGAACGCAGGGCCGGACAGGCGAGCACCAGAACGGCAATCCCCGCCCTGCCCCCAACGGGGGCCGGGTCGCCATCCACCCCCCCCCCCACACGCGGCACATCTGCCAGCCCACCGGCGATCCGCGCCGCACGCGGAt                                                                                                                         >  2:241186/1‑238 (MQ=255)
                                                                                                                      tcttCCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCGCCGCGCTTTGTCCATCCATGCCCTGCCCGA‑CCATCATGCCAAAACGGCGGTTCCGCGACTGGTGGTCGGTACAGGTAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                        <  2:82946/239‑1 (MQ=255)
                                                                                                                            gTATGGCGGTAGCCTCCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACACCCAGGCCCTGACCGA‑GCATCATGCAAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACAAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                  <  2:1161013/239‑1 (MQ=255)
                                                                                                                             tATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                     <  1:735643/219‑1 (MQ=255)
                                                                                                                             tATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCCACCGGGTTTCGGTACAGGTAAACACCCAATCGCCAAGCCCCGCCCTGCCCCTAAAGGGGGCAGGGGCGGCACCCACCGCCCCCCCCCCACGCCACATTTCCGCCCGCCCACGGGGGCTCAGCGCCGCACGCGCCTTCGcccc                                                                                                                  >  2:110445/1‑238 (MQ=255)
                                                                                                                               ttGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTTACCGGCTTGGTGTCGGGACAGGGGAGCACCAAATCACCAAGCCCCCCCATGCCCCTAAAGGTGGCAGGGGCGGCCCCCCGCGCCCCCCCCACACGCCACATTTCAGCCCACCCCCCGGGGGTCCGCGCCCCACCCGCCTTCGCCCCcaa                                                                                                               >  2:920192/1‑239 (MQ=255)
                                                                                                                                     gcgCCTTCCCCCCCCTGACCACTCTCCTACTGCGCGCTTTGGACATCCGTGCCCTGCCCGA‑GCATCATGCCAAAACGGCGGTTTCGCGATTGTTGGCCGGTACGGAGAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGTTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGa                                                                                                          <  2:172381/236‑1 (MQ=255)
                                                                                                                                       gCCTCCCACCCCCTGACCAATCTTCCCCTGCGCGCTGTGTACATCCATGCCCTGACCGA‑GCACCAGGCAAAACCGGCGGTTCCGCGACTGGTTGTCGTTCCAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATAc                                                                                                       <  2:400869/239‑1 (MQ=255)
                                                                                                                                              ccccccTGCCCAACCTTCTCCTGCGCGCTTTGTACACCCATGCCCTGCCCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGa                                                                                               <  2:97778/235‑1 (MQ=255)
                                                                                                                                                    ggACCAATCTTCCCCCGCGCGCTTGGTACACCCATCCCCGGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATccc                                                                                         <  2:588266/238‑1 (MQ=255)
                                                                                                                                                      ccaaaTCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAACCGGCGGTTCCGCGACTGTTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCg                                                                                       <  1:682614/236‑1 (MQ=255)
                                                                                                                                                      aCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTTCGCGCCTGGTTTGCGGTACAGGGAAGCACCCGACCGCCAAGCCCCCCCCTGCCCCTAAAGGGGGCAGGGTCGGCACCCCGCGCCGCCCCCAGCCCCGCCATTTCCTCCCGCCCCCGGGCGGTCCCCGCCGCACGCGCATTCGCCCCAACACCGACACTGCC‑GGACCTcccc                                                                                        >  2:466317/1‑238 (MQ=255)
                                                                                                                                                             ttctCCTGCGCGCTTTGCACACCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAAt                                                                                <  1:449295/239‑1 (MQ=255)
                                                                                                                                                                  ctccgcgtTTTGTACACCCATCCCCTAACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACGGGTTGTCGGTGCAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCaa                                                                           <  2:957481/236‑1 (MQ=255)
                                                                                                                                                                     cgcgcTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTAAGCGACCGGTTGTCGGTACAGGGAAACACCAAAACGCCAAACCCCCCCCCGCCCCTAAAAGGGGGAGGGGCGGCACCCAGCGCCGCACCCAAACGCCACACTTCAGCCAGCCCCCCGGGGGTCAGCGCCGCGCGCGCATTCGCACCACAACCCACACCCTC‑CGGCAACCCCCCACCCACGGCAaaaa                                                                        >  2:438514/1‑239 (MQ=255)
                                                                                                                                                                                      tGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACTCGCCTTGTTGTCGGGACAAGGAAGCACCAAAACGCCAAGCCCCCCCCCGCCCATAAAGGTGGCAGGGGCGGCCCCCCCCGCCGCACCCACACGCGCCCTTTCAGACAGCCCCCGGGGGAACCACGCCCCACTCCCCCTCTCACCAACACCGCCCCCCCC‑CGACATCCCCCCCCCCACCGCAATAACCGTTTTCACCCcgcc                                                        >  2:474131/1‑238 (MQ=255)
                                                                                                                                                                                       gCCCTGACCGA‑GCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCaa                                                                           <  1:1186152/218‑1 (MQ=255)
                                                                                                                                                                                         ccTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCA‑‑GTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAg                                                  <  1:244001/239‑1 (MQ=255)
                                                                                                                                                                                             aCCGA‑GCATCATGCCAAAACGGCGGTTACGCGGCTGGGTGTCGGTACAGGTAAGGACCAGATCGCCAAGCCCCGCCCACCCCCTAAAAGGGGCAGGGGCGGGACCCCCCGCCGCACCAAAACGCGCCATTTCCGCCACCCCCCCGGGGGTAACAGCCGTCAGCGCCTTCGCACCCAAACCCATAACGCCGGGCCACCCCCCCCCCACCAGC‑ATAACGCTTTTCTCCCCgccgcccacc                                                 >  2:645411/1‑234 (MQ=255)
                                                                                                                                                                                                cagggacagtGCCAAAACGGCAGGTAGGCGAGTGGGTGGCGGTACGGGCAAGCACCAGCCCGCAAACCCCCGCCATGCCCATAAAGGTGGCAGGGGCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                               <  2:1066083/229‑1 (MQ=255)
                                                                                                                                                                                                             aaacCGGCGGTTACGCGACTGGTTGCCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAAt                                 <  1:637378/239‑1 (MQ=255)
                                                                                                                                                                                                                   gcggTTCCGCGACGGGTTGTCGGTACAGGGAAGCACCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGa                            <  2:557134/238‑1 (MQ=255)
                                                                                                                                                                                                                    gggTGACGCGACGGGTGGGCGGAACAGGTAACCCCCAGATCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGAt                           <  1:466317/237‑1 (MQ=255)
                                                                                                                                                                                                                    cgtttACGCGACTGGTTGTCGGTAAAGGTAAGCACCAGACCGCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGAt                           <  1:115899/235‑1 (MQ=255)
                                                                                                                                                                                                                              aCTGGGTGTCGGGACAGGGAAGCACCAGATCGCCAAGCCCCCCCATGCCCATAAAAGGGGGAGGGGCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGACAGCCCCCGGGTGATCAGCGCCGTACGCGCACTCGCACCAAAACCCAAACCGTC‑GGACAACCCCCCCCCCATGGCAATAACGCTTTTCACCGCCCCGCCAAGCTGCACCCCCAGTAAATCCGGATTGCTGGAAAc                 >  2:1000458/1‑238 (MQ=255)
                                                                                                                                                                                                                                       cGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtt       >  1:692224/1‑239 (MQ=255)
                                                                                                                                                                                                                                        ggTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACttt      >  1:663065/1‑239 (MQ=255)
                                                                                                                                                                                                                                        ggTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACttt      >  1:573915/1‑239 (MQ=255)
                                                                                                                                                                                                                                            cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc  >  2:133958/1‑239 (MQ=255)
                                                                                                                                                                                                                                             aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc  >  2:937412/1‑238 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           
GCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCC  >  NC_000913/3782670‑3783142

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