Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,913 A→C 17.1% V250G (GTG→GGG)  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,9130AC17.1% 119.2 / 3.2 70V250G (GTG→GGG) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (37/19);  new base C (2/11);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782679‑3783150
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
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)
                                                                                                                                                                                                                                           cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca  >  1:176448/1‑239 (MQ=255)
                                                                                                                                                                                                                                           cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca  >  2:348616/1‑239 (MQ=255)
                                                                                                                                                                                                                                           cACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca  >  1:865866/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
AGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCG‑AAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGA‑GCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTC‑GGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA  >  NC_000913/3782679‑3783150

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: ‑