Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,907 T→G 27.1% 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,9070TG27.1% 84.8 / 8.9 74T252P (ACC→CCC) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base G (19/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 1.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCC  >  NC_000913/3782669‑3783141
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
gCACGACCTAGTAAAGTCAATTCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                              <  2:23438/239‑1 (MQ=255)
gCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                              <  1:983328/239‑1 (MQ=255)
gCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                              <  1:387423/239‑1 (MQ=255)
         aGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGCACGGTTGTCGGGACAGGTGAGCAccca                                                                                                                                                                                                                                     >  2:387423/1‑237 (MQ=255)
         aGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAg                                                                                                                                                                                                                                     <  2:667302/239‑1 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTACCGACTGGGTGTCGGTGCAGGGAAGCACCaaaaagcca                                                                                                                                                                                                                              >  2:530998/1‑231 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTGGTTGTTGGTACAGGGAAGCACCAAATCCCCa                                                                                                                                                                                                                              >  1:433923/1‑239 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGCTAAGCGACCGGGTGTTGGGACAGGTAAGCACCAAAACGCCa                                                                                                                                                                                                                              >  2:1206744/1‑239 (MQ=255)
                    tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTTGGGACCGGCAAGCACCAGAACGCCaaac                                                                                                                                                                                                                           >  1:82946/1‑236 (MQ=255)
                    tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTCCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGcc                                                                                                                                                                                                                          <  1:30353/239‑1 (MQ=255)
                         ccTCGCGCGCTTTTTTTCCGAAATATAATTCTTGATAAATTTCCTCGTTTATTGGCTTTTCACCGCCGAAGCGATGCGCCATTTCGCGTACTTCTTTCGTATTGCGGTAGCCTTCCCCCACCTGCCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGcc                                                                                                                                                                                                                     <  2:966700/239‑1 (MQ=255)
                         ccTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGACTGGTTGTCCGTACAGGGAAGAACCAGAACGCCCAAcccccc                                                                                                                                                                                                                      >  1:501866/1‑236 (MQ=255)
                         ccTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTGTCACGCCTCGCTGTCGGCTCACGGCAGCACCTAACAGCCAAGcccccc                                                                                                                                                                                                                      >  2:273988/1‑236 (MQ=255)
                           ccgggcgcgTTTTTTCCGCAATATAATACTTGAGAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                   <  1:767203/235‑1 (MQ=255)
                              cgcgcgTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCCGAACGCCAAGCCCCGCCATga                                                                                                                                                                                                                 >  1:97778/1‑237 (MQ=255)
                                                        ttGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCTACTGGTTGTTGGTACAGGTAAGCCCCAAAACGCCCAGCCCCGCCCTGCCCATAAAGGTGGCACGGTCGGCAcca                                                                                                                                                                                       >  1:557134/1‑237 (MQ=255)
                                                        ttGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCGCTGGTTGTCGGTACAGGGAAGAACCAGATCGCCCAGCCCCGCCATGCCCCTAAAGGTGGCAGGGGTGGCAccc                                                                                                                                                                                       >  1:985111/1‑238 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGGCTGGTTGTCGGGAAAGGGAAGGACCAAATCGCCAAACCCCGCCATGCCCACAAAGGGGGGAGGGTCGGCCCCCCgc                                                                                                                                                                                    >  1:186201/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTGGTTGTCGGGACAGGGAAGAACCAGATCGCCAAGCCCCGCCCTGCCCATAAAAGTGGCAGGGTGGGCACCCCgc                                                                                                                                                                                    >  1:1161013/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCACAACGCCCAGCCCCCCCATGCCCCCAAAGGTGGCAGGGTCGGGACCCac                                                                                                                                                                                     >  2:637378/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAGCACCAGATCCCCCAGCCCCCCCCTGCCCATAAAGGGGGCCGGGGCGGGCCCCAgc                                                                                                                                                                                    >  1:882197/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAACACCAAAACCCCAAGCCCCCCCCTGCCCATAAAAGAGGGAGGGGCGGGACCCAg                                                                                                                                                                                     >  2:682614/1‑238 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTAAAGGTAAGCACCAAAACGCCAAACCCCGCCATGCCCCTAAAAGTGGCCGGGGCGGCACCCac                                                                                                                                                                                     >  1:743100/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGTTGTCGGTACAAGTAAGCCCCAAATCCCCCAGCCCCCCCCCGCCCCTAAAAGTGGCAGGGGCGGCccccct                                                                                                                                                                                     >  1:281495/1‑236 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACTGGTTGTCGGTTCAGGTAAGGACCAAATCGCCAAACCCCCCCCTGCCCCTAAAGGGGGCAGGGGCGGCACCCAg                                                                                                                                                                                     >  2:732277/1‑238 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGGTACAGGTAAGCACCAGAACGCCCAACCCCGCCATTCCCCTAAAGGGGGGAGGGTCGGCATCCac                                                                                                                                                                                     >  1:711061/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCCGCTGGTTGTCGGGAAAGGTCAACACCAAAACCGCCAGCCCCCCCCAGCCCCTAAAGGTGGGCGGGGCGGCACCCac                                                                                                                                                                                     >  1:526344/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCCACTGGTTGTGTGTAAAAGGCAGAGCCAGAAAGCTAACCCCCCCCCAGCCCCTAAACGTCGTATGGCAGGCACCCaa                                                                                                                                                                                     >  2:1125273/1‑237 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACACGACTTGTTTTCGGGACAGGGAAGCACCAAAACGCCAAACCCCCCCATGGCCCTAAAGGGGGCAGGGGCGGGACCCAgc                                                                                                                                                                                    >  1:427033/1‑239 (MQ=255)
                                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAACAGACTGGTTGTCGGTTAAGGGAAGCACCAGAACGCCAAACCCCCCCCAGCCCCTAAAGgggggcggcggcaccccgac                                                                                                                                                                                     >  2:244001/1‑222 (MQ=255)
                                                            tAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCGCTGGTTGTCGGTACAGGTAAACACCAGAACGCCCAACCCCGCCCTGCCCATAAAGGGGGCAGGGGCGGCACCCAccg                                                                                                                                                                                   >  1:588266/1‑238 (MQ=255)
                                                            gaaaTTCCCGCGGTTAGTGCCTTTCAAACCCGAAAGTGAGGCCCCTGTTCGCTGACTTTTTTCTTGTTGCGGTCGCCCCCCCCCTCCTGCCCAGTCTGCTCCGGCGGGTTTTGTTCTCCCGTCCCTGGCCCAGGCCTCACGACAAACCGGCGTTTCCGCGCCGGGGTGCGGGTACAGGTAACCCCAAAATCACCACGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAgcgc                                                                                                                                                                                  <  2:248860/238‑1 (MQ=255)
                                                                         ttATTGGCATTTCACCGCCGAAGCGATGCGCCAGTCCGCGGACTTCTTTCGTATGGCGGTAGCCTTCCACCACCGGACCAATCTCCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                        <  2:214218/220‑1 (MQ=255)
                                                                                   ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGGAAGCACCAGATCGCCAAACCCCCCCCCGCCCCTAAAGGTGGCAGGGGCGGGACCCAGcgccccccctggctcttttacacat                                                                                                                                                            >  1:214218/1‑217 (MQ=255)
                                                                                    tCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                         <  1:1103245/239‑1 (MQ=255)
                                                                                           cGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGGCCGGGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                                                                  <  1:1125273/239‑1 (MQ=255)
                                                                                           cGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGCCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa                                                                                                                                                  <  1:530998/239‑1 (MQ=255)
                                                                                                      cgcCAGTTCGCGGACTTCTTTCGTATGGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCa                                                                                                                                        <  1:1224288/239‑1 (MQ=255)
                                                                                                            ttCGCGGACTTTTTTCTTATTGCGGTAGCCTCCCCCCCCCTCACCATTCTTCTCCTGCGCTCTTTGTCCATCCTTGCCCTGCCCGAGCATCATGCCAAACCGGCGGTTACGCGCCTGGTTGTGGGTACAGGTAAGCCCCAGAGCGCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCg                                                                                                                                  <  2:108712/239‑1 (MQ=255)
                                                                                                              cgcgGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGCAAGCACCAAAACGCCAAACCCCCCCCCGCCCCTAAAGGGGGCAGGGGCGGGCCCCCACCCCGCCCCCAGACGCCGCATTTTCGCCAGCCCCCGGGGGGTCAGCCCCgg                                                                                                                                 >  2:419411/1‑237 (MQ=255)
                                                                                                                      tcttTCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                        <  1:1206744/239‑1 (MQ=255)
                                                                                                                      tcttTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                        <  1:732277/239‑1 (MQ=255)
                                                                                                                      tcttTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAACCGAACGCAGGGCCGGACAGGCGAGCACCAGAACGGCAATCCCCGCCCTGCCCCCAACGGGGGCCGGGTCGCCATCCACCCCCCCCCCCACACGCGGCACATCTGCCAGCCCACCGGCGATCCGCGCCGCACGCGGAt                                                                                                                         >  2:241186/1‑238 (MQ=255)
                                                                                                                      tcttCCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCGCCGCGCTTTGTCCATCCATGCCCTGCCCGACCATCATGCCAAAACGGCGGTTCCGCGACTGGTGGTCGGTACAGGTAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                        <  2:82946/239‑1 (MQ=255)
                                                                                                                            gTATGGCGGTAGCCTCCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACACCCAGGCCCTGACCGAGCATCATGCAAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACAAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                  <  2:1161013/239‑1 (MQ=255)
                                                                                                                             tATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                     <  1:735643/219‑1 (MQ=255)
                                                                                                                             tATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACCGGGTTTCGGTACAGGTAAACACCCAATCGCCAAGCCCCGCCCTGCCCCTAAAGGGGGCAGGGGCGGCACCCACCGCCCCCCCCCCACGCCACATTTCCGCCCGCCCACGGGGGCTCAGCGCCGCACGCGCCTTCGcccc                                                                                                                  >  2:110445/1‑238 (MQ=255)
                                                                                                                               ttGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTACCGGCTTGGTGTCGGGACAGGGGAGCACCAAATCACCAAGCCCCCCCATGCCCCTAAAGGTGGCAGGGGCGGCCCCCCGCGCCCCCCCCACACGCCACATTTCAGCCCACCCCCCGGGGGTCCGCGCCCCACCCGCCTTCGCCCCcaa                                                                                                               >  2:920192/1‑239 (MQ=255)
                                                                                                                                  cGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGCCTGGTTGTCGGTACAGGTAAGCACCAAAACGCCCAGACCCCCCCCGCCCCTAAAAGGGGCAGGGGCGGGCCCCACCGCCGCCCCCAAACCCGGCCATTTAGCCCGCCCCCGCGGGCGCAGCGCCgccaccgccttcccccccaaac                                                                                                             >  2:579834/1‑217 (MQ=255)
                                                                                                                                     gcgCCTTCCCCCCCCTGACCACTCTCCTACTGCGCGCTTTGGACATCCGTGCCCTGCCCGAGCATCATGCCAAAACGGCGGTTTCGCGATTGTTGGCCGGTACGGAGAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGTTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGa                                                                                                          <  2:172381/236‑1 (MQ=255)
                                                                                                                                       gCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTGGGTGTCGGTAAAGGTAAGAACCAAAACGCCAAACCCCCCCCTGCCCATAAAGGTGGGCGGGGGGGCCCCCCACGCCGCCCCCAGAACCGGCACTTCCATCAGCCCCCGGGGGATACGCGCCGGCCgcgccctgcctcttataaacacc                                                                                                        >  2:735643/1‑219 (MQ=255)
                                                                                                                                       gCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTTGTTTTCGGTACAGGTAAGCACCACATCGCCCAGACCCCCCCAGCCCATTAAAGGGGCAGGGGGCGCCCCCCGCCCCGCCCCCAGCAGCGCCACTTCCGCCCAGCCCCGGGGGGTCAGCCCCGGACacgccactcccccccaacccaca                                                                                                        >  2:688544/1‑219 (MQ=255)
                                                                                                                                       gCCTCCCACCCCCTGACCAATCTTCCCCTGCGCGCTGTGTACATCCATGCCCTGACCGAGCACCAGGCAAAACCGGCGGTTCCGCGACTGGTTGTCGTTCCAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATAc                                                                                                       <  2:400869/239‑1 (MQ=255)
                                                                                                                                              ccccccTGCCCAACCTTCTCCTGCGCGCTTTGTACACCCATGCCCTGCCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                                                               <  2:97778/235‑1 (MQ=255)
                                                                                                                                                    ggACCAATCTTCCCCCGCGCGCTTGGTACACCCATCCCCGGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATccc                                                                                         <  2:588266/238‑1 (MQ=255)
                                                                                                                                                      ccaaaTCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTCCGCGACTGTTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCg                                                                                       <  1:682614/236‑1 (MQ=255)
                                                                                                                                                      aCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCGCCTGGTTTGCGGTACAGGGAAGCACCCGACCGCCAAGCCCCCCCCTGCCCCTAAAGGGGGCAGGGTCGGCACCCCGCGCCGCCCCCAGCCCCGCCATTTCCTCCCGCCCCCGGGCGGTCCCCGCCGCACGCGCATTCGCCCCAACACCGACACTGCCGGACCTcccc                                                                                        >  2:466317/1‑238 (MQ=255)
                                                                                                                                                             ttctCCTGCGCGCTTTGCACACCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAAt                                                                                <  1:449295/239‑1 (MQ=255)
                                                                                                                                                                  ctccgcgtTTTGTACACCCATCCCCTAACCGAGCATCATGCCAAAACGGCGGTTACGCGACGGGTTGTCGGTGCAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCaa                                                                           <  2:957481/236‑1 (MQ=255)
                                                                                                                                                                     cgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGACCGGTTGTCGGTACAGGGAAACACCAAAACGCCAAACCCCCCCCCGCCCCTAAAAGGGGGAGGGGCGGCACCCAGCGCCGCACCCAAACGCCACACTTCAGCCAGCCCCCCGGGGGTCAGCGCCGCGCGCGCATTCGCACCACAACCCACACCCTCCGGCAACCCCCCACCCACGGCAaaaa                                                                        >  2:438514/1‑239 (MQ=255)
                                                                                                                                                                                 acccttccccttcccGAGCACCATGCCAAAACGGCGTTTCCGCGACGGTTTGTCGGCACGGGTAAGCACCACTCCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcg                                                            <  2:882197/226‑1 (MQ=255)
                                                                                                                                                                                      tGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACTCGCCTTGTTGTCGGGACAAGGAAGCACCAAAACGCCAAGCCCCCCCCCGCCCATAAAGGTGGCAGGGGCGGCCCCCCCCGCCGCACCCACACGCGCCCTTTCAGACAGCCCCCGGGGGAACCACGCCCCACTCCCCCTCTCACCAACACCGCCCCCCCCCGACATCCCCCCCCCCACCGCAATAACCGTTTTCACCCcgcc                                                        >  2:474131/1‑238 (MQ=255)
                                                                                                                                                                                       gCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCaa                                                                           <  1:1186152/218‑1 (MQ=255)
                                                                                                                                                                                         ccTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCA‑‑GTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAg                                                  <  1:244001/239‑1 (MQ=255)
                                                                                                                                                                                             aCCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGGTGTCGGTACAGGTAAGGACCAGATCGCCAAGCCCCGCCCACCCCCTAAAAGGGGCAGGGGCGGGACCCCCCGCCGCACCAAAACGCGCCATTTCCGCCACCCCCCCGGGGGTAACAGCCGTCAGCGCCTTCGCACCCAAACCCATAACGCCGGGCCACCCCCCCCCCACCAGCATAACGCTTTTCTCCCCgccgcccacc                                                 >  2:645411/1‑234 (MQ=255)
                                                                                                                                                                                                cagggacagtGCCAAAACGGCAGGTAGGCGAGTGGGTGGCGGTACGGGCAAGCACCAGCCCGCAAACCCCCGCCATGCCCATAAAGGTGGCAGGGGCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                               <  2:1066083/229‑1 (MQ=255)
                                                                                                                                                                                                 aGCATCATGCCAAAACGGCGGTTACGCGACCGGTTGTCGGTACAGGGAAACACCAAATCGCCCAGCCCCCCCCTGCCCCCAAAGGGGGCCGGGGCGGCACCCCGCCCCCCACCCAAACGCCACATTTTAGCCCGCCCCCCGGGGATCAGCGCCCTACCCGCCTTCCCCCCCAAACCCAAACCCGCGGGCATCCCCCCCCCCACGGCCAAAACGTTTTTcctttttctttttcacctctc                                            >  2:1186152/1‑219 (MQ=255)
                                                                                                                                                                                                             aaacCGGCGGTTACGCGACTGGTTGCCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAAt                                 <  1:637378/239‑1 (MQ=255)
                                                                                                                                                                                                                   gcggTTCCGCGACGGGTTGTCGGTACAGGGAAGCACCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGa                            <  2:557134/238‑1 (MQ=255)
                                                                                                                                                                                                                    gggTGACGCGACGGGTGGGCGGAACAGGTAACCCCCAGATCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGAt                           <  1:466317/237‑1 (MQ=255)
                                                                                                                                                                                                                    cgtttACGCGACTGGTTGTCGGTAAAGGTAAGCACCAGACCGCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGAt                           <  1:115899/235‑1 (MQ=255)
                                                                                                                                                                                                                              aCTGGGTGTCGGGACAGGGAAGCACCAGATCGCCAAGCCCCCCCATGCCCATAAAAGGGGGAGGGGCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGACAGCCCCCGGGTGATCAGCGCCGTACGCGCACTCGCACCAAAACCCAAACCGTCGGACAACCCCCCCCCCATGGCAATAACGCTTTTCACCGCCCCGCCAAGCTGCACCCCCAGTAAATCCGGATTGCTGGAAAc                 >  2:1000458/1‑238 (MQ=255)
                                                                                                                                                                                                                                       cGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtt       >  1:692224/1‑239 (MQ=255)
                                                                                                                                                                                                                                        ggTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACttt      >  1:573915/1‑239 (MQ=255)
                                                                                                                                                                                                                                        ggTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACttt      >  1:663065/1‑239 (MQ=255)
                                                                                                                                                                                                                                            cAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc  >  2:133958/1‑239 (MQ=255)
                                                                                                                                                                                                                                             aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc  >  2:937412/1‑238 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
GCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCC  >  NC_000913/3782669‑3783141

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: